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The transmission of thermal neutrons through magnetized iron has been measured in its 
dependence upon the percentage deviation « from saturation and upon the thickness d of the 
sample. In agreement with the theory of Halpern and Holstein it was found that the percentage 
increase of transmission caused by magnetization is given by (m*p*d*/2)f(A/ed), where n is 
the number of iron atoms per unit volume. Writing for the scattering cross section of neutrons 
with parallel or antiparallel orientation of their spin with respect to the field, respectively 
ootp, we find p=2.0+0.1X10-* cm?*. From the determined value of the length A=3.2+0.3 
X< 10-* cm, the linear dimensions of the microcrystals can be determined to be = 1.4 10~* cm. 
For a thickness d=3.8 cm a transmission effect of almost 8 percent was observed if the mag- 
netization was brought to within 2.5 per mille of its saturation value; more than twice this 
effect can be expected from the same thickness at complete saturation. 


I. INTRODUCTION 


T is well known that the scattering cross 
section of a slow neutron in a magnetized 
substance depends upon the orientation of its 
spin with respect to the magnetization. This is 
due to the interaction of the magnetic moment 
of the neutron with that of the atoms in the sub- 
stance.! A rigorous theory of the phenomenon is 
complicated by various factors: It requires a 
basic assumption concerning the magnetic inter- 
action between a neutron and an electron and a 
good knowledge of the form factor for the mag- 
netic shells in the atom. It is also necessary to 
take into account the interference of the neutron 
waves, scattered from different atoms. 
For the purpose of the present paper it is 
sufficient, however, to note that the cross section 
1F. Bloch, Phys. Rev. 50, 259 (1936); 51, 994 (1937). 
J. Schwinger, Phys. Rev. 51, 544 (1937). O. Halpern and 
- H. Johnson, Phys. Rev. 55, 898 (1939). O. Halpern, 
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-- _— and M. H. Johnson, Phys. Rev. 59, 981 
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for a neutron of given velocity v with its spin 
parallel or antiparallel to the direction of mag- 
netization can always be written in the form 


(1) 


respectively, where oo stands for its value in the 
unmagnetized substance and where it is merely 
the theoretical evaluation of p in its dependence 
on v which is affected by the previously men- 
tioned complications. 

Consider now a monochromatic collimated 
beam of unpolarized neutrons with fixed primary 
intensity, which has passed through a layer of 
substance with thickness d: Let J be its intensity 
if the substance is unmagnetized, J+AJ if the 
substance is homogeneously magnetized. An ele- 
mentary consideration shows then that 


AI/I=cosh npd—1, 


ot=oot+P, 


oa =oo— p, 


(2) 


where n is the number of atoms per unit volume. 
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For npd<1 this simplifies to 
AI/I=}n*p'd’. (2a) 


Because of the smallness of p, Eq. (2a) repre- 
sents for all practical thicknesses d a sufficiently 
good approximation ; it seems desirable to verify 
experimentally the quadratic dependence of the 
“single transmission effect’? AJ/ZJ upon d and 
thus to determine the interesting quantity p. 
This, however, meets with a characteristic diffi- 
culty. All neutron polarization experiments, in- 
cluding the ones reported here, have so far been 
carried out on polycrystalline iron where the 
condition of validity of (2) and (2a), that the 
magnetization of the substance be homogeneous, 
can never be reached ideally. Because of magnetic 
anisotropy the direction of magnetization will be 
slightly different in the various microcrystals of a 
polycrystalline sample, depending upon the orien- 
tation of their crystal axis. It would in principle 
require an infinite magnetizing field to turn the 
magnetization of each microcrystal completely 
into the direction of the field and thus to reach 
perfect saturation and thereby magnetic homo- 
geneity. The first experimental evidence of a 
strong dependence of the polarization effects 
upon the degree of saturation of the iron was 
found by Powers.’ 

It was recently pointed out by Halpern and 
Holstein* that even small deviations from satura- 
tion can be understood to cause appreciable 
deviations from the ideal effect, described by (2). 
The magnetic moment of a neutron, passing 
through the polycrystal, will undergo a rapid 
precession around the direction of the magnetic 
induction B, which varies slightly from micro- 
crystal to microcrystal, and on the average this 
will result in a depolarization, counteracting the 
polarization effect described by (2). It is con- 
venient to introduce the percentage deviation 
from saturation 


e=(M,—M)/Maz, (3) 


where M is the total magnetization of the poly- 
crystal and M,, the value which it would have if 
the magnetization of all the constituent single 
crystals were perfectly parallel to the magnetizing 
field. The whole transition from complete de- 


? P. Powers, Phys. Rev. 54, 827 (1938). 
30. Halpern and T. Holstein, Phys. Rev. 59, 960 (1941). 
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polarization to full polarization will be found to 
occur in a state of almost complete saturation, 
and it is therefore sufficient to consider only the 
situation where e<1. In the limiting case «1 
and npd<1, which is the only one of practical 
importance, the theory of Halpern and Holstein 
yields instead of (2a) the formula 


AI/I=}-n*p*d?: f(d/ed), (4) 


where X is a length, correlated to the linear di- 
mensions of the single microcrystals and where 


1 
fls)=23*(e-4!4-—1). (5) 
x 


f is a monotonic function of its argument with 
the property f(0)=0, f(«©)=1; it represents the 
reduction factor of the single transmission effect 
due to depolarization and reaches unity only 
asymptotically as ¢ vanishes, i.e., according to 
(3), as complete saturation is reached. 

In the relation between \ and the linear dimen- 
sions of the microcrystals there enters the dis- 
tance which a neutron trayels during one full 
Larmor revolution. Let / be that distance, divided 
by 2z, i.e., the length of path of a neutron while 
its magnetic moment undergoes a precession by 
one radian. It is given by 


l=v/Bg, (6) 


where g is the ratio of the magnetic moment of 
the neutron to its angular momentum, B the 
magnetic induction inside a microcrystal, and 
the velocity of the neutron. Assuming v to be 
the thermal velocity at room temperature, i.e., 
v=2.5X10° cm/sec., taking B= 20,000 gauss and 
the measured‘ value g=1.86X10** c.g.s. one 
obtains from (6) 


1=6.7X10~ cm. (7) 


According to Halpern and Holstein one ob- 
tains a simple relation between \ and a length 6 
of the order of the linear dimensions of the micro- 
crystals in the two limiting cases where 6 is 
small or large compared to /. One has 


A\=P/5, (8) 
A= 6/2, (8a) 


for 6<l: 
for 6>1: 


4L. Alvarez and F. Bloch, Phys. Rev. 57, 111 (1940). 
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and for 6=1: AS6. (8b)® 


It is evident, from what has been said before, 
that a quantitative study of the single trans- 
mission effect requires not only the measurement 
of the percentage change of neutron transmission 
AI/I upon magnetization of the sample, but an 
equally accurate measurement of the percentage 
deviation from saturation, reached in the mag- 
netized state, i.e., of the quantity e, entering in 
Eq. (4). In the work presented here, we have 
determined AJ/I for a series of values of ¢ and 
for two different thicknesses d of the same sample 
of hot rolled steel with the threefold objective : 


5 It is clear that Eqs. (8) and (8a) are somewhat qualita- 
tive, insofar as 6, although doubtlessly of the order of the 
linear dimensions of the microcrystals, actually has to be 
obtained as a very complicated average value, not only 
over the sizes, but also over the varying geometrical shapes 
of the microcrystals; furthermore, as pointed out by 
Halpern and Holstein, there enters in the rigorous calcu- 
lation the demagnetizing field from the surfaces of the 
microcrystals. It may be noted that the rigorous and 
general expression of the length A entering in (4) can 
actually be written down as the average value of a certain 
double integral over the volume of a single crystal; 
assuming random orientation of the microcrystals it then 
appears that, due to the demagnetizing field, \ depends 
also upon the angle @ between the direction of the traversing 
neutrons and that of the magnetization. In the case 
corresponding to (8) one obtains . 


A~(1+ } sin* 6). 


In the case corresponding to (8a), the length A is simply 
correlated to the total area of the surfaces of all micro- 
crystals. In fact, if we denote by a the mean value of this 
area per unit volume of the sample, one obtains in this case 
1 cost @\~" 
cman — 7 ede ‘ 
r gq ( 9+2 cos? @ 3 ) 


In view of the fact, however, that in our arrangement we 
had always 6=x/2 and that the shape and size of the 
microcrystals are not known the simple formulae (8) and 
(8a) are quite sufficient for our purposes. 


POLARIZATION 


Fic. 1. Arrangement of apparatus. 


























1. To verify the depolarization theory of 
Halpern and Holstein. 

2. To determine the ‘‘polarization cross sec- 
tion” p, defined by (1). 

3. To obtain information about the charac- 
teristic length A, entering in (4) and thereby 
through (8) about the linear dimension 6 of the 
microcrystals in the sample under consideration. 


II. MEASUREMENT OF THE SINGLE 
TRANSMISSION EFFECT 


The method and arrangement to determine 
the single transmission effect were essentially the 
same as used by Alvarez and Bloch‘ in their 
determination of the neutron moment. The slow 
neutrons (see Fig. 1) emerged from a paraffin 
howitzer directly in front of the target chamber 
of the Stanford cyclotron® where a Be target was 
bombarded with 2.5-Mev deuterons. In a double 
screening wall of water-filled tin cans for the 
absorption of neutrons and of bricks for the 
absorption of y-rays, there was in front of the 
howitzer a circular hole of 2-in. diameter to 
obtain a collimated neutron beam. This beam 
passed through a square plate of hot rolled steel 
between the poles of an electromagnet at right 
angles to the magnetizing field and from there 
through a Cd tube of 2-in. diameter and 26}-in. 
length which provided further collimation to the 
detecting BF; filled ionization chamber. In order 

* The construction of this 27}-in. cyclotron by Bloch, 
Bradbury, Staub, and Stephens was made possible by 
grants from the Rockefeller Foundation, the Research 
Corporation, and from various private sources. It was 


completed in the autumn of 1941, giving under normal 
operation a beam of about 10 wamp. of 2.5-Mev deuterons. 
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to prevent neutrons, scattered on the material 
of the magnet from entering the chamber, Cd 
diaphragms were placed across the neutron beam 
both on the front side and on the back side of 
the magnet. Except for a circular hole of 2-in. 
diameter in front, the ionization chamber was 
completely shielded with Cd and embedded in 
a large block of paraffin, surrounded by a layer 
of water-filled tin cans, leaving only its front 
side open. All these precautions were necessary 
in order to obtain conditions where the majority 
of the recorded neutrons were of thermal energy 
and had actually passed through the strongly 
magnetized region of the steel sample ; even with 
a thickness of steel traversed of as much as 3.8 
cm we found that 67 percent of the neutrons 
recorded were thermal and had passed through 
the steel. This was ascertained by shielding just 
the steel sample with Cd and counting as back- 
ground the remaining neutrons, partly coming 
from elsewhere and partly consisting of those 
collimated neutrons above the thermal region 
which were still able to pass the Cd shield. 

The boron trifluoride ionization chamber which 
was similar to the one used by Powers? had an 
active volume of 2-in. diameter and 6-in. length. 
It was filled with 2 atmos. of BF;. The aluminum 
electrodes consisted of three concentric cylinders. 
The most inner and outer electrodes were con- 
nected to a potential of 3500 v and the inter- 
mediate one to the grid of a 38 tube. Porcelain 
was used as an insulator and during a full year 
did not seem to cause any trouble. Both leads 
were brought out of the chamber by spark plugs. 
The bottom of the chamber was sealed by an 
ordinary rubber gasket. 

The boron-disintegration pulses in the ioniza- 
tion chamber were amplified in a Wynn-Williams 
type amplifier and, through a scale of 8 or 16 
_ scaling circuit, were recorded by a Cenco counter. 
Since the collection time for the ions is by no 
means very small compared to the characteristic 
time of the amplifier, the number of pulses 
recorded would strongly depend on the voltage 
applied to the chamber. This effect can be re- 
duced if one sets the bias of the recording device 
so that it is below the size of the great majority 
of the pulses which, due to the imperfect collec- 
tion, consist of a,rather broad group tailing off 
toward lower energies. The runs were taken with 


good steady conditions of the cyclotron, and we 
always operated the beam at such an intensity as 
to give us a rate of about one count per second 
on the Cenco counter. In order to be independent 
of fluctuations of the neutron intensity, we used 
the monitor of Fryer and Staub,’ consisting of 
another small BF; ionization chamber and an 
electronic circuit which periodically recorded in- 
tervals of an equal amount of charge having been 
collected in this chamber. The monitor chamber 
was surrounded with paraffin and placed at a 
distance of about one foot from the Be target of 
the cyclotron. If the monitor would be ideal and 
no changes in the recording apparatus would 
occur, one should obviously record the same 
number of neutrons per monitor interval, inde- 
pendent of the intensity of the cyclotron beam. 
Actually, we found in the course of our measure- 
ments that this number fluctuated by about 1 
percent; this may have been due to a slight de- 
pendence of the monitor performance upon the 
position of the target spot or to variations of the 
line voltage and, therefore, of the collecting 
voltage of the ionization chamber (which would 
affect the rate of recorded pulses). In order to 
eliminate the influence of these slow variations 
upon the accuracy of our measurements, the 
magnetizing current was turned on and off in 
alternate monitor intervals, and the percentage 
difference of the total number of counts, ob- 
tained under the two different conditions was 
then taken as a measure of the transmission 
effect. 

By a suitable choice of the condenser, dis- 
charged by the ionization current in the monitor 
chamber, the monitor intervals were made of an 
approximate duration of 14 min. At the end of 
each monitor interval conditions were auto- 
matically changed through a telephone selector 
relay, which at the same time switched alter- 
natingly over from one to the other of the two 
Cenco counters, thus allowing one .counter to 
read all counts taken with and the other counter 
all those without magnetizing current. Measure- 
ments were performed on two samples of hot 
rolled steel of identical material. They were 
made of machined square plates of 5X5 cm, one 
sample having a total thickness of 1.91 cm, the 


7 E. Fryer and H. Staub, Rev. Sci. Inst. 13, 187 (1942). 
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other of 3.81 cm. Two opposite sides of the square 
were clamped between the poles of the magnet, 
its face perpendicular to the neutron beam and 
carefully centered with respect to it. The runs 
were taken by setting the magnetizing current 
to a certain value and then taking counts with 
this current alternatingly on and off; the setting 
of the current was frequently changed to ascer- 
tain that independent runs with the same value 
of the current gave statistically consistent results. 
A sufficiently large number of counts was taken 
for each setting to obtain the difference of the 
neutron intensity with and without magnetizing 
current to within an accuracy of about } percent 





and 3} percent of the intensity for the thin and 
the thick sample, respectively. The following 
tabulation gives the results for the two samples 
and for various values of the magnetizing current 
up to the maximum value of 15 amperes which 
we were allowed to pass permanently through our 
magnet. The transmission effect AJ/I and its 
probable error tabulated here were obtained by 
dividing the directly observed difference of in- 
tensity by I(1—b), where 6 is the percentage 
background, taking into account that the back- 
ground is independent of the presence or absence 
of the magnetizing current. 


Thin sample (1.91 cm). Percentage background b=6.95+0.1 percent. 


0.50 
0.43.13 


Current (amperes) 
AI/I (percent) 


AI/I (percent) 3.33 4.25 


0.70 
1.76+.25 


Current (amperes) 5.0 9.0 
3.27+4.25 


1.0 2.0 3.0 
1.08+.19 1.65+.24 2.78+.26 
15.0 
3.67.25 


Thick sample (3.81 cm). Percentage background b = 32.8+0.81 percent. 


Current (amperes) 1.0 1.2 
—0.06+.72 
Current (amperes) 6.9 9.0 
6.39+.71 


AlI/I (percent) 0.80+.76 


AI/I (percent) 4.51+.69 


It is of course important, when dealing with 
such small effects, that any small source of error 
be eliminated. Both before and between our runs 
various tests and blank runs were undertaken to 
ascertain carefully the reality of the observed 
effects. The test probably most severe consisted 
of runs where the iron sample was substituted 
by a brass sample with a thickness so as to trans- 
mit the same number of neutrons as the iron and 
under otherwise identical geometrical and physi- 
cal conditions; the fact that we thus obtained 
well within the accuracy of our transmission 
effects no difference in the recorded counts, 
taken with and without the magnetizing current, 
shows particularly that the stray field of the 
magnet does not disturb the recording apparatus 
since it would be even stronger with brass than 
with iron between the poles. 


Ill. MEASUREMENT OF THE MAGNETIZATION 


As pointed out in Section I, it is not the de- 
pendence of the single transmission effect upon 
the magnetizing current, which is of physical 
interest, but its dependence upon the magnetiza- 
tion of the sample. After carrying out the meas- 
urements, described in the previous paragraph. 


2.0 3.0 4.0 
2.79+.69 3.24+.76 3.27.77 
11.25 15.0 
7.22+.75 7.70+.69 





we thus had to devise methods by which the 
small deviations from saturation reached in our 
experiments could be determined with fairly 
good accuracy. 

We define the magnetization by 


M=(B—H)/4r, (0) 


whereby its determination rests upon that of the 
difference of the magnetic induction B and the 
magnetic field H in a region where one can expect 
both to be uniform. The geometrical shape of the 
square plates used in the transmission measure- 
ments is rather inconvenient for such a determi- 
nation ; besides, since their linear dimensions are 
comparable to those of. the poles of the magnet, 
B and H cannot be expected to be very uniform 
all over the sample. We had therefore first to 
determine ¢ in its dependence upon H on a sample 
of the identical material, but of such size and 
shape as to give more favorable conditions for 
the measurement of the magnetization. Such 
conditions were found in the so-called “isthmus 
method”’ which we have used, where a long and 
thin cylindrical piece of the material under test, 
the “isthmus,” is mounted across the poles of 
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the magnet and where B and H are measured in 
its central region. The arrangement is schemati- 
cally indicated in Fig. 2. The cylindrical isthmus 
is wound with a winding w:, consisting of a 
certain number of turns. On a coaxial brass 
casing is wound a second winding we, consisting 
of the same number of turns. If we assume B 
and H to be parallel to the cylinder axis within 
the cylindrical sections, enclosed by w; and we, 
their changes AB and AH upon a sudden change 
of the magnetizing current can be determined 
from the ballistic throw of a galvanometer G, 
connected in series with the windings. Indeed, 
let A be the wound area of the isthmus, a; and a 
the wound areas of the air space surrounded by 
the windings w; and we, respectively ; let, further, 
Dy, be the ballistic deflection of the galvanometer 
if connected in series with w; and D, the deflec- 
tion if G, w;, and we are all three connected in 
series, the latter two such that their induced 
electromotive forces oppose each other. Taking 
the galvanometer readings upon variation of the 
magnetizing current between the same two limits, 
one has then 


Dz 
~ k(ag—A—ay) 





(10) 


1 a, 
AB=— (v.-——».). (10a) 
kA a2—A—a, 


The constant k depends upon the characteristics 
of the galvanometer and other elements of the 
circuit; we determined it in the usual way by 
calibration with an air coil. Its actual value, 
however, does not enter in the ultimate determi- 
nation of e. 
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Fic. 2. Arrangement of windings and sample for 
the determination of H and B 
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The requirement of accuracy in the determi- 
nation of the quantities (10) and (10a) is evi- 
dently quite stringent since their difference 
AB—AH, which determines the change of mag- 
netization, becomes increasingly smaller upon 
approach to saturation. Changing, for example, 
the magnetic H field in the neighborhood of 3000 
gauss by 300 gauss, we found the corresponding 
change of B—H to be about 10 gauss, or } per 
mille of its saturation value. In order to ascertain 
this change to within 5 gauss or } per mille of 
saturation, it requires therefore a determination 
of the deflections D,; and D: and of the wound 
areas A and dz such that in-the difference of 
(10) and (10a) they do not involve an error of 
more than 5/300 or about 1} percent, i.e., each 
of these quantities must be known with an 
accuracy of about 1 percent; a; which itself 
represents only a small correction, due to the 
finite thickness of wire and insulation of winding 
wi, has to be known with correspondingly less 
accuracy. While it is true that the determination 
of these quantities with the required accuracy 
demands great care, it is not prohibitive, and 
we feel confident that we have succeeded. The 
determination of the wound areas requires merely 
well-machined cylindrical suffaces, the diameter 
of which can be easily measured with a microm- 
eter ; the correction for wire thickness and insula- 
tion must, of course, be applied after remeasuring 
the diameter when the wire has been wound on. 
A sufficient accuracy of the galvanometer reading 
can be ascertained by patient repetition. 

The only remaining source of error can arise 
if B and H are appreciably non-uniform in the 
investigated central region of the isthmus. While 
it is true that in the absence of the isthmus the 
magnetic field in this region is sufficiently uni- 
form and that its distortions by the isthmus be- 
come negligible if the isthmus is thin enough, 
one has nevertheless to make sure that this uni- 
formity is maintained under the conditions of the 
experiment. We have considered as criterion for 
sufficient thinness of the isthmus that making 
the isthmus still thinner does not change the 
result. After carrying out the measurements with 
an isthmus of 5.1-cm length when its own diam- 
eter was 0.961 cm and that of the brass casing 
1.430 cm, we have repeated them once more with 
another isthmus of the same length and the two 
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Fic. 3. Change of B—H versus 1/H*. The actual value of B—H is obtained by 
adding to the value of A(B—H) the constant value 19,800 gauss. 


diameters, respectively 0.638 cm and 0.946 cm. 
Since the results, obtained for the approach to 
saturation, agreed in the two cases within the 
experimental error, it seems safe to assume that 
the error involved from non-uniformity was 
negligible. > 

There is another independent and rather severe 
test for the correct determination of the wound 
areas and for the uniformity of H and B. If the 
magnetizing current is more and more increased, 
and the corresponding increases of B and H are 
computed from (10) and (10a), it must turn out 
that while both B and H keep on increasing, 
their difference must converge towards a finite 
value, which through (9) determines the satura- 
tion value of the magnetization M. Both non- 
uniformities and errors in the wound areas would 
manifest themselves in an unlimited increase or 
decrease of the computed value of M. Unless 
there is an accidental cancellation of these errors, 
the very fact of the computed value of M con- 
verging towards a finite saturation value indi- 
cates the correctness of the performed meas- 
urement. 

From the theoretical side one can predict the 
law of approach to saturation as long as one 
assumes that in the magnetization process of a 
macroscopic sample the magnetic moment of 


each microcrystal is changed only in its orienta- 
tion but not in its magnitude. One obtains then 
for strong fields 


M=M,-—a/H?, (11) 


where a depends upon the material* and is a 
slowly varying function of H. 

This law has been verified with great accuracy 
by Czerlinsky ;* the small deviations which he 
has observed for very strong magnetic fields in- 
dicate indeed a slow increase of a with H. Within 
the experimental error we were able for high 
magnetic fields to represent our results by a law 
of the form (11) with a constant a, indicating 
that for the material and the fields used by us, 
the above-mentioned deviations were of no im- 
portance. In evaluating the percentage deviation 
e of Eq. (3), we have used the deviation of the 
magnetization from that value, obtained by (11), 
upon extrapolation to H= © with constant a."® 


ismus (Springer, 


® R. Becker and W. Doring, apie seme ly 3 
. Primakoff, Phys. 


Berlin, 1939), p. 167. T. Holstein and 
Rev. 59, 388 (1941). 

*E. Czerlinsky, Ann. d. Physik 13, 89 (1932). 

1° Tt should perhaps be mentioned that it is essential to 
use (11) for extrapolation; the increase of magnitude of 
the microcrystalline magnetization, which will actually 
occur in very strong fields and which is omitted in that 
formula, has no effect upon the depolarization and on 
the value of ¢, to be used in the sense of the theory. Our 
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In Fig. 3 we have plotted A(B—H) and its per- 
centage deviation from saturation ¢ in their de- 
pendence on the magnetic field H, as obtained 
by the isthmus method for hot rolled steel. Al- 
though deviations of the order of one-tenth of 
a percent are to be considered significant in this 
plot, we do not claim nearly such an accuracy 
for the absolute determination of the extrapolated 
value of the magnetic saturation. 

Once we have established the relationship be- 











TABLE I. Percentage deviation (e) from saturation. 
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Thin sample (1.91 cm) Thick sample (3.81 cm) 
Current H € Current H € 
(amp.) (gauss) (percent) | (amp.) (gauss) (percent) 

0.50 592 2.50 1.0 570 2.56 
0.70 795 1.70 1.2 678 2.01 
1.0 1100 1:16 2.0 1060 1.21 
2.0 1690 0.63 3.0 1360 0.87 
3.0 2140 0.41 4.0 1590 0.69 
5.0 2710 0.25 6.9 2060 0.44 
9.0 3370 0.17 9.0 2280 0.38 
15.0 3980 0.12 11.25 2480 0.30 
15.0 2770 0.24 





















tween ¢€ and H by the isthmus method, the 
problem of assigning the proper value of « to a 
given magnetization current in the single trans- 
mission experiments, described in Section II, 
becomes comparatively simple. It reduces to the 
determination of the magnetizing field H for a 
given value of the current, the percentage error 
in ¢ thus involved being merely of the order of 
that of H. To determine the field in the central 
region of the plate, through which the neutrons 
pass in the transmission experiment, we used a 
thin rectangular search coil, placed right along- 
side the face of the plate. Its wound area was 
again effectively determined by calibrating the 
ballistic throw of the galvanometer against that 
produced by the known field at the center of a 
current solenoid. The position of the search coil 
was not very critical since, moving it over the 
whole area through which neutron passage occurs, 
we did not find the field to vary by more than 
about 3 percent; with the same accuracy we 
found the field along the entrance and exit face 
of the plate to be the same. In addition, we had 
to convince ourselves that between these two 























extrapolation with constant a is very nearly correct, since 
variations would make themselves noticeable only in con- 
siderably higher fields. 
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faces there were no appreciable deviations of H; 
since the thicker sample consisted actually of 
three, the thinner of two separate plates, this 
could be done by spacing these plates just enough 
to insert other rectangular search coils between 
them. It was thus easy to determine any possible 
difference of H between the plates and along their 
outside faces. Having found no such difference 
within the accuracy of our measurements of a 
few percent, we see no reason to doubt that the 
field H was sufficiently uniform throughout the 
passage region of the neutrons, and that for a 
given magnetizing current we had that percent- 
age deviation from saturation which from the 
previous isthmus measurements corresponded to 
the field measured for this current. 

Table I gives the values for H and the corre- 
sponding ¢ for those currents for which the trans- 
mission effect was measured. 

While measuring the magnetizing field for 
these samples, we measured at the same time 
the induction by winding the plates in the central 
region and observing the induced current pulse in 
the galvanometer. The values for the magnetiza- 
tion thus obtained agreed with those of the 
isthmus measurements for percentage deviations 
from saturation above about 0.5 percent. For 


TABLE II. Dependence of single transmission 
effect AI/I upon e. 














Thin sample (1.91 cm) Thick sample (3.81 cm) 
e (percent) Al/I (percent) « (percent) Al/I (percent) 
2.50 0.43+.13 2.56 0.80+.76 
1.70 1.76+.25 2.01 —0.06+.72 
1.16 1.08+.19 1.21 2.79+.69 
0.63 1.65+.24 0.87 3.24+.76 
0.41 2.78+.26 0.69 3.27.77 
0.25 3.334.25 0.44 4.51+.69 
0.17 3.274.25 | 0.38 6.39+.71 
0.12 3.67 +.25 0.30 7.22+.75 
0.24 7.70+.69 
| 








smaller deviations we found systematic and ever- 
increasing differences, evidently arising from the 
fact that under the unfavorable geometrical con- 
ditions, presented by the square plates, the 
knowledge of the wound areas and the uniformity 
of B and H became insufficient to make the 
determination for the smaller deviations reliable. 
The agreement was, however, sufficient to indi- 
cate that there were no accidental differences in 
the material used for the isthmus and the square 
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Fic. 4. Single transmission effect for both plate thicknesses as a function of 
deviation from saturation. Heavy line represents best fitting theoretical curve 


with corresponding saturation effect E.. 


plates, and that we could indeed expect the 
same relation between e and H for both. 


IV. DISCUSSION OF THE RESULTS 


From the measurements, reported in Section II 
and Section III, we are now able to give the final 
results in Table II, which represents the de- 
pendence of the single transmission effect AJ/I 
upon the percentage deviation from saturation ¢ 
for our two samples. 

These results are plotted in Fig. 4; the ordi- 
nates are the single transmission effects E= AJ/IJ, 
the abscissae the corresponding values of 1/e. 
Since the ratio of the thicknesses d of the two 
samples is just about 1 : 2, we have chosen a 
representation in which for the thick sample the 
scale of the abscissae is reduced by a factor 2, 
that of the ordinates by a factor 4. 

One can see that by this transformation the 
points obtained for the two samples fall within 
the statistical error into the same pattern; this 
shows that in agreement with Eq. (4) the polar- 
ization effect is proportional to the square of 
the thickness and that the depolarization factor 
f depends merely upon the product ed. Quite 
independent of the special form (5) of the func- 
tion f, we consider this agreement as a first im- 


portant test for the correctness of the theory of 
Halpern and Holstein. 

The complete quantitative test of this theory 
is evidenced by the drawn out curve of Fig. 4. 
It represents the best fitting member of the two- 
parameter family of curves which one obtains 
from (4) and (5) by varying the two unknown 
constants p and X. The thus experimentally de- 
termined values of these two constants are with 
their estimated errors 


p=2.0+.1XK10-™ cm’, 
A=3.2+.3X10-* cm. 


(12) 
(13) 


The best theoretical evaluation of p is given by 
Hamermesh ;" actually he finds for the quantity 
w=np (n=the number of iron atoms per cc) 
w=0.090 cm from which it would follow with 
n=8.2X10"/cc 


p=1.1X10™ cm’; (14) 


this is smaller by a factor 1.8 than the observed 
value, and therefore from (4) would lead to 
about three times too small transmission effects. 
The results of Powers? for thin iron agree with 
our experimental determination of p, indicating 


11M. Hamermesh, Phys. Rev. 61. 17 (1942). 














that his iron was sufficiently magnetized to 
reach the saturation value. 

In comparing the theoretical and experimental 
values of p, it has to be remembered that the 
calculation of the former is based upon the not- 
too-well-known wave functions for the 3D-shell 
of Fe. A qualitative description of the disagree- 
ment in terms of the wave functions can be given 
by stating that all their linear dimensions would 
have to be contracted by a factor 13 in order to 
reestablish agreement. It is true that so big an 
error in the scale of wave functions, calculated by 
the Hartree method, is surprising and that it is 
by no means certain that this is the cause of the 
observed discrepancy. On the other hand, the 
neutron polarization depends on the atomic form 
factor of the 3D-shell alone and represents a 
rather more severe experimental test on the 
Hartree method than the spectroscopic and x-ray 
measurements by which it has been tested 
before. 

As another alternative, Halpern and Johnson” 
have suggested the possibility of the neutron 
having the spin 3, instead of the usually accepted 
value 3. In view of the difficulties which such an 
hypothesis would introduce in the interpretation 
of nuclear phenomena, the sole support from the 
above discrepancy in the single transmission of 
neutrons does not seem to us sufficient to make 
it appear probable. 


%0O. Halpern and M. H. Johnson, Phys. Rev. 57, 160 
(1940). 
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As mentioned in Section I, the determination 
(13) of \ is connected with the linear dimension § 
of the microcrystals. According to (8) and (8a) 
two values of 6 are possible for a given X, de- 
pending upon whether 6 is small or large com- 
pared to the length /, given in (7). Depending 
upon the choice we get 


from (8) 6=1.4X10-' cm, (15) 
from (8a) 56=6.4X10-* cm. (15a) 


The value (15a) is far too big to be plausible, 
and we consider the value (15) as a fair determi- 
nation of the microcrystalline dimensions by 
means of this method. 

In view of our results we foresee no essential 
experimental difficulties to overcome the rather 
strong depolarizing influences and thus to obtain 
considerably higher polarization effects than the 
ones obtained so far. It is true that in the case 
of thin iron we were able to obtain almost the 
complete effect, but for obvious reasons we were 
far from doing so with the thick sample. This 
could be helped, no doubt, by raising the mag- 
netizing field to the order of 10,000 gauss which 
would merely require a somewhat larger electro- 
magnet than the one at our disposal. Besides, 
however, it may well turn out that hot rolled 
steel is by no means the ideal material in the 
sense that other grades of iron will consist of 
smaller microcrystals or may be found easier to 
saturate, both desirable features to keep the 
necessary magnetizing field as low as possible. 
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The nuclear spin of the chlorine atom was determined by 
observing the line intensities in the band spectrum of the 
homonuclear diatomic molecule of that atom. The lines of 
even J values show a different intensity from those of 
odd J values, the ratio of intensities being given by 
(1+1)/I where J is the nuclear spin of the atom in question. 
In the band absorption spectrum of chlorine the intensity 
ratio of lines of the Cl*®CI* molecule has been measured 
and compared with the value found by Elliott. These 
values are in substantial agreement and yield a value of 
5/2 for the nuclear spin of Cl*. In this spectrum no lines 
of the molecule Cl°7CI*7 were observed due to the low con- 


1. INTRODUCTION 


T is a well-known fact that homonuclear dia- 
tomic molecules and certain linear polyatomic 
molecules exhibit the phenomenon of alternating 
intensities in their band spectra. This intensity 
ratio arises from the fact that the statistical 
weights of rotational levels with (+) symmetry 
are different from the statistical weights of levels 
with (—) symmetry. This ratio can be shown 

to be 

(I+1)/I 


for homonuclear diatomic moleculeswhere J is the 
nuclear spin of the nucleus of the element in 
the molecule. For polyatomic molecules the in- 
tensity ratio is more complicated. The nuclear 
spin of a nucleus may therefore be determined by 
observing and measuring the alternating intensity 
ratio in the band spectrum of its homonuclear 
diatomic molecule. Elliott' analyzed the ':—*0 
band system of Cl, in absorption and made a 
determination of the alternating intensity ratio 
of the molecule Cl®C1l*. The value of this ratio 
was 1.35 from which he concluded that the most 
probable value of the nuclear spin of Cl®* is 5/2. 
In several bands he was able to identify the lines 
of the molecule Cl*®CI*? but due to the small 
concentration of Cl*? (24.2 percent) he was not 
able to identify lines of the homonuclear mole- 


ons Now at Case School of Applied Science, Cleveland, 
io. 
(19a) Elliott, Proc. Roy. Soc. 123, 629 (1929); 127, 638 
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centration of Cl*’ in the gas. By the method of thermal 
diffusion 200 cm’, S.T.P., of chlorine gas was obtained 
containing 45 percent Cl*?. The absorption spectrum of 
this gas showed lines of the Cl*7Cl*? molecule in the 112 
band of the '£-+*0 band system of Cly. For this band the 
alternating intensity ratio for the Cl*7Cl*7 molecule (1.28) 
was the same, within the experimental error, as for the 
Cl*CI® molecule (1.26), indicating that the nuclear spin 
of Cl*? is probably also 5/2. In interpreting the intensity 
ratios one must take into account the overlapping of lines 
and general background present. 






cule Cl’Cl*’. Elliott found no evidence of alter- 
nating intensities in the spectrum of CI®Cl*’. 
This agrees with theoretical predictions. 

It is now possible, having obtained chlorine 
enriched in the heavy isotope, to obtain the band 
spectrum of Cl*"Cl* and make a determination 
of the nuclear spin of Cl*. 


2. EXPERIMENTAL 


The spectrum of Cl, was obtained in absorp- 
tion using a 21-foot concave reflection grating 
ruled 30,000 lines to the inch. An Eagle mounting 
was used which gave a dispersion of about 1.25A 
per millimeter in the first order. 

Only a limited amount, approximately 200 cm’ 
S.T.P., of “enriched” chlorine containing 45 per- 
cent Cl*’ was available. For this reason an ab- 
sorption tube was constructed to give the greatest 
absorption for this small volume of gas. The 
optical arrangement is shown in Fig. 1. The 
absorption tube was 2.5 meters long and 1.2 cm 
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Fic. 1. Optical arrangement of apparatus. 
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inside diameter. After cleaning thoroughly with 
nitric acid and washing, windows were sealed on 
the previously ground ends of the tube. Chlorine 
was admitted to the evacuated tube. The pres- 
sure in the tube was measured by means of a 
mercury manometer protected by a thin layer of 
concentrated sulphuric acid. The tube was filled 
with Cl, to a pressure of 50 cm of mercury. 

Light from a 100-candlepower Point-O-Lite 
lamp was focused on the small mirror E which 
reflected the light into the tube. A second 
aluminized mirror at the far end reflected the 
light back down the tube past mirror E, through 
cylindrical lens G, and onto the slit of the grating 
spectrograph. The plane glass window at the 
grating end of the tube was placed at a slight 
angle to prevent light from being directly re- 
flected into the slit. This arrangement necessarily 
sacrifices light intensity in order to obtain suffi- 
cient absorption. An exposure time of 48 hours 
was required using a slit width of 0.03 mm. 

The standard photometric technique using a 
step filter for plate calibration was employed. 
The step filter was calibrated by means of a 
photronic cell and a quartz monochromator. 
A Koch and Goos photographically recording 
microphotometer was used to reduce the spec- 
tographic plate into an intensity curve. Each 
line was resolved graphically from its neighbors 
and its intensity calculated from plate calibration. 





TABLE I. Wave-lengths, etc., of band lines in Cl*CI*. 
Band 1-+12. P and R branches. 








Ain A vin cm~ 








15 17 5164.77 19356.57 0.019 
16 18 5165.86 19352.47 0.014 
17 19 5167.03 19348.11 0.017 
18 20 5168.26 19343.50 

19 21 5169.53 19338.74 0.016 
20 22 5170.88 19333.69 0.019 
21 23 $172.29 19328.42 0.020 
22 24 5173.76 19322.93 0.013 
23 25 5175.33 19317.06 0.015 
24 26 5176.92 19311.13 0.013 
25 27 5178.62 19304.80 0.014 
26 28 5180.39 19298.22 0.013 
27 29 5182.19 19291.49 0.014 
28 30 5184.08 19284.46 

29 31 5186.03 19277.22 0.017 
30 32 5188.04 19269.76 0.009 
31 33 5190.14 19361.97 0.016 
32 34 5192.21 19253.96 0.012 
33 35 5194.51 19245.73 0.017 
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3. RESULTS AND DISCUSSION 
a. Analysis of Spectrum of “Enriched’”’ Chlorine 


The strongest discrete absorption of chlorine is 
the band progression 1-11, 1-12, and 1-13, 
which lies between 5130A and 5230A. The elec- 
tronic transition is '2—>+0 and consequently each 
band consists of single P and R branches only. 


TABLE II. Wave-lengths, etc., of band lines in CCI’. 
Band 1-12. P and R branches. 











Jp Jr Ain A vin cm™ a(max.) 
11 13 5163.66 19360.72 0.020 
12 14 5164.40 19357.97 0.020 
13 15 5165.31 19354.55 0.020 
14 16 5166.25 19351.03 0.024 
15 17 5167.25 19347.29 0.026 
16 18 5168.27 19343.47 

17 19 5169.40 19339.24 - 0.029 
18 20 5170.58 19334.83 0.023 
19 21 5171.86 19330.03 0.032 
20 22 5173.14 19325.23 0.026 
21 23 5174.53 19320.04 0.019 
22 24 5175.94 19314.82 0.024 
23 25 5177.44 19309.20 0.017 
24 26 5179.02 19303.33 0.016 
25 27 5180.64 19297.26 0.016 
26 28 5182.38 19290.81 0.015 
27 29 5184.09 19284.44 

28 30 5185.92 19277.62 

29 31 5187.78 19270.73 0.017 
30 32 5189.77 19263.34 0.026 
31 33 5191.79 19255.82 0.022 
32 34 5193.88 19248.09 0.026 











TABLE III. Wave-lengths, etc., of band lines in Cl*7Cl*". 
Band 1-—+12. P and R branches. 








Jp Jr Ain A vin cm a(max.) 














7 9 5163.85 19360.00 0.010 

8 10 5164.25 19358.52 0.012 

9 11 5164.91 19356.01 0.015 

10 12 5165.57 19353.57 0.013 
5166.41 19350.43 





13 15 5167.83 19345.12 0.013 
14 16 5168.76 19341.25 0.011 
15 17 5169.65 19338.39 0.016 
16 18 5170.73 19334.26 0.009 
17 19 5171.68 19330.71 0.0012 
18 20 $172.79 19326.48 0.009 
19 21 5174.12 19321.58 0.010 
20 22 5175.33 19377.06 

21 23 5176.77 19311.68 0.013 
22 24 5177.89 19307.44 0.008 
23 25 $179.52 19301.46 0.010 
24 26 5181.00 19296.00 0.005 
25 27 5182.67 19289.73 0.008 
26 28 5184.27 19283.75 0.005 
27 29 

28 30 5187.58 19271.43 

29 31 5189.53 19264.20 0.007 
30 32 5191.35 19357.47 0.008 
31 33 5193.38 19249.93 0.010 
32 34 5195.26 19242.98 0.008 
























NUCLEAR SPIN OF Cl? 


In band 1-12 the P and R branches fall one on 
top of the other, the lines of one branch lying 
on the lines of the other whose J values are two 
removed. Hence even J values fall on even and 
odd on odd and the alternating intensity phe- 
nomenon is retained, not obliterated, as would be 
the case if even J value fell on odd J value. 
From a spectrogram taken with the enriched 
chlorine it has been possible to identify 24 lines 
which belong to the molecule Cl*’CI*’ in the 


TABLE IV. Rotation constants of band 1-+12 of Cl*’Cl’. 








B” B 
0.1115 cm" 0.2228 cm—! 











1—12 band. Due to the increased complexity and 
over-lapping of lines it has been impossible to 
identify any lines in the 1-11 and 1-+13 bands 
as belonging to the molecule Cl*7Cl*’. 

Tables I, II, III give wave-lengths, etc., for 
the lines of the three molecules. The rotational 
constants are given in Table IV. The general 
appearance of the spectra of normal and enriched 
chlorine may be seen in Fig. 2. 

The maximum value of the absorption coeffi- 
cient, a(max.), for a given line is calculated from 
the expression 

T= Toe *??, 


where J is the intensity at the center of the 
absorption line, p is the chlorine concentration 
in gram molecules per liter, and x is the thickness 
of chlorine in cm. Elliott? studied the relations of 
slit widths and line shapes to the line intensities 
and found no appreciable error was introduced 
in neglecting their effect. 


b. Alternating Intensity Measurements 


Elliott’s measurements on the alternating in- 
tensity ratio of Cl*CI* in band 1-412 of ‘“‘normal”’ 
chlorine were repeated by the author. The value 
obtained was 1.31 as compared to Elliott's 
value of 1.34 for this band. The relative intensity 
of lines of the Cl*CI* molecule to lines of the 
ClI*Cl*? molecule was found to be 1.32 as against 
Elliott’s value of 1.35 for the same ratio. These 
values are in agreement with each other within 
experimental error. 


? A. Elliott, Proc. Roy. Soc. 127, 638 (1930). 
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Fic. 2. Absorption spectrum of Cl, (band 1-+12 of the 
12-—+*0 system). a is the spectrum of Cl; containing 45 per- 
cent Cl*’, b is the spectrum of Cl: containing 24.2 percent 
Cl*’, The rotational quantum numbers of the P branch are 
given for each molecule. 


The alternating intensity ratios for the two 
homonuclear molecules were obtained from the 
spectrum of the “‘enriched”’ chlorine. These ratios 
were found to be 1.26 for Cl®Cl*® and 1.28 for 
Cl*"Cl*’, This ratio may be calculated in two 
ways. The first way consists of comparing the 
average intensity of all the strong lines to the 
average intensity of the weak lines. The second 
method consists of finding the intensity ratio of 
each line to its immediate neighbors and then 
averaging over all lines. These two methods give 
the same result in this particular case. 

As was pointed out by Elliott the experimental 
values of the ratios should be regarded as lower 
limits since any background due to overlapping 
lines, etc., will make the measured ratio smaller 
than the true value. A background of less than 
20 percent will raise Elliott’s value 1.34 to 1.40 
which is the theoretical value for a nuclear 
spin of 5/2. 

One would expect this background effect to be 
greater in the spectrum of the “‘enriched”’ chlorine 
since the intensity of all the lines is relatively 
less. This seems to be in accord with the measured 
values of 1.26 and 1.28 as the alternating in- 
tensity ratios in the spectrum of “enriched”’ 
chlorine. However since the intensity of the 
spectrum of Cl*’Cl*’ is comparable to the in- 
tensity of the spectrum of Cl*CI* the relative 
effect of the background may be neglected. It 
would appear then that the nuclear spin of Cl*’ is 
the same as the nuclear spin of Cl*. The most 
probable value of the nuclear spin of Cl*® is 5/2 
as determined by Elliott? and by the author from 
the spectrum of ‘“‘normal”’ chlorine. 

The author is indebted to Professor W. W. 
Watson, who suggested this problem, for his 
guidance and helpful suggestions. 
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Investigations have been made on the natural conduc- 
tivity, and on the conductivity induced by gamma-rays in 
isooctane and liquid oxygen. The characteristics of the two 
types of conductivity are so very different that it may be 
concluded that the natural conductivity is not due to stray 
radiation or to cosmic rays. The curves for the natural 
conductivity are such that the current increases much 
faster than the field. Evidence is presented supporting 
Plumley’s theory that the observed conductivity is due to 


(Received June 21, 1943) 








dissociation of the liquids by the field, and against the 
view that it has a thermionic origin at the cathode, or the 
view that it is due to radiation. The prediction that the 
logarithm of the current should be a linear function of 
the square root of the field is verified. The much larger 
currents induced by gamma-rays rise at first rapidly and 
then more slowly in a manner indicating a slow approach 
toward saturation. The results fit Jaffé’s formula satis- 
factorily. 





INTRODUCTION 


ELATIVELY few investigations on the elec- 
trical conductivity of highly insulating 
liquids have been made. In comparison with the 
immense amount of work which has been done 
on the electrical conductivity occurring in, and 
imparted to, gases, it is evident that much re- 
mains to be done in accumulating information 
about the corresponding phenomena in highly 
insulating liquids. The first publication in this 
field was that of P. Curie,' who found that the 
electrical conductivity of petroleum ether, carbon 
tetrachloride, carbon disulphide, and benzene 
was increased by exposure to gamma- or x-rays, 
an investigation suggested, of course, by the 
effect of these rays on air. Later, other investi- 
gators studied the current-voltage characteristics 
of hexane and similar liquids when exposed to 
gamma-rays,’ the mobilities,’ the coefficient of re- 
combination,‘ and the coefficient of diffusion of 
the ions. Jaffé? made an important contribution 
to the understanding of the current-voltage 
characteristics of hexane when irradiated by 
gamma-rays. In his 1908 paper he considered the 


*From a dissertation submitted in partial fulfillment of 
the requirements for the Ph.D. degree at Washington 
University. 

tNow at the Radiation Laboratory, Massachusetts 
Institute of Technology, Cambridge, Massachusetts. 

1P. Curie, Comptes rendus 134, 420 (1902). 

2G. Jaffé, J. de Phys. 5, 263 (1906); Ann. d. Physik 25, 
257 (1908); Le Radium 10, 126 (1913); Ann. d. Physik 42, 
303 (1913). 

31. Adamczewski, Acad. Polonaise Sci. et Lettres, Bull. 
5~7A, 217 (1934); Acta Phys. Polonica 3, 235 (1934); Ann. 
de physique 8, 309 (1937). 

*H. J. Van der Bijl, Ann. d. Physik 39, 1, 170 (1912); 
C. Bialobrzeski, J. de phys. et rad. 7, 329 (1936). 








current to be accounted for as the sum of two 
separate currents; one of these rising to a satura- 
tion value like the current in a gas, while the 
second is an ohmic current. In 1913 he discarded 
this view in favor of a theory according to which 
the ions are produced in very densely packed 
columns isolated from each other. These columns 
are the tracks of the photoelectrons ejected by 
the gamma-rays. An important feature of this 
theory was that it gave a quantitative relation- 
ship between the current and the field, a relation- 
ship which expressed the great difficulty of 
securing saturation, as a result of the ease with 
which recombination occurs in the tightly packed 
columns of ions. This theory has been extended 
successfully to cases of ionization in compressed 
gases. The most recent application of Jaffé’s 
theory is that of Mohler and Taylor,® and later 
Taylor’ alone, who used it successfully in their 
work on ionization in carbon disulphide by 
x-rays. While there is no reason to suppose that 
Jaffé’s theory is inadequate, it is certainly de- 
sirable to extend the scope of investigations on 
the conductivity imparted to liquids by ionizing 
radiations. 

When all ionizing agents such as x-rays and 
radium are removed, it is found that there still 
exists a residual or natural conductivity in these 
liquids. However, there is no general agreement 
as to its characteristics. In the first place, the 


5 J. W. Broxon and G. T. Meredith, Phys. Rev. 54, 9 
(1938); J. Clay, Rev. Mod. Phys. 11, 123 (1939). 

*F. L. Mohler and L. S. Taylor, J. Research Nat. Bur. 
Stand. 13, 659 (1934). 
(1936) S. Taylor, J. Research Nat. Bur. Stand. 17, 557 
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residual conductivity of a supposedly pure insu- 
lating liquid can be reduced by a factor of a 
thousand or more by treating it in such a way 
as to remove traces of water. Such treatment 
may be prolonged refluxing with an agent like 
sodium or phosphorus pentoxide, or exposure to 
a high electric field. When a reasonably low 
limiting value of the conductivity has been 
reached, the next question is to find how it 
originates. Some authors* believe that the re- 
sidual conductivity should be attributed to 
cosmic and similar radiations. Other physicists® 
interpret their investigations to mean that the 
conductivity is due to a thermionic emission from 
the cathode combined with a Schottky effect. 
This has been criticized by LePage and Du- 
Bridge.'° The most recent work is that of 
Plumley" who interprets his experiments to sup- 
port the potential dissociation theory, which was 
originally proposed by Onsager” for very weak 
electrolytes. According to Plumley’s view there 
are even in the purest hydrocarbons such as 
hexane a very small number of ions resulting 
from the spontaneous dissociation of molecules. 
A favorable orientation of the molecules with 
respect to the field increases the number of 
dissociations, the number increasing rapidly with 
the strength of the field. The investigations, 
while often supporting the point of view of the 
various authors, are seldom complete enough to 
decide definitely against alternative theories. 
For example, in the experiments for which the 
thermionic point of view was suggested, there 
was no systematic investigation of how the cur- 
rent depended on the volume of the liquid be- 
tween the electrodes. For a given field, the 
thermionic theory would imply no dependence 
on the volume of the liquid, whereas all other 
theories would require for a given field a larger 
current, the greater the amount of liquid in- 
volved. However, it must be recorded that all 
who have worked on the residual conductivity 
of highly insulating liquids have found great 


8G. Jaffé, Ann. d. Physik 28, 326 (1909); M. C. Bia- 
lobrzeski, J. de phys. et rad. 7, 329 (1936); A. Rogozinski, 
tothe 128S (1937); A. Rogozinski, Phys. Rev. 60, 148 
* E. B. Baker and H. A. Boltz, — Rev. 51, 275 (1937). 
'W. R. LePage and L. A. DuBridge, Phys. Rev. 58, 
61 (1940). 
"H. J. Plumley, Phys. Rev. 59, 200 (1941). 
2L. J. Onsager, J. Chem. Phys. 2, 509 (1934). 
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Fic. 1. Experimental cell I: (1) amber plug, (2) brass 
rod, (3) spring, (4) brass tube, (5) graduated brass disk, 
(6) 20-ampere Buss fuse wire, (7) mixture of paraffin and 
stopcock greases, (8) packing, (9) metal flange, (10) silvered 
glass jar, (11) sylphon, (12) brass electrode, (13) outlet 
for the liquid, (14) inlet for the liquid. 


difficulty in securing reproducible results. Fre- 
quently, for no apparent reason at all, the value 
of the observed current will change by a factor 
of two or so. Getting the liquid into a state 
where it will retain reasonably constant charac- 
teristics is often a tedious trial and error process. 

This investigation was undertaken to study 
both the natural and induced conductivities 
under as varied conditions as possible in the 
hope that some of the ambiguities now existing 
as to interpretation could be removed. The 
liquids chosen for investigation were isooctane 
and liquid oxygen. The reason for choosing iso- 
octane was that we were fortunate to secure a 
supply of very pure liquid through the courtesy 
of Dr. G. W. Waters of the Shell Oil Company 
at Wood River, Illinois. The conductivity was 
investigated not only as a function of the field 
but also as a function of the spacing between the 
electrodes, or what amounts to much the same 
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Fic. 2. Experimental cell II: (1) fine nickel wire, (2) 
brass cylinder, (3) soft ion bar, (4) brass rod, (5) brass 
electrode, (6) fine nickel screen, (7) copper electrode, 
(8) — ring, (9) outlet for the liquid, (10) inlet for the 
iquid. 


thing, the volume of the liquid. Measurements 
were also carried out at different temperatures. 


APPARATUS 


The apparatus and methods employed may be 
conveniently discussed under several headings. 


1. The Conductivity Cells 


In all, four cells were used in this investigation. 
Two of them had parallel plate electrodes, the 
spacing between which could be adjusted so that 
the effect of distance between the electrodes on 
the observed current could be investigated. Two 
others had concentric cylindrical electrodes. One 
of these cells was easy to shield with lead for 
investigating a possible cosmic-ray origin, and 
the other was specially designed for investigating 
the effect of temperature. 
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Cell I: A diagram of this cell is shown in Fig. 1; 
the dimensions of interest are the radius of the 
electrodes, 2.4 cm, the pitch of the screw, 1/24th 
of an inch. By means of the sylphon it is possible 
to vary the distances between the electrodes 
from zero to 6 mm. The insulation of the upper 
electrode, which is the one connected to the 
electrometer, consists of two pieces of amber. 
The upper piece is made vacuum tight by 
Sealstix. The glass part of the apparatus is made 
from standard Pyrex glass pipe. The junction 
between its flange and the brass disk is made 
vacuum tight by a gasket of fuse wire fitted into 
the groove in the flange in place of the more 
usual rubber gasket. Care was taken to ensure 
that the faces of the two electrodes were accu- 
rately parallel. The inside of the glass was coated 
with evaporated silver, which was electrically 
grounded to act as a shield for the upper elec- 
trode. The liquid did not rise above the level 
of the side tube (13). 

Cell II: This all glass cell is shown in Fig. 2. 
The upper electrode, in this cell connected to the 
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Fic. 3. Experimental cell III: (1) nickel screen, (2) monel 
cylinder, (3) dural rod, (4) blobs of solder, (5) brass 
cylinder, (6) guard ring, (7) Dennison’s ‘“‘number 4 express” 
— wax, (8) inlet for the liquid, (9) outlet for the 
iquid. 
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high voltage, is a brass disk 2.88 cm in diameter 
which can be moved up or down by manipulating 
a magnet in the vicinity of the soft iron bar (3). 
The lower electrode (connected to the elec- 
trometer) is a copper disk supported by a strong 
tungsten wire which is sealed into a glass spider 
at the end of the inner glass tube (6); on the 
surface of the glass tube (6) is a grounded nickel 
screen. The metal cylinder (2) prevents a corona 
discharge from the fine wire (1) from giving 
trouble. Great care was taken to ensure paral- 
lelism between the electrodes. 

Cell III: This cell, designed to study the effect 
of external radiation, is shown in Fig. 3. It has 
two concentric cylindrical electrodes, the outer 
one of monel 9.6 cm long, 1.11-cm internal 
diameter, and the inner one of dural 9.2 cm long 
and 0.635-cm outside diameter. Blobs of solder 
serve to center the monel tube. Care was taken 
to secure exact centering. Cylinder (5) serves as 
a guard ring. The end (7) is covered with a good 
grade of sealing wax to ensure high insulation. 
The outside of the tube, above the sealing wax, 
is painted with Aquadag which is grounded. 
When the liquid is circulating in the system, it 
flows in through the tube (8) and leaves the 
tube through (9). 

Cell IV: The cell (Fig. 4) is designed to be in- 
serted in a thermos jar to allow investigations to 
be made at different temperatures. The dural 
rod (6) is 8.62 cm long and 1.68 cm in diameter. 
The monel cylinder (7) has an inside diameter of 
2.42 cm. There are guard rings both inside and 
outside the cylinder (4) which is evacuated. At 
(1) there is amber insulation. Continuous circu- 
lation of the liquid is not possible with this cell. 
The liquid is distilled into the cell through (8) ; 
it can be tipped back into the distilling flask 
and the process repeated. 


2. Treatment of the Isooctane 


In order to secure low and reasonably steady 
residual currents, it is necessary to remove the 
constituents (whatever they may be) which are 
associated with a high conductivity (say 10-* 
mho/cm) which is to be found even in the purest 
liquid obtainable commercially. The method 
finally adopted was to reflux the liquid over 
metallic sodium for a long time in a closed system 
including the conductivity cell. A simple design 
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of the distilling and refluxing system is shown in 
Fig. 5. By increasing the heat input to the dis- 
tilling flask, it was possible to arrange for the 
newly distilled liquid to pass through the con- 
ductivity cell back into the flask in a continuous 
cycle. In this way any suspended matter that 
might have been present in the cell initially was 
eventually removed. The conductivity sometimes 
obtained in this work is the lowest on record, 
3X10-*° mho/cm. (Since the relationship be- 
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Fic. 4. Experimental cell IV: (1) amber plug, (2) Aqua- 
dag film, (3) to diffusion pump, (4) Hanovia “Liquid Bright 
Platinum” paint serving as a guard ring, (5) brass cylinder, 
(6) dural rod, (7) aeedl cylinder, (8) glass tube connected 
to a two-way stopcock. 


tween current and voltage is not linear, no 
precise meaning can be attributed to this unit; 
it is, however, a short and convenient unit for 
comparing results.) 


3. Use of Liquid Oxygen 


Investigations of liquid oxygen were carried 
out only with cell IV. The apparatus was very 
carefully washed and then dried thoroughly. 














The tube (8) was connected through a suitable 
two-way stopcock which allowed one first to 
produce a vacuum in the apparatus and then to 
admit oxygen through a filter from a tank at a 
pressure of about 12 Ib./in.* above atmospheric 
pressure. On immersing cell IV in a thermos 
bottle of liquid air it was possible to fill it with 
liquid oxygen. It is believed that, though traces 
of water may be responsible for the conductivity 
in isooctane above the limiting value finally 
reached, they cannot be a cause of trouble in 
liquid oxygen for traces of water would be frozen 
to the walls and unavailable to form ions. 


4. Measurement of Current 


The currents measured ranged all the way from 
10-" down to 10-!* ampere. They were measured 
by means of the Western Electric D-96475 elec- 
trometer tube mounted in a vacuum. The grid 
of the tube was grounded through one of a 
set of S. S. White resistors whose values ranged 
from 4.06 X10" to 2.4 10° ohms. The tube was 
used in the well-known balanced circuit set-up 
described by Penick." A sensitivity of 200,000- 
mm divisions per volt was obtained, so that with 
a resistor of 4.06 X 10" ohms, a deflection of 1 mm 
indicated a current of 1.2X10-'7 amp. This 
sensitivity was seldom used except when the 
current to be measured was extremely low. For 
high currents care was taken to use a resistor 








Fic. 5. Disti'lating system: (1) to Wood’s metal seal off 
device, (2) indentations, (3) entrance for the liquid and 
sodium, (4) glass tube for equalizing the pressures between 
the cell and the reservoir, (5) inlet for the liquid from the 
cell, (6) outlet for the liquid to the cell, (7) constriction, 
(8) capillary. 


1831D). B. Penick, Rev. Sci. Inst. 6, 115 (1935). 
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between the grid and the cathode which limited 
the voltage change on the grid to less than 0.05 
volt. This was necessary to secure linearity in 
the system. 


5. High Voltage Supply 


The investigation included measurements on 
the current through the liquids with voltages as 
high as 6000 volts. Slight fluctuations in the 
voltage applied to the cells can produce effects 
which will completely mask a small current. 
Constant voltage can be secured by means of a 
sufficient number of B batteries. In this investi- 
gation voltages as high as 1485 volts were pro- 
vided by means of B batteries. It was im- 
practicable to use B batteries for higher voltages, 
An electronic high voltage supply was designed 
and constructed which had the following charac- 
teristics: (1) Constant voltages up to 5500 volts 
could be produced at will. The constancy was 
comparable to that of a bank of B batteries, 
Tests showed that the output voltage changed 
by less than 0.1 volt in a time interval of ten 
minutes. The fluctuations were often consider- 
ably less than this. (2) The device gave high 
voltages of either sign with no difference in the 
constancy of the output. (3) No detectable ripple 
could be observed. An account of the high voltage 
supply will be published elsewhere. 


EXPERIMENTAL RESULTS 
Induced Conductivity (Isooctane) 
1. Dependence on Field and Volume 


The conductivity induced by gamma-rays in 
highly purified isooctane was investigated in all 
the cells. Only with cells I and II was it possible 
to study the way in which the current depends 
on the spacing for any given electric field. (It was 
possible of course to investigate the dependence 
on field with all four cells.) Typical results for 
the current obtained with isooctane in cell Il 
irradiated by gamma-rays are shown in Fig. 6. 
It is clear that the current is almost, but not 
quite, independent of the direction of the electric 
field. This asymmetry is a spurious effect and is 
to be attributed to the small contribution to the 
observed current from the ‘‘dead”’ region around 
the back of the electrode which is connected to 
the electrometer. The direction of the current in 
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this region is not affected by the sign of the 
voltage on the other electrode; hence, the small 
asymmetry. This interpretation was verified by 
using one cell with a lead shield so arranged that 
the direct ionization of the ‘‘dead” region by 
the radium is greatly reduced; it was found that 
this shielding reduced the small asymmetry very 
greatly. The current increases noticeably with 
the field showing that saturation calls for much 
higher fields. Results of this kind were obtained 
by Jaffé? with hexane many years ago. 

He proposed an interpretation according to 
which the ionization due to the photoelectrons 
ejected here and there in the liquid by gamma- 
rays takes place in very short columns densely 
packed with ions. The difficulty of obtaining 
saturation is due to the strong initial recombina- 
tion within the columns. Jaffé developed a 
formula expressing the value of the current at 
any field in terms of the mobility, diffusion, and 
recombination coefficients of the ions. The for- 
mula is 


1 1 aN 1 b*w?X? 

iE) « 

Ie te 7.85bw X 2D? 
where J, is the current per unit area at a field X, 
I,, the saturation current at infinite field, N the 
number of ion pairs produced per cm length of a 
column, 6 the initial diameter of the column 
within which the ions are found, and a, w, 
and D are the recombination, mobility, and 
diffusion coefficients, respectively, of the ions. 
F(b’w?X?/2D*) is a complicated function of the 
argument which approaches unity when the argu- 
ment is sufficiently large. In practice this is 
secured by using a high enough field so that 
Eq. (1) becomes 


1 1 1 
~=—(1+5-). (2) 
I De P 4 


Here S is aN /7.85bw. The validity of this theory 
may be tested by plotting 1/J4 against 1/X. 
If the plot results in a curve which is practically 
straight for all values of X above a certain value 
then one can infer that the theory applies. 
Moreover, the intercept of this line (when ex- 
trapolated) on the ordinate through the origin 
gives the value of 1/J,. It is thus possible, 
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Fic. 6. Variation of induced current with field strength 
in isooctane irradiated by constant radiation from 2000 ug 
radium. Curves A, B, C, D, E, F, and G represent results 
when d=0.848, 0.636, 0.371, 0.265, 0.159, 0.106, and 
0.053 cm. Open circle and cross represent the respective 
results for positive and negative voltages applied to the 
high voltage electrode of cell II. 


assuming the validity of Jaffé’s theory, to calcu- 
late the value of the saturation current. The 
slope of the line in such a graph has the value 
S/T. 

The. directly observed experimental results 
collected in Fig. 6 are replotted in Fig. 7. It is 
clear that these “‘Jaffé plots’’ have straight por- 
tions which can be extrapolated back to the 
ordinate at the origin which corresponds to an 
infinite field. The characteristics of this family 
of curves are presented in Table I. Here d is the 
spacing between the electrodes. The fourth 
column is quite significant for it shows that J,, is 
almost exactly proportional to d, as it should be 
if the total number of ion pairs is proportional to 
the available volume. On the other hand, the 
fifth column shows that the currents observed 
for an experimentally obtainable field X = 20 
kv/cm are not quite proportional to the spacing 
between the electrode. It is to be expected that 
with the larger spacings, the currents would 
deviate further from the saturation value, as is 
found to be the case. 


2. Dependence on the Intensity of Radiation 


In order to determine the intensity variation 
of the current induced in the liquid at room tem- 
perature, current-intensity runs were taken for 
different fields in a cell of concentric cylinder 
type, whose construction is similar to that of 
cell III. The intensity of radiation was varied by 
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Fic. 7. Inverse current and inverse stress (field strength) 
relationship for isooctane irradiated by constant radiation 
of 2000 wg radium. Curves A, B,.C, D, and E represent the 
results when d=0.053, 0.106, 0.159, 0.265, 0.371 cm, re- 
spectively. Open circle and cross represent the results when 


positive and negative voltages are applied to the high 
voltage electrode of cell II. 


using different calibrated gamma-ray standards 
supplied by the Bureau of Standards. Figure 8 
shows what one obtains when the intensity of 
radiation is selected and the current measured 
at various fields. It is clear that the present re- 
sults agree with those obtained by Adamczewski," 
who measured the ionization current through 
hexane with various intensities of x-rays. 

With a fixed value of field strength the induced 
current was measured as a function of the in- 
tensity of radiation. A typical curve is shown in 
Fig. 9. It is a rather remarkable phenomenon 
that the relationship between the induced current 
and intensity of radiation is quite linear, a fact 
which, in principle, enables us to use the ioniza- 
tion in such a liquid to indicate the intensity 
of radiation. If the radiation is sufficiently strong 
and if the liquid inside the ionization chamber is 
fairly pure, so that the induced current is much 
greater than the residual current, the author 
believes that this kind of measurement may be 
developed into a dependable and accurate method 
of measuring radiations. 


“|, Adamczewski, Acta Phys. Polonica 3, 246 (1934). 
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3. Dependence on Temperatures 


The effect of temperature on the currents 
through isooctane as a function of the field was 
studied by means of cell IV. Figure 10 shows the 
results obtained when the field was set at some 
selected value and the temperature varied. 
Figure 11 shows what one obtains when the tem- 
perature is selected and the current measured at 
various fields. It is clear that for any given field 
the currents are larger the higher the tempera- 
tures. This can be explained qualitatively by 
the fact that viscosity decreases with increasing 
temperature and when viscosity decreases, the 
ions in the column are more readily separated 
from each other so that closer approach to the 
saturation current is obtained. 

The results in Fig. 11, together with some 
others not shown, are replotted as Jaffé plots in 
Fig. 12. It is interesting to note that all the lines, 
when extrapolated, hit the ordinate at the origin 
at substantially the same place, indicating that 
for each curve we have the same J,. This is as 
it should be, for with the radium placed in a 
fixed position with respect to the cell, one would 
expect the same total number of ion pairs to be 
produced whatever the temperature may be. 
(It.is true that a change in density with tempera- 
ture could affect the absorption of gamma-rays 
slightly and this would change the total number 
of ion pairs produced per unit volume. However, 
this is probably too small an effect to make a 
noticeable change in J.) 

From Fig. 12 we find an empirical relation be- 
tween the value of S and T to be S=const./T' 8, 
Since S=aN/7.85bw it is likely that N and b are 
but little affected by the temperature so that 
the principal factors varying with the tempera- 
ture are a and w, the recombination and mobility 
coefficients, respectively. Thus measurements of 
the slopes of the straight lines in the Jaffé plots 


TABLE I. Characteristics of curves of Fig. 7; 
d is the space between electrodes. 














10-11 Tm X107! In 
d Teo IaXi07# d Teo S X10" 
0.053 6.6 1.51 2.85 0.550 10.9 
0.106 3.3 3.02 2.845 0.465 7.3 
0.159 2.2 4.53 2.845 0.412 6.2 
0.265 1.33 7.51 2.83 0.397 4.03 
0.371 0.96 10.04 2.81 0.405 2.83 
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can be used to indicate the way in which a/w 
changes with the temperature. 

In cells containing electrodes of different ma- 
terials, the magnitude of the current was quite 
independent of which electrode was the cathode, a 
result entirely in accord with the view that the 
conductivity due to gamma-rays is a volume 
effect in the liquid and not a surface effect at 
the electrodes. 

No time lag in the value of the observed cur- 
rent was noticed when the field was changed. 
This result is not in agreement with that found 
by Taylor’ for ionization of carbon disulphide by 
x-rays. Our statement is subject to the proviso 
that conclusions cannot be drawn as to time 
intervals less than about three minutes as the 
resistance condenser factor of our set-up was of 
this order. 

It may be concluded that the phenomena ob- 
served when isooctane is irradiated by gamma- 
rays may be attributed to the formation of ions in 
the liquid and that the results are described ade- 
quately by Jaffé’s theory. 


Induced Conductivity (Liquid Oxygen) 


Because of the experimental restrictions im- 
posed by the low temperature of liquid oxygen, it 
was impracticable to use either cell I or II to 
determine the effect of changing the spacing. 


“or 


ps — 


I, IN 10°"? AMP7CM2 








B 
Cc 
2) 
i r i i 1 1 a 
0 I 2 3 4 5 6 7 
X IN KV/CM 


_ Fic. 8. Variation of induced current with field strength 
in isooctane at room temperature (7 =297°K). Curves A, 
B, C, and D represent the respective results when the 
amounts of radium are 100, 50, 20, and 10 ug. 
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Fic. 9, Variation of induced current with intensity 
of radiation at X = 5830 v/cm. 


Cell IV was used for studies on liquid oxygen. 
The direct experimental measurements are 
plotted in Fig. 13. (To facilitate comparison, the 
light dotted line shows the curve for isooctane at 
273°K in the same cell. The curves are adjusted 
arbitrarily to coincide at 10 kv/cm.) The curve 
for liquid oxygen indicates that saturation is 
much harder to attain in it than in isooctane. 
The corresponding Jaffé plot is shown in Fig. 14. 
The intercept on the ordinate at the origin is so 
small that it would be almost useless to use the 
curve to determine the value of J... However, the 
value of J, is of the order 30 times more than J, 
at 10 kv/cm, whereas, in the case of isooctane, 
the value of J,, is only of the order 2 times more 
than J, at 10 kv/cm (for T=273°K). If the 
273°K curve for isooctane is transferred to 
Fig. 14 and adjusted so that J, is the same for 
both, we find that the Jaffé plot for isooctane is 
the lower curve. The S for liquid oxygen is thus 
much larger than the S for isooctane. 

The high value of S can be attributed princi- 
pally to the factors aN/w. N may well be con- 
siderably larger for liquid oxygen than for 
isooctane because the oxygen molecules offer a 
larger cross section to the gamma-ray photo- 
electrons than do the isooctane molecules. Little 
can be said as to the quotient a/w since no 
measurements have been made to determine these 
constants in liquid oxygen. 


° Natural Conductivity (Isooctane) 
1. Dependence on Field and Volume 


The investigation of the natural conductivity 
in isooctane and similar liquids, that is, the 
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Fic. 10. Variation of induced current with absolute tem- 
perature in isooctane irradiated by constant radiation from 
4000 zg radium. Curves A, B, and C represent the results 
when X=14,850, 8088, and 1348 v/cm, respectively. 
Open circle represents the results when positive voltages 
are applied to the high voltage electrode of cell IV, and 
cross represents the results when negative voltages are 
applied to the high voltage electrode. 


residual conductivity to be found in the absence 
of ionizing agents such as radium, is full of ex- 
perimental difficulties. The currents to be meas- 
ured are so small that it is difficult to secure 
reproducible results. This is also accentuated by 
the need for measurements at high voltages as 
well as at low voltages. Slight variations in the 
absolute value of the high voltage, which may be 
quite insignificant in a relative sense, often pro- 
duce electrostatic induction effects of a magni- 
tude comparable with the current being measured. 
There is no agreement among previous investi- 
gators as to the nature and origin of the natural 
conductivity. It has been attributed to (1) radia- 
tion, probably cosmic rays, (2) thermionic emis- 
sion from the cathode, (3) cold emission from the 
cathode,'® and (4) dissociation induced by the 
field. In order to effect a discrimination between 
these views it is desirable to study the conduc- 
tivity as a function of several variables in one and 
the same investigation. Thus for a given field, the 
current should be independent of the spacing 
between the electrodes, if the explanation be 
either (2) or (3). The supporters of these theories 
have not made this particular test. 
The dependence of the current through iso- 


1% R. W. Dornte, Ind. Eng. Chem. 32, 1529 (1940). 
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Fic. 11. Variation of induced current with field strength 
in isooctane irradiated by constant radiation from 4000 yg 
radium. Curves A, B, and C represent the results when 
T=315.3, 260.6, and 210°K, respectively. Open circle 
represents results when positive voltages are applied to 
the high voltage electrode of cell IV, and cross represents 
results when negative voltages are applied to the high 
voltage electrode. 
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Fic. 12. Variation of inverse current and inverse stress 
(field strength) curves with the absolute temperature in 
isooctane irradiated by constant radiation from 4000 yg 
radium. Curves A, B, C, D, E, and F represent the results 
when 7 =313.3, 273, 260.6, 230, 210, and 191°K, respec- 
tively. Open circle and cross represent the results when 
positive and negative voltages are applied to the high 
voltage electrode of cell IV. 


octane on spacing was investigated in cells I and 
II. Although the same amount of effort was spent 
on the experiments involving these two cells, the 
ultimate conductivity (at the same field) ob- 
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tained after prolonged refluxing and circulation 
was about 20 times less in II than in I. This is 
typical of much of the work in the field. Many 
cells with refluxing and circulating systems have 
been tried out. Often two systems, substantially 
identical, yield widely different ultimate con- 
ductivities. Figures 15 and 16 give directly the 
values of the currents per unit area as a function 
of both the spacing and the field. It is to be noted 
that the current is independent of the sign of the 
potential applied to the high voltage electrode. 
This was also found later with cells III and IV 
where the electrodes are made of different metals. 
Had the current been due to a thermionic emis- 
sion at the cathode, it would have depended very 
greatly on the nature of the cathode. We may 
conclude, therefore, that these experiments give 
no support to explanations (2) nor indeed to (3). 
Another characteristic to be noted is that the 
curves for the natural conductivity are entirely 
different in shape from the curves for the con- 
ductivity induced by gamma-rays. This result 
makes it extremely unlikely that one can attri- 
bute the residual conductivity to feeble radia- 
tions of the gamma-ray type. Figure 17 shows 
that the currents observed in cell I are strictly 
proportional to the spacing between the elec- 
trodes, or what amounts to the same thing, to the 
volume of the liquid between the electrodes. 
Figure 18 shows how the currents in cell II de- 
pend on the spacing. While it is clear that these 
currents also definitely increase with the spacing, 
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Fic. 13. Variation of induced current with field strength. 
Curve A (dotted curve) represents the results when iso- 
octane was irradiated by 4000 yg radium at 7 =273°K. 
Curve B (solid curve) represents the results when liquid 
oxygen was irradiated by 4500 yg radium at 7=90°K. 
Open circle and cross represent the respective results for 
positive and negative voltages applied to the high voltage 
electrode of cell IV. 
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they are not strictly proportional to the spacing. 
It may be that this is to be attributed to the 
great difficulty in making measurements in this 
region. (Over one-half of points in Fig. 18 
correspond to J, less than 10-' amp./cm*.) On 
the other hand, since the currents in Fig. 18 are 
so much smaller than those in Fig. 17 it may well 
be that the way in which the lines pass above the 
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Fic. 14. Inverse current and inverse stress (field strength) 
curves. Curve A (dotted curve) represents the results 
when isooctane was irradiated by 4000 ug radium at 
T =273°K. Curve B (solid curve) represents the results 
when liquid oxygen was irradiated by 4500 yg radium at 
T=90°K. Open circle and cross represent the respective 
results for positive and negative voltages applied to the 
high voltage electrode of cell IV. 
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Fic. 15. Variations of residual current and field strength 
curves for isooctane. Curves A, B, and C represent the 
results when d=0.318, 0.424, and 0.530 cm. Open circle 
and cross represent the results when positive and negative 
— are applied to the high voltage electrode of 
cell I. 
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Fic. 16. Variation of residual current and field strength 
curves for isooctane. Curves A, B, C, and D represent the 
results when d=0.053, 0.106, 0.265, and 0.636 cm, re- 
spectively. Open circle and cross represent the results 
when positive and negative voltages are applied to the 
high voltage electrode of cell IT. 


origin indicates the superposition of two effects, 
a small one (possibly due to radiation) whose 
characteristic curve cuts the zero ordinates above 
the origin and a larger one which is proportional 
to the spacing. In Fig. 17 the larger effect is 
possibly so enhanced that the small effect is 
swamped. It is perhaps unprofitable to speculate 
further on the point. Clearly more experimental 
work should be done. 

Plumley" has applied the potential dissociation 
theory to explain residual conductivity. Ac- 
cording to this view a very small number of the 
molecules even in the purest hydrocarbons are 
dissociated into ions. Other molecules on the 
verge of dissociation will be dissociated by a field 
if their orientation be favorable. Plumley gave a 
formula for the number of dissociated molecules, 
which, therefore, also leads to an expression for 
the current in terms of the temperature T and 
the field X. The formula is 


I4=C exp [(2¢e/kT) (eX /300D)!]. (3) 


Here C is proportional to the number of dis- 
sociated molecules at zero field. The exponential 
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term describes the way in which a field X in- 
creases the degree of dissociation ; k is Boltzmann's 
constant, e the unit charge, and D the dielectric 
constant. This theory can be tested by plotting 
logy Ia against X?. If it is correct, a straight line 
should be obtained and its slope should be 
2e!/2.3kT300'!D!. The ‘‘2” enters the equation 
through a particular assumption as to the nature 
of the force between the ions. If a straight line 
with a slope differing from the calculated value 
be obtained from the experimental results, it may 
be inferred that while the theory, in the main, is 
satisfactory, the correct form of the force be- 
tween the ions has yet to be found. 

In Figs. 19 and 20 we have re-plotted the re- 
sults given in Figs. 15 and 16 for cells I and II, 
respectively. Insofar as the plots are straight 
lines one may infer that the potential dissociation 
theory is satisfactory. The slopes of the lines in 
Fig. 19 are substantially the same, a result which 
we should expect for there is nothing in the 
theoretical expression for slope which depends on 
the spacing between the electrodes. In these 
observations, however, the variation in spacing 
is not very great. When this variation is made 
large, the curves of Fig. 20 are obtained; here the 
slopes are very different for the different spacings, 
This cannot be accounted for by Plumley’s 
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Fic. 17. Variation of residual current with electrode 


spacing at constant field strength. Curves A, B, and C 
represent the results when X =10, 8, and 4 kv/cm, re- 


spectively. 
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theory. The theoretical value of the slope 
2e8/2.3kT300'D# is 0.0092 at 297°K. For cell I 
(Fig. 19) the experimental value is 0.0215, while 
for cell II (Fig. 20) it ranges from 0.0105 to 


0.0183. 


2. Dependence on Treatment of the Isooctane 


Observations were made on the residual cur- 
rent through cell II during a period of six weeks. 
This cell together with its distilling flask forms a 
closed system. The first observations were made 
after circulation for about twenty-four hours. 
Results about three times larger in magnitude 
than those represented in Figs. 16 and 18 were 
obtained. During the course of repeated distilla- 
tion and circulation the current values slowly 
diminished to the limiting values shown in Figs. 
16 and 18. The slope of the log I, against X! lines 
also diminished by roughly 20 percent. 

When the liquid had been “‘purified” by pro- 
longed distillation over sodium, exposure to a 
high electrostatic field for an hour did not pro- 
duce any further reduction. 


3. Time Lag 


Some investigators report hysteresis effects in 
the current when the field is charged. No time 
lags of this sort were found in this investigation 
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Fic. 18. Variation of residual current with electrode 
spacing at constant field strength. Curves A, B, and C 
represent the results when X=20, 14, and 8 kv/cm, 
respectively. 
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Fic. 19. Relation between logio J4 and X# for isooctane. 
Curves A, B, and C represent the results when d=0.318, 
0.424, and 0.530 cm. Open circle and cross represent the 
results when positive and negative voltages are applied to 
the high voltage electrode of cell I. 
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Fic. 20. Relation between logio J4 and X# for isooctane. 
Curves A, B, C, and D represent the results when d=0.053, 
0.106, 0.265, and 0.636 cm, respectively. Open circle and 
cross represent the results for positive and negative volt- 
ages applied to the high voltage electrode of cell II. 


larger than the purely instrumental ones associ- 
ated with the resistance capacitance character- 
istics of the circuit. 


4. Effect of External Radiation 


It was possible to surround cell III with 3 
inches of lead and an extra inch on the top. No 
measurable change in the current at any field 
could be detected. On replacing cell III by a 
Geiger-Miiller tube it was fourid that its response 
was cut down by 22 percent by the presence of 
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Fic. 21. logioJ« against 1/7 curves for isooctane. 
Curves A, B, and C represent the results when X = 14,850, 
4010, and 1215 v/cm. Open circle and cross represent the 
respective results for positive and negative voltages applied 
to the high voltage electrode of cell IV. 


the lead. It may, therefore, be inferred that little 
or none of the current through the isooctane in 
cell III could be attributed to cosmic radiations. 
Other investigators* have published results in 
support of the view that the residual conductivity 
.in their experiments is due to radiation which can 
be reduced by shielding. The data are summarized 
in Table II. 

It is difficult to understand why no evidence of 
an effect due to stray radiations was found in this 
investigation. The ‘“‘conductivity” of the liquid 
was as low as that used by the other investi- 
gators. There is a difference which may perhaps 
be important. The volume of the cell used in this 
work was far smaller than in the cells used by the 
other investigators, while the field was larger. It 
is conceivable that had the current at fields below 
1 kv/cm been large enough to measure accu- 
rately, a dependence on shielding would have 
been found. 

Assuming that the residual conductivity of air 
is due to the formation of 3 ion pairs per cc per 
sec., and making reasonable assumption as to the 
method of predicting the ionization in isooctane, 
as suggested by Mohler and Taylor’s work,® it 
may be concluded that the current to be expected 
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in isooctane is 6X10-'* amp./cm*. This is of the 
order of magnitude of the observed current below 
about 1 kv/cm, but unfortunately it is not easy 
to make decisive measurements on the effect of 
shielding when the current is so small. To secure 
further evidence in this direction it would be 
helpful to use a cell containing a much larger 
volume of liquid. 


5. Effect of Temperature 


Equation (3) indicates theoretically how J, 
should depend on the temperature. According to 
this equation the relation between log J, and 
1/T should be linear. The experimental results 
are shown in Fig. 21. The most noticeable thing 
about these curves is the change in slope for the 
lower part of the temperature range. This has the 
effect that the current at low temperatures is 
larger than theory predicts. 

It should be pointed out here that in the 
temperature range from 318°K to 263°K the 
slopes of logy Z4 against 1/7 are of the same 
order of magnitude as those for toluene measured 
by LePage and DuBridge.'® The only difference 
is that the present slopes of these plots decrease 
with a decrease of field, while those obtained by 
LePage and DuBridge increase with decreasing 
fields, as actually required by the potential 
dissociation theory. Again this difference can be 
reconciled by assuming that the observed current 
is composed of two effects, one due to cosmic 
radiation and the other due to the potential 
dissociation. The current measured by LePage 
and DuBridge is about 50,000 times larger than 
the current measured in the present investigation 
(at the same X and the same 7) ; thus, one would 


TABLE II. Evidence for ionization of liquids 
by natural radiation. 

















Jaffé Bialobrzeski Rogozinski Pao 
v 200 800 — 5.98 
X 1500 900 -- 19000 
ry 0.15 2 — 0.317 
ro 2.00 4 _ 0.555 
I, 216 600 135 190000 
t 3 5 7 8 
R 0.66 0.61 0.28 0.00 








In this table the following symbols are used: v =the effect volume of 
the liquid in cubic centimeters; X =the maximum field strength in 
volt/cm used in the investigation; r: =radius of the inner electrode in 
cm; r2=radius of the outer electrode in cm; J» =the current density 
in number of ions/cc «sec. at X; ¢ =thickness of the lead used in cm; 
R =effect of the external radiation. 
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expect that, in their measurements, the po- 
tential dissociation effect overrides the other to 
show a result in accordance with the theory, while 
in our measurement the two effects may be of the 
‘same order of magnitude. By calculating the 
number of ions per cubic centimeter per second at 
T =273°K from Fig. 21, we found these values to 
be of the same order of magnitude as the value 
calculated for the cosmic ray alone, which is 
about 4240 ions/cc: sec. This tends to support the 
argument given above. 

The relationship between J, and X at different 
temperatures is shown in Figs. 22 and 23. As T 
decreases from 309°K to 195°K, the slope de- 
creases almost 13 percent. Theoretically, if Eq. (3) 
applies, the slope should increase about 38 
percent. 

The assumption made in this method of testing 
Eq. (3) is that every ion pair produced finds its 
way to the electrodes, which is the same thing as 
saying that we have a saturation current. 
Plumley" has established a criterion for testing 
whether or not saturation conditions prevail. On 
applying this test it would appear that we actu- 
ally do have saturation. However, as the temper- 
ature is lowered the mobility of the ions de- 
creases while the coefficient of recombination 
increases. Both factors act so as to diminish the 
degree of saturation as the temperature is 
lowered. This explanation is still to be confirmed. 


30 ° 











x IN KvAem 





Fic. 22. Variation of residual current with field strength 
at constant temperature. Curves A, B, C, and D represent 
the results when 7'=309.5, 296.7, 273, and 195°K, re- 
spectively. Open circle and cross represent the results for 
positive and negative voltages applied to the high voltage 
electrode of cell IV. 
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Fic. 23. logio J4 against X! curves for isooctane. Curves 
A, B, C, and D represent the results when 7 = 309.5, 
296.7, 273, and 195°K, respectively. Open circle and cross 
represent the results when positive and negative voltages 
are applied to the high voltage electrode of cell IV. 








Natural Conductivity (Liquid Oxygen) 






The current passing though liquid oxygen as a 
function of the field in the absence of radium is 
shown in Fig. 24. Cell IV was used so that these 
results are directly comparable with those for 
isooctane in Fig. 22. The outstanding results here 
are that, when the same cell is used, the current in 
liquid oxygen is much less than that in isooctane 
for fields below about 5 kv/cm; for higher fields, 
say 15 kv/cm, the currents are of the same order 
of magnitude. As a matter of interest, we tabulate 
the conductivity of isooctane obtained from 
curve C, Fig. 22, and that of liquid oxygen in 
Table ITI. 




















TABLE III. Comparison of natural conductivities of iso- 
octane and liquid oxygen. 






























Isooctane (JT =273°K) Liquid oxygen (T =90°K) 
xX Ga xX Ga 

1348 7.3 1214 0.252 

4050 7.9 4010 0.307 

6760 11.0 6760 2.495 
9460 12.5 9460 4.52 
12140 17.9 12140 7.61 
14860 26.1 13520 9.84 
14860 14.70 








X is the field strength in volt/cm; oa is the average of positive and 
negative conductivities in 10- mho/cm. 
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Fic. 24. Variation of the residual current with field 
strength in liquid oxygen. Open circle and cross represent 
the respective results for positive and negative voltages 
applied to the high voltage electrode of cell IV. 


In order to compare the slopes of the straight 
lines representing log J4 against X! for liquid 
oxygen and isooctane the results shown in Fig. 24 
are replotted. In Fig. 25 curve A represents the 
results for isooctane at 7=273°K, curve B for 
liquid oxygen at T=90°K. It is noted that the 
liquid oxygen has a greater slope as predicted by 
the potential dissociation theory because of its 
lower temperature. The experimental slope for 
curve B is 3.56 10-?, while the theoretical value 
is 3.03 X 107. 

In working with isooctane at ordinary temper- 
atures there is always the slight possibility that 
the conductivity is due to the ions into which the 
last irremovable traces of water dissociate and 
not to charged molecules or fragments of mole- 
cules of isooctane. Such an explanation for the 
behavior of liquid oxygen cannot be sustained for 
it seems certain that any water that might 
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Fic. 25. Relation between logio Ja and X*. Curve A 
represents the results for isooctane at T=273°K and 
curve B represents the results for liquid oxygen at 
T=90°K. Open circle and cross represent the respective 
results for positive and negative voltages applied to the 
high voltage electrode of cell IV. 


happen to get into the apparatus would be 
firmly frozen to the solid surfaces. 
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in the Short Infra-Red Spectrum 


S. Breen,* J. N. HumpHRey, AND L. R. INGERSOLL 
Department of Physics, University of Wisconsin, Madison, Wisconsin 
(Received May 22, 1943) 


The magneto-optical activity of thin crystals of nickel sulphate a-hexahydrate has been 
measured by means of the spectroradiometric method in the spectral range 0.74 to 1.9% and 
for ordinary and low (—125°C) temperatures. The curves showing the dispersion of the 
rotation have marked anomalies in the region of \ 0.694 and particularly 1.164, matching the 
ones found in the optical rotation for these wave-lengths. The ratio of the rotations at low 
temperatures and room temperature is itself an anomalous curve. The effects can be qualita- 
tively explained on the basis of the inverse Zeeman effect. 





HE interesting anomalous optical properties 

of thin crystals of nickel sulphate a-hexa- 
hydrate have been studied at some length, par- 
ticularly as regards optical activity, circular 
dichroism, and absorption in the wave-length 
range 0.25y-2.34 at ordinary’ and low® tempera- 
tures. The magnetic rotation has also been 
measured in the visible* and ultraviolet‘ regions. 
The present work is a study of the magnetic 
rotatory power (Faraday rotation) in the short 


infra-red, particularly in the region of the ab- 
sorption band at A 1.164, where the optical 
activity shows such a striking anomaly. This has 
been done for both ordinary and low tem- 


peratures. 


EXPERIMENTAL METHODS 


The spectroradiometric method of measuring 
magnetic rotatory dispersion has been described’ 
several times and need be only briefly mentioned. 
Light from a special tungsten strip filament lamp 


*Now at Radiation Laboratory, Massachusetts Insti- 
tute of Technology, Cambridge, Massachusetts. 

1N. Underwood, F. G. Slack, and E. B. Nelson, Phys. 
Rev. 54, 355 (1938) (see this article for earlier references) ; 
F. G. Slack and Philip Rudnick, Phil. Mag. 28, 241 (1939) ; 
L. R. Ingersoll, P. Rudnick, F. G. Slack, and N. Under- 
wood, Phys. Rev. 57, 1145 (1940). 

*P. Rudnick, F. G. Slack, and J. J. O’Connor, Phys. 
Rev. 58, 1003 (1940); W. C. Knopf, Jr., and W. C. Gil- 
more, Jr., J. Opt. Soc. Am. 32, 619 (1942); P. Rudnick 
and L. R. Ingersoll, J. Opt. Soc. Am. 32, 622 (1942). 

*F. G. Slack, R. T. Lageman, and N. Underwood, Phys. 
Rev. 54, 358 (1938). 

*J. J. O'Connor, C. Beck, and N. Underwood, Phys. 


Rev. 60, 443 (1941). 

5L. R. Ingersoll, J. On. Soc. Am. 5, 157 (1921); Phil. 
Mag. [6] 18, 74 (1909). For description of the thermopile- 
galvanometer system in use at present see L. R. In ll 
and W. R. Winch, Phys. Rev. 44, 401 (1933); L. R. 


Ingersoll, J. Opt. Soc. Am. 27, 412 (1937). 
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was concentrated by a spherical aluminized glass 
mirror into a narrow conical beam, polarized by 
passage through a double-image prism (only one 
image being used) with principal plane at 45° 
azimuth. This light traversed the conically bored 
polepieces of a'large electromagnet, converging 
to a focus on the crystal specimen between the 
poles. It was then in turn focused by a second 
concave mirror, after passage through a large 
Wollaston calcite double-image prism (whose 
principal planes were vertical and horizontal) on 
the double slit of the spectrometer. The two 
slits were vertically above one another, and one 
of the two images from the Wollaston fell on 
each slit. After passage through the spectrometer, 
which used a 9-cm high heavy glass 45° prism 
with Wadsworth mirror mounting and two con- 
cave mirrors of 40-cm focus, the two beams fell 
on the two arms of a vacuum differential thermo- 
pile which was connected with a Kipp and Sons 
ZC galvanometer. When the system was balanced 
(i.e., the two beams of equal intensity), with no 
magnetic field, the galvanometer spot stood at 
zero, but exciting the magnet, with the crystal 
between the polepieces, caused a small Faraday 
rotation of the plane of polarization of the beam 
which fell on the Wollaston prism and the re- 
sulting unbalance gave a galvanometer deflec- 
tion. To interpret this deflection in degrees or 
minutes of rotation, it was compared with that 
produced by the mechanical rotation of the 
polarizer through 1°. Measurements could be 
made in general for wave-lengths between 0.5u 
and 2.3 although the results in the present case 
are confined to the region 0.74—-1.9u. 
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Fic. 1. Crystal holder and method of cooling. C, crystal; 
P, magnet polepieces ; J, Inconel caps. 


The method of keeping the crystal specimen 
at a suitable low temperature requires a little 
description. Preliminary experiments involving 
various ways of cooling the crystal by more or 
less intimate contact with liquid air in general 
proved unsuccessful, the chief difficulty en- 
countered being the film of condensed moisture 
on the specimen or windows or both. The problem 
was finally solved by mounting the crystal 
(Fig. 1) on a flat brass bar, the end of which 
extended into a flask of liquid air and past which 
was flowing a stream of cold evaporated air. To 
avoid the disturbing effect due to the proximity 
of the relatively warm polepieces, these were 
sheathed with specially designed Inconel caps. 
The resulting temperature, which was measured 
by two small calibrated thermocouples in contact 
with the bar just above and below the crystal, 
could be maintained very steady at —125°C, 
and this was quite satisfactory for our purpose. 
The room temperature measurements were made 
at 25°C. Suitable water-cooling arrangements for 
the magnet helped in keeping this, as well as 
the lower temperature, satisfactorily constant. 

All rotation measurements were made for a 
field of approximately 13,300 oersteds which was 
used throughout the experiments. The rotations 
were determined for reversal of this field in terms 
of the galvanometer deflections produced, as 
already described. 

The cleaved crystal laminae varied in thickness 
from 0.4 to 0.9 mm; some eight different speci- 
mens were used altogether. Their diameter was 
only slightly larger than the 5-mm diaphragm 
hole over which they were mounted. The optic 
axis was, of course, perpendicular to the crystal 
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surface. Both right and left crystals were used, 
Since it might conceivably be urged that the 
anomaly found in the magnetic rotation was a 
spurious effect connected with the great anomaly 
in the optical rotation® at \ 1.16, some measure. 
ments were made on a combination of the two 
sorts of crystals which gave a net optical rotation 
of almost zero. All results on right- and left-hand 
crystals and on the combination gave magnetic 
rotations of the same sign and essentially the 
same magnitude. 


RESULTS AND DISCUSSION 


The results are shown in Figs. 2 and 3. Too 
much reliance should not be placed on the abso- 
lute values of the ordinates as almost the entire 
emphasis in this work was on the relative values, 
i.e., the dispersion of the rotation (i.e., of the 
Verdet constant), or form of the curves. A striking 
anomaly, somewhat resembling that found for 
the optical rotation® at 1.16, is seen to hold 
also for the magnetic rotation. (There is also a 
minor anomaly undoubtedly connected with the 
absorption band at 0.69u.) There is this im- 
portant difference, however, that while in the 
natural or optical rotation the temperature effect 
is almost independent of wave-length,’ in the 
present case the ratio of the two rotations is 
itself an ‘‘anomalous’” function of the wave- 
length, as shown in Fig. 3. 

In attempting to explain these results we at 
once find ourselves in difficulty because of a lack 
of knowledge of the refractive index, or rather 
indices, of this material in this spectral region. 
As a matter of fact we have already expended a 
great deal of effort in trying to clear up this 
matter but the results to date are not such as will 
justify their use. It will apparently be necessary 
to attack the problem from another angle. 

In the meantime, however, a fairly satisfactory 
qualitative explanation can be offered for the 
form of the curves of Figs. 2 and 3. If the 
explanation of the Faraday rotation is based on 
the inverse Zeeman effect and if it is further 
assumed that the intensities of the absorption 
bands are unequal, an anomaly® of the type of 

6 See reference 1, Ingersoll, Rudnick, Slack, and Under- 
wood, p. 1149. 

7 See reference 2, Rudnick and Ingersoll, p. 624. 


8See Wilhelm Schiitz, Wien-Harms Handbuch 16, 8 
(1936). 
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Fig. 2 results. To see if it would be possible to 
explain, on a reasonable basis, the form of the 
ratio curve as given in Fig. 3, we used the general 
equation for magnetic rotation,® written for a 
single term in the summation as 
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w being the frequency of the light. 
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Fic. 2. Dispersion of the magnetic rotation of crystal- 
line nickel sulphate at room temperature and low tem- 
perature. 


® Reference 8, p. 86. 
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Fic. 3. Solid line, observed ratio of the magnetic rotations 
of crystalline nickel sulphate at —125°C and 25°C. Dotted 


line, theoretical ratio calculated from Eq. (1). 


Lacking more exact information on the varia- 
tion of m, the terms preceding the bracket have 
been treated as constant. We then find that if 
we assume B/A = —0.01 and C/A = —5.00 and 
if we furthermore apply an over-all factor of 
about 3, there results a curve (see dotted line in 
Fig. 3) which has some qualitative resemblance 
to the one found by experiment. This is perhaps 
as far as it is profitable to go until more can be 
learned about the dispersion of the refractive 
index. 

We wish to acknowledge our indebtedness to 
Professor Philip Rudnick and to Messrs. Francis 
Davis, Edward Strait, W. J. Pearce, and others 
who have aided in this investigation, and par- 
ticularly to the Wisconsin Alumni Research 
Foundation for financial assistance. 
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The quantum theory of damping developed by two of us (Heitler and Peng) is applied to 
the production of mesons by proton-proton collisions. For this purpose the modification of the 
meson theory proposed by Mller and Rosenfeld is used. A primary radiation consisting of 
protons with a suitable energy spectrum is assumed, and it is shown that the rate of meson 
production is so high that nearly all mesons are produced in a top layer of the atmosphere of 
thickness 15-30 cm H;0. The variation‘of the meson intensity with energy, height, and geo- 
magnetic latitude is found to be in good agreement with the experiments. The transverse 
mesons, which have a very short lifetime, are seen to give, by decay, a satisfactory account 
of the soft component in the high atmosphere. A number of other effects (meson showers, 
transformation into neutretto’s) are discussed in Sections VI and VII. 





I. INTRODUCTION 


NTIL recently it has not been possible to 
apply Yukawa’s meson theory of the 
nuclear forces to cosmic-ray mesons and thus to 
establish the identity of the particles predicted 
by Yukawa with the cosmic-ray mesons. The 
reason for this deep-rooted difficulty is the 
following: The interaction between a meson and 
a nuclear particle is, in contrast to the electron- 
light interaction, a strong one, and becomes in- 
creasingly stronger at high energies. This makes 
a proper treatment of the reaction forces exerted 
by the meson field on the nuclear particles im- 
perative. However, as is well known, a treatment 
of the radiation reaction is intimately connected 
with the divergence difficulties occurring in every 
quantized field theory. To remove this difficulty 
two distinctly different sets of ideas have been 
put forward recently. Their difference can best be 
understood by remembering Lorentz’s expansion 
of the reaction force which a light wave emitted 
by an electron exerts on the electron. This 
reaction force can be expanded according to 
powers of the electronic radius r: The first term.is 
proportional to the acceleration and to 7 thus 
diverging for a point particle. This term is 
usually thought to be included in the inertia of 
the particle. The second term, the usual damping 
term, is independent of r and proportional to the 
time derivative of the acceleration. Higher terms 
are proportional to positive powers of r and are, 
as a rule, neglected. In the meson case it is the 
charge-and-spin degrees of freedom of the nuclear 
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particle which are coupled strongly with the 
meson field. We expect, therefore, that the reac- 
tion force will produce a twofold effect : (i) a large 
inertia to be attributed to these degrees of 
freedom, (ii) a large damping. In the first set of 
theories! mentioned above attention is concen- 
trated on the first effect. It is clear that in order 
to make the inertia term finite, a finite particle 
radius has to be introduced which makes a 
relativistic treatment of the nuclear particle so 
far impossible. In the second kind of theory no 
physical reality is attributed to the first term of 
the reaction force at all. The particle is strictly 
considered as a point particle. By suitable sub- 
traction the diverging inertia term of the reaction 
force is made to vanish (and so are the other 
diverging integrals occurring in the theory). The 
only finite part of the reaction force is then the 
damping term. Along this line Dirac’s? new 
quantum-electrodynamics (confined so far to the 
electromagnetic field) is based. Independently of 
Dirac but in the same spirit, though less general, 
two of us (Heitler and Peng)* have made an 
attempt at ‘“‘guessing”’ the correct equations of 
1 The so-called ‘‘strong coupling theory” put forward by 
Wentzel [Helv. Phys. Acta. 13, 269 (1940) ], Oppenheimer 
and Schwinger [Phys. Rev. 60, 150 (1941) ], and recently 
Pauli and Dancoff [Phys. Rev. 62, 85 (1942) ]. This theory 
is closely connected with former suggestions made 
Bhabha, Ma, and Heitler, in which it was assumed that a 


proton and neutron can exist in excited charge and spin 
states. 

2 Dirac, Proc. Roy. Soc. 180, 1 (1942) and lectures given 
at the Dublin Institute for Advanced Studies, to appear 
shortly in the Communication of the Dublin Institute for 
Advance Studies. 

a a _ and Peng, Proc. Camb. Phil. Soc. 38, 2% 
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quantum electrodynamics by omitting consist- 
ently and systematically all the diverging features 
from the quantized field theories. A simple new 
set of equations has been obtained which, how- 
ever, can so far only be used for non-static 
problems (transition probabilities) whilst the 
static field appears to be a problem of a higher 
order of difficulty. This paper*® will, ,in the 
following, be referred to as J. In contrast to the 
older expansion method used previously for the 
calculation of transition probabilities the new 
theory includes a damping force. In cases where a 
classical analogue exists it can be shown that the 
theory is equivalent to a treatment of the second 
classical damping term. The transition proba- 
bilities are to be calculated from a set of simul- 
taneous inhomogeneous integral equations. The 
application to multiple processes in J shows that 
the results are at least ‘‘reasonable.’’ We hope to 
be able to show that our theory can, at any rate 
as a good approximation, be derived from Dirac’s 
quantum-electrodynamics.‘ 

It should be possible to decide between the two 
kinds of theories experimentally. The cross sec- 
tion for the scattering of mesons by a proton 
depends upon the energy in a different way in the 
two theories. More precise measurements of the 
scattering cross section especially as a function of 
the energy are therefore very desirable. An ex- 
perimental decision would virtually decide 
whether the elementary particles have a finite 
radius or are point particles (cum grano salis). 

Further applications of our theory to the 
creation of mesons by light quanta and through 
proton-proton collisions are made in two papers,°® 
the results of which are the basis of the present 
work (referred to as JJ and III, respectively). 
They are summarized as far as they are needed in 
Section II. The aim of this paper is to show that 
the theory gives a satisfactory account—as far as 
we can see—of all the chief cosmic-ray phenomena 
connected with mesons, including their creation, 
their diffusion through the atmosphere, meson 
showers, and the transformation into neutrettos. 

* This does not mean that our results should be regarded 
as a proof for Dirac’s theory, even if our equations should 
prove to be in quantitative agreement with the experi- 
ments. Indeed many of the special features of Dirac’s 
a, only become apparent when static problems are 


*Hamilton and Peng; Heitler and Peng, to appear 
shortly in the Proceedings of the Royal Irish Academy. 
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We shall assume that mesons are created by a 
primary radiation consisting of protons, in ac- 
cordance with the geomagnetic evidence. 

It will be seen that the rate of production of 
mesons is so great that practically all the mesons 
are produced in a very thin top layer of the 
atmosphere. This is indeed what recent experi- 
ments® have shown to be the case. 

A certain ambiguity arises from the fact that at 
present it is not quite known which of the various 
modifications of the meson theory is to be used. 
As we endeavor to give a connected account of 
cosmic-ray and nuclear phenomena we choose that 
form of the meson theory which gives the best 
account of the nuclear forces: This is undoubtedly 
the form of the theory proposed by Mller and 
Rosenfeld.? Thus we assume that vector and 
pseudoscalar mesons exist (with equal coupling 
constants for transverse and pseudoscalar mesons) 
and that neutrettos (neutral mesons) also exist, 
whose coupling constants are half of those of 
charged mesons. 

In accordance with this theory and with the 
fact that only pseudoscalar mesons are found at 
sea level we assume that vector mesons have a 
much smaller lifetime than pseudoscalar mesons. 
A lifetime of 10-* sec. at rest suffices to explain 
that all vector mesons decay practically at the 
point where they are created. The meson pro- 
ducing layer of the atmosphere will. be seen to 
have a thickness of less than 15-30 cm H,0O. 
Most of the electrons arising from the decay of 
the transverse mesons are, therefore, also pro- 
duced in a very thin top layer of the atmosphere. 
We shall see in Section V that the number of these 
decay electrons and their energy spectrum is just 
of the right magnitude and type to produce, by 

cascade multiplication, the soft component in the 
high atmosphere. The observed intensity curve is 
in good agreement with the calculated one. Thus 
the soft component can be explained as tertiary 
products of the incoming primary protons. 


II. PRODUCTION OF MESONS BY PROTON- 
PROTON COLLISIONS. RANGE OF 
FAST PROTONS 


If a fast proton collides with another nuclear 
particle it may emit a meson in analogy to the 





6 Schein, Jesse, and Wollan, Phys. Rev. 59, 615 (1941). 
7 Maller and Rosenfeld, Kgl. Danske Vid. Sels. Math.- 
Fys. Medd. 17 (1940). 
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“bremsstrahlung” emitted when a fast charged 
particle is deflected in a Coulomb field. The rate 
of production has been calculated in’ JJI by 
applying the method of Weizsacker-Williams 
which has proved useful in the electromagnetic 
case: The field of the fast proton is expanded into 
a Fourier series and is, at sufficiently high energies, 
equivalent (i) to a beam of light quanta arising 
from the expansion of the Lorentz-transformed 
Coulomb field and (ii) to a beam of mesons 
arising from the Lorentz-transformed nuclear 
field. Both types of quanta interact with the 
nuclear particle at rest and can produce a meson, 
(i) by the process, 


hv+P—N+ Yt (1) 


and (ii) by scattering. It has been shown in JJI 
that the contribution of (i) is entirely negligible 
compared with that of (ii) except if the proton 
has an energy >10'* ev which is not of great 
interest in the present work. 

Thus the cross section for meson production is 
obtained by multiplying the number of mesons of 
a given energy occurring in the spectrum of the 
transformed nuclear field by the cross section of 
scattering of mesons by the nuclear particle at 
rest. 

Throughout this paper we shall use ‘‘natural 
meson’”’ units, putting 4=c=y (meson mass) = 1. 
Energies are then measured in units of yc? ~90 
Mev and cross sections in units of (h/yc)? 
=4.3X10-?* cm?. Furthermore, we introduce a 
unit for the thickness of matter traversed, 
namely, that thickness which a fast particle with 
unit charge has to traverse in order to lose an 
energy by ionization equal to wc?=1. This unit 
thickness is about 45 cm H,O. Thicknesses meas- 
ured in these units are denoted by x. The 
coupling constants of longitudinal, transverse, 
and pseudoscalar mesons with a nuclear particle 
are denoted by g?, f?, f’? (dimensionless), re- 
spectively, and their numerical values taken 
from the theory of nuclear forces (cf. reference 7) : 


g?=0.054, ft=f’?=0.13. (2) 


The coupling constants for neutrettos are halves 
of these values (fo? =f’? = 0.065, go?=0.027). For 
the meson mass we assume the value 75 of the 


proton mass. 
In the “equivalent spectrum” of the field of the 
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fast proton with energy E, charged and neutral 
mesons of all polarizations occur. The number of 
longitudinal mesons, however, is negligible. The 
number of transverse and pseudoscalar mesons 
with energy between ¢ and e+de was found in JJ] 
to be 


€)de d €)de d 

, P gy =—fDw, ogg =—fD,. (3) 
TE TE 

Here the D’s are rather complicated functions of 
¢/E, involving Hankel functions, but when 
plotted against «/E it turns out that they can 
quite well be represented by the following simple 
functions (with errors less than 30 percent) - 


Di =165, Dye =50(€/E)}, 
Du+Dys=200, 4Di+Dys=115. 


Only these combinations of the D’s will occur, 
D:, comprises both transverse polarizations. The 
number of neutrettos occurring in the spectrum is 
half of (3). 

The above method and therefore the expres- 
sions (3) are only accurate if: (i) €>1, (ii) EM 
(proton mass). For smaller proton energies the 
Weizsicker-Williams method fails completely, 
(For a more detailed discussion of the validity of 
the method cf. JIJ.) In the following we shall 
only be interested in meson energies e > 1/f<3 for 
which (i) is fairly well satisfied. We shall use 
(3) for energies E down to values ~M. This is 
certainly crude, but cannot involve very large 
errors because the slower protons are not very 
effective in producing mesons, on account of their 
small energy. If E< Ma proton may still produce 
mesons. A guidance as to the order of magnitude 
of this effect can be obtained in the following 
way: If the damping is neglected altogether and 
if EM the rate of meson production can quite 
easily be calculated directly by the old methods. 
The result can be compared with that obtained 
by using the Weizsacker-Williams method, also, 
of course, neglecting the damping, which means 
¢<1/f<3. The result is that the actual rate of 
meson production is about 10 times smaller than 
that obtained by using (3). We therefore expect 
that the rate of meson production drops rapidly for 


§ Cf. for example: Nordheim and Nordheim, Phys. Rev. 
54, 254 (1938). 
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E<M, and we shall neglect the contribution 
from protons with E<M. 

The cross section for scattering of a meson by 
a nuclear particle has also been derived in JJJ.° 
Since the exact expressions would be very com- 
plicated, we use only their asymptotic forms for 
é>1/f. « may then still be either <M (‘‘non- 
relativistic case’’) or > M (“‘extreme relativistic 
case”). Since a meson with given charge and 
polarization may be ‘‘scattered” into a meson 
with different charge (a charged meson may be 
transformed into a neutretto and vice versa) and 
different polarization, we write the result in the 
form of a matrix attributing rows to the primary 
and columns to the secondary particles. At- 
tributing rows and columns to the various 
polarizations as indicated in the formulae we 
found for the scattering cross sections: 


(1 +--+ + + +)long. 


transv. >charged 


pseud. 
sll (eKM) (5a) 


transv. neutral 





long. 


long. 


transv. charged 


167 


eM 


pseud. 


(e>M) (5b) 


pseud., 


} - | transv. neutral 








8 6 9» a 


The rows and columns marked transv. are to be 
understood as giving the transition probabilities 
into any one of the two transverse polarizations 
of the secondary meson and thus have actually to 
be understood as a submatrix: for instance 


1. 
(; ‘) instead of } in (Sb). At the places marked 
44 


by dots in (5) expressions occur which are of a 
smaller order of magnitude containing high 


* Earlier work assuming charged mesons only: Heitler, 
Proc. Camb. Phil. Soc. 37, 291 (1941), and Wilson, idid., 
301. (In these two papers only the non-relativistic case is 
treated.) Exact relativistic expression (for charged vector 
mesons only) by S. T. Ma, sbid., in the press. Cf. also 
Fierz, Helv. Phys. Acta. 14, 105 (1941); Landau, J. Phys. 
U.S. S. R. 2, 483 (1940). 


negative powers of e. We shall use (5a) and (5b) 
for energies e«< M and e2> M, respectively. 

In the rows and columns marked transv. of 
(Sb) an additional diagonal term has been’ 
omitted which increases with e but is multiplied 
by such a small factor that it becomes only 
appreciable for «21000. There are also other 
processes, not considered in this paper at all, 
which become important for such extremely high 
energies. Throughout this paper it must be kept 
in mind that our formulae are no longer valid 
whenever such high energies are involved. 

The most remarkable feature of (5) is the 
occurrence of selection rules. A meson has the 
tendency to conserve its charge as well as—to 
some extent—its polarization, and this becomes a 
strict law in the limit of high energies. The oc- 
currence of these selection rules is a special 
feature of our theory and entirely due to the 
damping. Nosuchselection rules occur if damping 
is neglected (which is permissible only for ¢ < 1/f) 
and our matrices (5) would be filled everywhere 
with expressions of the same order of magnitude. 
For the interpretation of cosmic rays it will be of 
particular importance that charged mesons can- 
not be transformed into neutrettos in a collision 
with a nuclear particle, except at very small 


energies (e~1/f) (cf. Section VII). 


Another remarkable fact is that the energy 
dependence changes from a 1/e* to a 1/e law for 
e>M. This fact will be fundamental for our 
understanding of cosmic rays. 

Multiplying (5) by (3) we obtain the cross 
section for the production of a meson of given 
energy and polarization. If «<M it is seen from 
(5a) that the pseudoscalar and the transverse 
parts of the equivalent spectrum combine to 
produce either transverse or pseudoscalar mesons. 
Thus the cross sections for the production of a 
meson with energy ¢ become 


8 de 
mel ~(Dirt+Dp), (6a) 
. (e<M) 


(6b) 


4 de 
®,,.(e)de = — f? —(Di+Dys). 
3. Ce 


Equation (6a) is already summed over both 
transverse polarizations. If «>M the pseudo- 
scalar part of the equivalent spectrum produces 
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only pseudoscalar mesons and the transverse part 
only transverse mesons. Thus 


de 
#1,(¢)de= (8/M) f* —Du,, (6c) 
(e>M) 
Pps(€)de= (16/M) f? —Dp.. (6d) 


In (6c, d) € is not the energy of the meson pro- 
duced but the energy lost by the moving proton, 
i.e., the energy of the meson produced plus the 
recoil energy of the nuclear particle originally at 
rest. The probability for the meson to take up an 
energy ¢’ leaving the energy e—e’ to the recoil 
particle is simply de’ /e. This is true if, as in fact is 
the case, the angular distribution of scattering is 
uniform in a Lorentz system where the meson and 
nuclear particle are colliding with opposite and 
equal momenta; thus the cross section for pro- 
ducing a meson of energy ¢’ is in the extreme 
relativistic case: 


BE q 
de! f $,,(-)—246,,(e)de’, (6'c) 
, € 


E de 
ae { ,.(€)—3,,(e’)de’. (for Ee’) (6'd) 
, € 


€ 


The energy distribution remains the same but the 
number of mesons is only 4 or 3? of what it would 
be if all the energy would be taken up by the 
meson. The rest of the energy is taken up by the 
recoil particles. Their energy distribution is also 
given by (6c, d) and their number is (6c, d) 
multiplied by $ and 3, respectively. These recoil 
particles are further capable of producing mesons 
and recoil particles. The process repeats itself 
until the energy has degenerated to a sufficiently 
low value to make further meson production 
impossible. The process very much resembles the 
cascade multiplication of electrons but is not so 
pronounced because the energy distribution (6) 
favors low energies more strongly than in the 
case of the bremsstrahlung emitted by a fast 
electron. Thus the energy degenerates more 
quickly than in the electromagnetic case. A de- 
tailed treatment of this cascade process lies 
outside the scope of this paper which is only 
intended to give a first orientation. We can take 
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account of it in a crude way by using (6c, d) 
instead of (6’c, d) also for the number of mesons 
produced. By doing so we represent the energy 
loss of the primary proton correctly. Also the 
total energy content of the charged mesons is 
correct because the energy given by the recoil 
particles to neutrettos is compensated by the 
energy of those charged mesons which are pro- 
duced by the recoil particles from neutrettos. 
Using (6c, d), we overrate the number of fast 
mesons by a factor 2 or 3, respectively, but we 
also underrate the number of slower mesons. 
Thus the error committed by using (6c, d) is a 
slight distortion of the energy spectrum in favor 
of high energies, whilst the total number of 
mesons will be somewhat bigger than what we 
obtain in this way. 

We multiply (6) by the number of nuclear 
particles contained in a cylinder of unit length, 
measured in our x units, and with cross section 
equal to our unit cross section. For water or air 
this figure is 1.18, but is not very different for 
other materials (for Pb it is about 1.3). 

We then obtain the number of mesons with 
energy ¢ produced by a proton in traveling the 
distance dx (using (2) and (4)): 





de ) 
1,(e)dedx = 82—dx (7a) 
t 
( (e<M), 
de 
Pps(€)dedx = 41—dx (7b) 
€ y 
de ) 
®,(€)dedx = 21—dx (7c) 
e 
» (e>M). 
e\'de 
®,,(€)\dedx = 12.3(— —dx (7d) 
E é y, 


For neutrettos the above figures are to be halved. 
All these formulae are only valid for «>1/fS3. 

From (7) it is seen immediately that a very 
fast proton produces more mesons with energy 
between 1/f and M than with e> M, but far the 
greater part of the energy is contained in the fast 
mesons. 

The mesons with e> M are all emitted in the 
forward direction, within a very small angle. 
This is not so for the mesons with e<M. The 
latter play a not very important role, except in 
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TABLE I. Range of fast proton. 








Eo 20 50 100 300 1000 = 10,000 
XB». M 0.18 0.45 0.83 2.0 4.1 28 

















the top part of the atmosphere. For the calcula- 
tion of intensities, the error will not be very large 
if we disregard any angular dependence and 
assume that all particles are always emitted in the 
forward direction. This we shall do throughout this 
paper. For this reason, and because we use in 
many places asymptotic laws instead of exact 
ones, and finally because we have only taken a 
crude account of the cascade process mentioned 
above, we must not expect too high an accuracy 
for our calculations. On the average, errors will 
be of the order of magnitude of, say, a factor 2. 
From (7) the energy loss of a proton can be 
obtained immediately. Multiplying by % to ac- 
count for the energy lost by producing also 
neutral particles and summing over-all polariza- 
tions we find (for E> M, of course) 
dE ad 
-—= 5 


dx polar. Wf 


eb(e)de=43 log O0.3E. (8) 


The energy loss is very high. It depends in a 
similar way on E to the ordinary energy loss by 
ionization but is roughly a hundred times larger. 
Per cm Pb (about {th of our x units) a proton 
with an energy of 310° ev loses an energy of 
2X 10° ev. 

Accordingly the range of a fast proton is very 
small. Of course, (8) is only valid for E>M; 
therefore, we can calculate only the distance a 
proton travels until it is slowed down to an 
energy ~ M. We find 


Fo dE 
XE, u= f az /(-—) 
M dx 


1 

ge eee (9) 
where ‘‘/i”’ is the integral-logarithm. For practi- 
cal purposes /ix can well be replaced by x/log x 
(this is exact for large x). For x£o,m we thus find 
for the range the values given in Table I. 
In our units the thickness of the atmosphere is 22. 
Thus a proton needs an energy of more than 7000, 
i.e., 7X10" ev in order to penetrate through the 
whole atmosphere and still retain an energy M. 
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The majority of the protons entering the atmos- 
phere at a latitude of 50° (E~22—50) lose the 
effective part of their energy in distances of 
0.2-0.5 (9-23 cm H,0). 

The theory does not tell us how quickly a 
proton loses its energy after having been slowed 
down to an energy M. Although meson produc- 
tion is then negligible compared with its rate at 
higher energies, the energy loss may still be much 
greater than that due to ionization. If our above 
estimate is correct (rate of meson production #5 
of that at higher energies), the energy loss would 
be 10 times greater than that due to ionization 
and the range of a proton with E~ M about 2 of 
our units (1 m H,0). 

Below we shall need the distance a proton 
travels in order to lose energy from E, to E: 


Eo 
E= - 
43 log 0.3E) 43 log 0.3E 





(10) 


in which we have replaced the /i(x) function by 
x/log x. 


Ill. PRODUCTION AND DIFFUSION OF PSEUDO- 
SCALAR MESONS IN THE ATMOSPHERE 


Since the primary protons have an extremely 
short range the majority of the mesons is pro- 
duced in a thin top layer of the atmosphere. The 
transverse mesons are expected to decay almost 
at once, and only the pseudoscalar mesons will 
travel through the atmosphere. Their absorption 
is due to two factors: (i) energy loss by ioniza- 
tion (ii) B-decay. The latter depends upon the 
distance the meson travels but not on the amount 
of matter traversed. We assume for simplicity 
that the density of the atmosphere at a depth x 
below the top is proportional to x. The probability 
of a meson at a depth x and with energy e« 
decaying while traveling the distance dx is there- 
fore dxb/ex, where } is inversely proportional to 
the lifetime + of a meson at rest. We choose r 
from the results of Rossi and Hall,"° who meas- 
ured the ratio r/p. For »=185, +r becomes 
2X7X10~* sec. 6 is then 


b=— X (distance in cm corresponding to one 
cr x unit at sea level) X (height of at- 
mosphere in x units) 


or 56=13. 
1 Rossi and Hall, Phys. Rev. 59, 223 (1941). 
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Let f(e,x)de be the number of mesons with 
energy ¢« at depth x and S(e, x)dx the number of 
mesons produced at depth x in the distance dx. 
Then f satisfies the diffusion equation 


—a———f+S, (11) 


This differs from the diffusion equation con- 
sidered and solved previously by Euler and 
Heisenberg" by the inclusion of the source func- 
tion S (which is determined by our theory). The 
term df/de accounts for the energy loss by 
ionization. Equation (11) can be solved by 
introducing new variables 


n=e+x,x 


instead of ¢€ and x. The solution with the bound- 
ary condition f=0 for x=0 is 


— b/n pz —¢ —b/n 
so, 2)=("—) J so.0(") dé, 
x 0 g 


n=e+x. (12) 


S is determined by our theory if we know the 
number of protons of each energy and at each 
depth x. The latter depends, of course, on the 
energy spectrum of the protons falling onto the 
top of the atmosphere. From the measurements 
at gteat depths it will be seen below that the 
primary spectrum is within certain comparatively 
wide energy regions a simple power law. For our 
purpose it is convenient to include also the 
logarithmic term occurring in (8) in the expres- 
sion for the primary spectrum; thus we assume 
that the number of primary protons with energy 
larger than E is of the form A | E/43 log 0.3E}~=. 
There is, of course, not the slightest reason why 
a should be exactly a constant. The experiments 
only show that a does not vary much if E changes 
by a factor 10 or so. Indeed, the measurements at 
extreme depths indicate an increase of a with 
increasing E. We shall determine a from under- 
ground measurements and shall find a= 2.2 for E 
between 100 and 1000. For larger E, a is bigger. 
Consequently we expect a to be smaller for 
E<100. The latter energy region will be of 
importance for the upper regions of the atmos- 


- "Euler and Heisenberg, Ergeb. d. exakt. Naturwiss. 17 
( 1938). 











phere and for the intensity curve of the soft 
component. a@ can be determined then from the 
shape of the Regener-Pfotzer-curve. We found 
a=1.3 satisfactory: for E<100. Thus the chief 
phenomena of cosmic radiation will be explained 
by the following crude but simple primary 
differential spectrum: 








a 
dF(E) =— (for E> 100), 

dE (E/43 log 0.3E)?2 

‘ (13) 
dF(E)=— (for E<100), 

dE (E/43 log 0.3E)'* 

2.2 £43 log 30\°* 

pettiasist! jaben imi, Oe TY 


143 100 


the connection bétween B and A being determined 
by the continuity of the primary spectrum. In 
reality a will change gradually from smaller to 
larger values as E increases. For any depth x 
larger than 2 only the high energy part of (13) 
will be important. The number of protons with 
energy larger than E at depth x is, according to 
(10) and (13): 


F(E, x)=A (——_+s) 
43 log 0.3E 


E —1.3 
F(E, x) -3(——__+s) 
43 log 0.3E 
100 is 
-3(. +s) 
43 log 30 


100 “63 
A(——_ +s) , (14) 
43 log 30 


according to whether 





E 100 
a nae 
43 log 0.3E 43 log 30 


The source function S(e, x) for pseudoscalar 
mesons is then, by (7b, d), 


1 
S(é, x)=41—F(M, x) (e<M), (15a) 
€ 


1 f’se\' OF(E, x) 
Sle, x) =12.3— f (=) ooo 
ev, \E OE 


(e>M). (15b) 

















3) 


ar 


b) 
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Fic. 1. Energy spectrum of pseudoscalar mesons at sea 
level, theoretical and experimental. Dotted curve: theo- 
retical spectrum at a height of 4000 m (50 cm Hg). 


Thus f(e, x) becomes: 


n—x b/n pz £o/9 
fle, x) -11("—) f aa 
x 0 


(n— or 


(n—E<M), (16a) 


fle, o=123() “far f af] 


g" A F(E, &) 
(n—£)2+" OE 











(n—§>M), (16b) 


n=et+x. 


In (16) the limits 0, x have to be replaced: in 
(16a) by n»n— M-- -x and in (16b) by 0---7—M if 
n—M lies between 0 and x. After the integration 
n is again to be replaced by e+x. In (16) the 
expression (14) for F(E, X) is to be inserted. 
Actually each integral (16) consists of two parts 
with two different values of a. It turns out, how- 
ever, that for most values of e, x, only one part of 
(14) is predominant, the other leading only to 
minor corrections; in particular for any depth 
x 22 only that part of F with a= 2.2 is important. 

The integrals (16) cannot be worked out 
exactly in closed form, but for most values of ¢, x 
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good approximations can be found. Equation 
(16a) is needed only for a few values of x and 
¢ («<M) and has been worked out numerically. 
In (16b) the logarithm occurring in F may be 
replaced by an average value. The chief contri- 
bution arises from small £ (corresponding to the 
fact that the protons have extremely small 
range). For any value of x>2, and even with 
reasonably good approximation for x=1, & can 
be neglected against 7 =e+x and the integration 
be extended to © instead of x. The exact condi- 
tion for this to be true is 


” 


>—__—_ (17) 
43 log 0.39 


which is satisfied for all but extremely small x or 
extremely large e. Equation (16b) then becomes, 
if only the high energy part of F, (a= 2.2) is used 





r(ii+é a—b/ b/9 
f(e, x) =12.3A CFO/IEC ”(*) es 


l(a) (a—}—/n) 
X (43 log 0.3n)2-!-""  n=e+x. (18) 


x 


For x <1 the contribution of (16b) has only been 
estimated but here (16a) is much more important 
than (16b). 

From (18) the asymptotic laws can immedi- 
ately be read off: 

(i) Tail end of the energy spectrum. Consider 
€>«x which for the lower part of the atmosphere 
also implies >b=13. 


fle, x)de« de-e-#(log 0.36)", (19) 


Since a= 2.2 the energy spectrum falls off a little 
less rapidly (because of the logarithmic term) 
than ¢?-?, 

(ii) Great depths. For underground measure- 
ments in dense materials the decay constant 0 
should be put equal to zero. f(e, x) is then a func- 
tion of e+x only. Integrating over e we find for 
the total number of mesons: 


I(x) = “fle, x)de x x~*(log 0.3x)*"". (20) 
us 


Thus the total intensity decreases like x~*, 
again apart from a logarithmic modification, and 
it is in this way that a= 2.2 was determined. 

In Fig. 1 the theoretical energy spectrum is 










































86 


compared with the measurements by Blackett!” 
for small and medium energies. It is seen that the 
shape of the theoretical curve is nearly the same 
as the experimental one. On the whole the 
theoretical spectrum falls off less rapidly than the 
experimental one. This was to be expected be- 
cause of the crude way in which we have ac- 
counted for the cascade production of mesons as 
explained in Section II. The experimental mini- 
mum at e= 20 does not appear in the theoretical 











Fic, 2. Theoretical number of pseudoscalar mesons and 
experimental number of penetrating particles as function 
of height. 


curve and is probably due to a special absorption 
process as was supposed by Blackett (production 
of proton pairs?). The maximum at e=10 is 
clearly indicated in the experiments as well as in 
the theoretical curve. It is difficult to trace the 
origin of this maximum as many factors combine 
to produce it. This maximum did not appear, 
however, in the calculations of Heisenberg and 
Euler, in which the exact source function was not 
used. For the low energy end (e~3) of f(e, x) also 
the contribution from those very few ‘very ener- 
getic protons which penetrate to the lower parts 
of the atmosphere and produce a large number of 
slow mesons is important. Hence the rise of the 
spectrum for small e. Actually the curve will bend 
down to zero for «<3, as indicated in the figure, 
but in order to reproduce this part of the curve 
correctly it would be necessary to use the exact 
cross section for scattering instead of the 
asymptotic one (5). 

There is another reason why we must expect to 

12 Blackett, Proc. Roy. Soc. 159, 1 (1937), 165, 11 (1938). 


Compare also: Jones, Rev. Mod. Phys. 11, 235 (1939); 
Hughes, Phys. Rev. 57, 592 (1940). 
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find few slow mesons for e <3, say. It will be seen 
in Section VII that a meson, once having reached 
an energy e <4, has a probability of 79 percent to 
be transformed into a neutretto before reaching 
the end of its range. This would remove a large 
fraction of mesons from the lower end of the 
spectrum. 

The high energy tail end of the spectrum, i.e., 
the decrease of f(e, x) <«*-? (modified by the 
logarithmic term) has already been adequately 
discussed by Euler and Heisenberg" and was 
found to agree with the measurements. In order 
to see how the energy spectrum varies with 
height we have also plotted the theoretical 
spectrum for x=15 (50 cm Hg). Slow mesons 
become relatively more predominant. 

In Fig. 2 the total number of mesons 


I(x) =f ft x)de 
us 


is plotted against height and compared with the 
recent measurements by Schein, Jesse, and 
Wollan." Here also the agreement is as good as 
can be expected. The maximum occurs at an 
extreme height of x=0.5 corresponding to only 
22 cm H;20 or 1.8 cm Hg. In the measurements at 
extreme heights presumably some of the primary 
protons also are included as ‘‘penetrating parti- 
cles,’’ but it is difficult to say how many protons 
have to be classed as “‘penetrating’’ because we 
do not know exactly the range of a proton with 
energy less than M. We estimate the range to be 
about 2 x units (=9 cm Pb) for E=M and ac- 
cordingly smaller for E<M. The theoretical 
curve of Fig. 2 refers to pseudoscalar mesons only. 
The two curves are normalized so as to agree at 
sea level. 


IV. LATITUDE EFFECT 


In the calculations of Section III no account 
has been taken yet of the fact that the primary 
energy spectrum is cut off for energies smaller 
than Ee depending on the geomagnetic latitude 
8. Es=22 for 8=50° and E»=150 for #=0° 
(equator). The correction can easily be made. 
From (16) the contribution from those protons 


which, at the top of the atmosphere had an 


18 Schein, Jesse, and Wollan, Phys. Rev. 59, 615 (1941). 
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energy < Ez has to be subtracted, i.e., for which 
t+E/43 log 0.3E < Es/43 log 0.3Es. 


It has turned out that for “normal latitudes’’ 
50°, say, the correction is negligible throughout 
the atmosphere. This means that no increase of 
intensity is to be expected for more northern lati- 
tudes, not only at sea level but also at any height 
above sea level, except the very top of the 
atmosphere, and even there the increase is small. 
This explains the well known knee of the latitude 
effect at 50°. 

The reason is easily to be understood. Protons 
stop producing mesons for E<M=10. 

Even the contribution from protons with 
energy between 10 and 22 is very small on ac- 
count of their smaller energy (in fact even 
smaller than assumed in this paper where the 
asymptotic law for meson production for E>10 
has been extended to E= 10). Thus the reason for 
the knee of the latitude effect is simply the fact 
that protons of energy less than Eso° cease to 
produce, or produce very few mesons. The same 
applies to the soft component. We shall see in 
Section V, that the soft component is probably 
due to the decay of transverse mesons in the high 
atmosphere, the former being also produced by 
the same primary protons. Thus the same knee of 
the latitude effect is to be expected for electrons, 
which is also found experimentally. The fact 
that the soft component shows the same latitude 
knee as the hard component supports the as- 
sumption that both are produced by the same 


TABLE II. Latitude effect. (Is0e— Jequ) /Ise. 








mesons soft component 
x 22 (sea level) 15 2 1 x=2 


theor. 0.13 0.13 0.18 0.33 0.62 
exp. 0.1 — — — 0.75 











primary radiation, i.e., by protons, as explained 
by our theory. 

For latitudes <50° an appreciable latitude 
effect is to be expected. We calculate the latitude 
defect Aof(e, x), i.e., the difference of f(€, x) between 
the North Pole and the latitude 3. As mentioned 
above Avf(e, x) is negligible for 8 > 50°. From (16) 
those contributions have to be subtracted for 


wae and Dymond, Proc. Roy. Soc. 171, 321 


which §+E/43 log 0.3E <Ee/43 log 0.3E». Thus 
Asf(e, x) is given by (for »—~>M): 


n—x\ 9 aby Foe 
safle, 12.3( ) f dt{ “aE 
0 n—& 


x 





(n—&)]}' &* OF(E, &) 
» (21) 
E J(y—&)*"* dE 
Es; Es 
+f=—____, 
43 log 0.3Ea; 43 log 0.3Es 
Es n 
"43 log 0.3E0 43 log 0.3 








5) 


n=e+x. 


(In the condition for é, — has been neglected 
compared with n). Asf is, of course, zero whenever 
in (21) an upper limit is smaller than the lower 
limit. A similar integral holds for 7»—& <M. 

We have worked out the integrals numerically. 
Here the change of a at E=100 is important. 
Integrating over e we find the latitude effect for 
the total number of mesons: 


Iw—Io= f Aof(e, x)de. 
/f 


The results are given in Table II. The latitude 
effect is fairly constant for low levels and only 
increases for x <2. The figure for sea level is in 
good agreement with the experiments, whilst at 
greater elevations no measurements seem to have 
been made yet for the hard component. For the 
figure concerning the soft component see 


Section V. 


V. THE TRANSVERSE MESONS AND 
THE SOFT COMPONENT 


In addition to pseudoscalar mesons the pri- 
mary protons produce a large number of trans- 
verse mesons which we assume to decay at once. 
If the mesons have energy €>1—which is nearly 
always the case—the decay electrons are emitted 
in the forward direction and take up with equal 
probability, any amount of energy between 0 and 
e. Thus a large number of slow and fast electrons 
are produced. The fast ones multiply by cascade 
multiplication and it will be seen that a satis- 
factory account of the soft component can be 
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given in this way. Let F(e, x) be the number of 
protons at depth x with energy larger than «. 
[F is given by (14). ] Then the number of elec- 
trons produced at x in the energy interval da is 
” de 
aa f F(e, x) ®,,-(e)—. 
3 € 
Let now C(a, &) be the cascade multiplication 
function (i.e., the number of electrons produced 


in a depth ¢ by one primary electron of energy 
3), the total number of electrons at depth x is 


Z.(x)= J ‘dt J - J , “ould 


X Fe, &)C(a,x—€). (22) 


Equation (22) can be worked out partly by 
making suitable approximations and partly nu- 
- merically. For #;,(e) the expressions (7a, c) and 
for F (14) are to be used. For the high atmosphere 
(x <4, say) practically only the low energy part 
of F(a=1.3) is important. For C we have used 
the figures given by Arley.'® 

% Arley, Proc. Roy. Soc. A168, 519 (1938). In the cascade 
theory different units for the thickness of matter traversed 


are usually used. For air or water one cascade unit is about 
§ of our x units. 
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To obtain an idea of the effect to be expected 
we consider large values of x. (x >5, say.) We can 
then make the following approximations: F(e, £) 
decreases rapidly with ¢ for any ¢e of importance. 
C(a, x—£) can be replaced by C(a, x) and the 
integral over — extended to infinity instead of to x. 
This amounts to assuming that all the mesons 
produced in the atmosphere are practically pro- 
duced at the very top and, decaying at once, 
produce a “primary” electron spectrum N(a)da 
say, which later produces cascade effects. N(3) is 
given by 


C-) de ) 
N(a)= f OO f F(e, dé. (23) 


We divide the range of energy 3 into 4 regions, 
and find for (a) (considering logarithmic terms 
as constants, and using the low energy part of F 
only) : 


da /43 log 0.3a\°* 
N(a)da = 33B— ——) (a> Es=22) 
3° a 
61 99 
rs aT: 
a’ 343 log 0.3a 





(Eo>a>M=10) 


142 
=Bds — +0.34} (10>3>1/f=3) 


a 
(24) 





=5.6Bda (a <3). 


Equation (24) plays the role of the primary 
electronic energy spectrum responsible for the 
soft component. N(a) decreases approximately 
like 1/a?, at any rate in the high energy region. 
But there are also a large number of soft electrons 
present. Apart from the latter N(a) is indeed 
very similar to the primary electron spectrum 
deduced first by Nordheim and Heitler'* from the 
intensity curve of the soft component on grounds 
of the cascade theory and later found by Bowen, 
Millikan, and Neher"’ by latitude measurements. 
We therefore expect that (24) should give rise to 
an intensity curve for the soft component very 
similar to the observed one. The effect of the low 
energy electrons with small range also present in 


16 Nordheim, Phys. Rev. 51, 1110 (1937); Heitler, Proc. 
Roy. Soc. 161, 261 (1937). 

17 Bowen, Millikan, and Neher, Phys. Rev. 52, 80 (1937) 
and 53, 855 (1938). 
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(24) will be that the intensity curve does not fall 
down very steeply near the top but becomes in 
fact rather flat. 

For the most interesting part of the atmosphere 
(x= 2, 3, maximum of the Regener-Hotzer-curve) 
the use of (24) is too crude an approximation. In 
fact quite a number of low energy mesons are 
still produced at levels x >2. We have therefore 
used (22) directly and evaluated the integrals. 
We find that on the average one primary proton 
produces about 1.5 electrons at x=2 or 3 in this 
way. 

In addition to the electrons produced by the 
decay of transverse mesons there are also elec- 
trons produced by the decaying pseudoscalar 
mesons. Their number is smaller but not negli- 
gible compared with the number of electrons due 
to transverse mesons. We easily find: 


tun faef anf <5 £)C(a,x—£), (25) 
0 0 a 2 


where f(e, —) is the energy spectrum of pseudo- 
scalar mesons as calculated in Section III. The 
result is for instance the following for x=2 and 
for one primary proton: Z;,=1.5, Z,.=0.6, 
Z ie+Zys= 2.1. Thus one primary proton produces 
in the average altogether about 2 electrons at the 
Regener-maximum, roughly 3 of which are due to 
pseudoscalar and 3 to transverse mesons. This is 
in very good agreement with the experiments. In 
Fig. 3 we have plotted the intensities of all the 
cosmic-ray components for the top part of the 
atmosphere (normalized for 100 incident protons). 
The agreement with the Regener-Hotzer-curve 
for the total intensity is as good as can be ex- 
pected. The position of the maximum is, however, 
slightly shifted towards greater heights, than is 
found experimentally. The reason is probably 
that we have been overrating the production of 
slow mesons (and therefore slow electrons) by 
protons of energy just above M. 

The latitude effect of the soft component is 
very large. We have calculated it for x=2 and 
included the result in Table III. The agreement 
with the measurements by Bowen, Millikan, and 
Neher"’ is good. The latter measurements refer to 
the total intensities, of which, however, the soft 
component is far the strongest. 

There is one experiment which is apparently in 
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contradiction to our assumption that the soft 
component originates from primary protons. 
Whilst the mesons in the lower atmosphere show 
an east-west asymmetry as big as the latitude 
effect proving that they originate from positive 
primaries no such east-west asymmetry has been 
found for the soft component in the high atmos- 
phere although the latter is also—and very 
strongly—latitude sensitive.'* If this result is 
taken at its face value it would mean that the 
soft component originates from equal numbers of 
positive and negative primary particles. We must 
remember, however, that the soft component is, 
if our picture is correct, created by the primary 
protons in a rather indirect way: The protons 
first produce mesons which decay into electrons, 
which then multiply by means of the cascade 
process. An east-west asymmetry can only be 
expected if in each of these three processes the 
particles are always strictly emitted in the 
forward direction with no appreciable angular 
straggling. This is surely the case if the secondary 
mesons have energies large compared with M, 
and if these mesons then decay while still having 
an energy large compared with one. 

In the high atmosphere, however, low energy 
particles are predominant. Indeed, as it appears 
from the calculations of this section, most of the 
electrons in the high atmosphere originate from 
mesons with energy less or not much greater than 
M. These mesons are bound to be emitted in all 
directions with only a slight favoring of forward 
directions. Even for the more energetic particles 
the predominance of forward directions will be 
impaired to some extent by the many processes 
which have to take place to produce the soft 
component. It might, therefore, well be that the 
east-west effect is completely masked by the 
large angular straggling of the soft component. 


TABLE III. Cosmic-ray intensities (100 primary particles) .* 








energetic protons 


soft comp. soft+hard mesons and neutrons 





x=2 x =2 sea level sea level 
theor. 208 220 1.7 1/500 
exp. — 225 4 1/2000 








*The column soft+hard (x =2) does not include the primary 
protons at this depth. Their number should be added to the theoretical 
figure 220 but is probably less than 20 (cf. Section VI). 


18 Johnson, Phys. Rev. 56, 219 (1939). 
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For a judgment of this theory it must be re- 
marked that, once the number of primary protons 
is normalized, and their energy spectrum given, 
the theory allows one to calculate the absolute 
intensities of all cosmic-ray components at each 
depth without further adjusting our intensities to 
those of the experiments. It is clear that the 
values of the intensities depend entirely on the 
constants occurring in the theory, especially on 
the nuclear coupling constants g, f. All these 
constants are either determined from nuclear 
facts (g, f) or by direct experiments (meson mass, 
decay constants). It is interesting to compare the 
intensities obtained by our theory with the 
experimental ones. Table III gives a selection 
of data. The absolute number of mesons at sea 
level is about twice of that calculated. This is 
quite what we had to expect as was remarked in 
Section II. (The last column refers to measure- 
ments by Jannossy, discussed in Section VI.) 


VI. PROTONS, NEUTRONS, MESON-SHOWERS 


The primary protons entering the atmosphere 
are very quickly slowed down to energy ~10. 
Since their energy loss is mainly due to producing 
charged mesons we must expect that half of 
these ‘‘protons’’ have in fact become neutrons 
after having travelled only a very small distance. 
In Fig. 3 we have plotted the number of these 
protons or neutrons with energy > M as a func- 
tion of depth for the primary energy spectrum 
(13), the number simply being given by the 
relation (14). It is seen that their intensity de- 
creases very rapidly. For a discussion of the 
experiments, however, we would rather like to 
know the number of protons (not neutrons) with 
all energies, including E< M. This number can 
only be guessed at present, since we do not 
know at what rate protons lose their energy 
once E is less than M. All we can say is that 
the number of protons drops very quickly to 
half the initial value (the other half have become 
neutrons). We might estimate that the energy 
loss of a proton becomes say 10 times smaller as 
E becomes less than M (cf. Section II1)—a figure 
which can only be a very crude guess. The in- 
tensity curve thus obtained is dotted in Fig. 3. 
The recoil protons and neutrons are not included 
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in these curves. The curve marked “‘total’’ does 
not include these protons except those at the 
very top of the atmosphere. 

After having reached an energy E<M the 
protons and neutrons have a much larger range 
than the distance that a proton or neutron can 
travel while having an energy > WM. It is there. 
fore not surprising that a large number of slow 
protons and neutrons are found in the upper 
atmosphere. Fluctuations of the range are large, 
since energy is lost in large portions, thus a small 
but appreciable fraction of them may even 
manage to travel through a considerable part of 
the atmosphere. In addition there are numerous 
fast and slow recoil protons. The experiments 
show a rapid increase of slow protons and 
neutrons with height, though their number igs 
still very small at a height of 4000 m compared 
for instance with the number of mesons. A quan- 
titative discussion is not yet possible. A certain 
fraction of these slow protons and neutrons are 
recoil products or are produced by photoelectric 
nuclear disintegrations, etc. 

In addition to these slow particles we must: 
expect, in any part of the lower atmosphere, a 
very small number of very energetic protons and 
neutrons, namely, the primaries themselves which 
have sufficient energy to travel such large dis- 
tances. Let us estimate their number, at sea level. 
The energy necessary to penetrate through the 
atmosphere is at least 7X10" ev. The primary 
energy spectrum (13) is valid for E up to 1000 
(10" ev). For higher energies we know from the 
measurements at great depths that a@ increases 
with increasing energy. Thus for the energies in 
question we may assume a= 3, for E> 1000, say. 
The number of primary protons at sea level is 
then given by A’x-*, x=22 where A’ is deter- 
mined by A [Egq. (13)] and the condition that 
the spectrum is continuous at E=1000, thus 
A’=23 if A, B, A’ are normalized for 100 inci- 
dent protons. Thus we find 1/500 energetic 
protons or neutrons at sea level for 100 primary 
protons at the top of the atmosphere. 

This figure may in fact be still too large, for 
the following reason: As was mentioned in 
Section II, our formula for the energy loss of a 
fast proton is no longer valid for energies 
>5 X10" ev, because then other modes of energy 
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TaBLE IV. Number of mesons (with «>10) produced by 
energetic proton or neutron. 
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loss become appreciable.'® Thus still higher 
energies may be required for a proton to pene- 
trate to sea level making their number still 
smaller. From measurements discussed below 
Janossy found that the number of protons or 
neutrons at sea level is about 1/12,000 of the 
total number of cosmic-ray particles at sea level. 
Putting the latter value to be about 6 (for 100 
incident protons) the experiments would give the 
number of these energetic protons and neutrons 
as 1/2000 which may be considered to be in 
reasonable agreement with our theoretical esti- 
mate 1/500. 

When these energetic particles travel through 
matter they will produce a number of mesons 
in succession thus producing a meson-shower. 
Showers consisting of penetrating particles have 
been observed by several authors.” The most ex- 
tensive measurements are due to Janossy: A radi- 
ation whose intensity was found to be 1/12,000th 
of the total cosmic-ray intensity at sea level was 
found to produce showers consisting of pene- 
trating particles which are certainly mesons with 
energy ¢€>10. The transition curve of these 
showers in lead reaches saturation at about 
5 cm Pb (about 1 x unit). About a third of the 
primary radiation producing these showers con- 
sists of neutral particles. Janossy found, how- 
ever, that for showers produced by the neutral 
primaries saturation is only reached after 10 cm 
of Pb. Each shower consists of an average 
number of 2-6 recorded penetrating particles. 
The actual number of particles in each shower 
may be larger. 

We can only give a very preliminary discussion 
of these experiments. The most obvious explana- 
tion is that the very energetic protons and 
neutrons expected from our theory are responsible 

% The limit of validity for the energy loss formula is 
higher than that for the scattering cross section because 
the former is an integral of the latter. 


*° Wataghin, de Souza Santos, and Pompeia, Phys. Rev. 
57, 61 and 339 (1940). Janossy, Proc. Roy. Soc. 179, 361 


(1942); Janossy, McCusker, and Rochester, Nature 148, 
660 (1941). 
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for these showers. There should be about equal 
numbers of protons and neutrons. Janossy found 
about a third to be neutral particles. Let E> M 
be the average energy of the protons and neu- 
trons and xz,y their range as given in Table I 
(i.e., the distance travelled until the energy is M). 

The number of mesons with energy >10 pro- 
duced is then according to (7): 


Eo dE E de 
eemeamear svn ce f —~2.1xBp, M; 
M —dE/dx 10 e 


N,,=21 


Nos ~ 15.5Eo-!xxp, M- 


Thus the total number of mesons produced is 
for various energies E» (from Table I) given by 
Table IV. ‘ 

These values cannot, however, be compared 
directly with the experiments. It has kindly 
been pointed out to us by Janossy, that the 
saturation point of the transition curve is not at 
all identical with the range of the primary radia- 
tion but marks the point where enough mesons 
are produced to be recorded as a shower. This 
figure is about 2-3. A further increase of the 
lead thickness will only increase the size of the 
showers but not the number of the showers. 
Experimentally, we find, therefore, that 2 or 3 
mesons are produced in the first 5 cm Pb or in 
the first x unit. This is in very good agreement 
with Eq. (26), the theoretical number of mesons 
produced in x= 1 being 2.1+15.5E,-!. We do not 
know E> exactly, but it certainly is rather large 
(>100). For greater thickness of Pb, Table IV 
shows that the size of the showers increases very 
much and even large showers are quite possible. 
On the other hand it is difficult to understand 
why neutrons should give rise to showers with a 
larger saturation thickness, than that for showers 
produced by protons. The processes considered 
in this paper lead to no explanation of this 
asymmetry. 

There are, however, a number of points, 
omitted so far in the present theory, which will 
have to be considered before a final discussion 
of the penetrating showers can be given: 

(i) As was mentioned in Section I, there are 
further modes of energy loss by a fast proton, 
expected from the present theory, including also 
a higher rate of production of transverse mesons, 
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which become important if the energy exceeds 
the value 5000, say. The primary protons would 
then hardly have a chance at all to reach sea 
level. The protons and neutrons responsible for 
the penetrating showers must be the fast recoil 
particles mentioned above which naturally have 
also smaller energies than the primaries. (If this 
argument holds the last column in Table III 
has no meaning.) Some of the processes taking 
place at extremely high energies depend on the 
Coulomb field of the proton which might perhaps 
indirectly account for the difference of the be- 
havior of the protons and neuti ‘s. 

(ii) Throughout this paper we have neglected 
the occurrence of multiple processes. It is indeed 
possible that a fast proton creates mesons, not 
only one by one in succession, but several mesons 
in one elementary act. We obtain the rate of 
occurrence of this event by multiplying the 
equivalent meson spectrum of a fast proton by 
the cross section for the splitting up of a meson 
into several mesons instead of by the scattering 
cross section (5). The splitting up processes have 
been calculated in J for the case when all energies 
concerned are smaller than M. In this case it 
has been found, indeed, that the splitting-up 
cross section for such multiple processes is always 
at least 10 times smaller than that for ordinary 
scattering (which is the reason why multiple 
processes were neglected in this paper). It may 
well be—but this could only be decided by 
further, rather difficult, calculations—that the 
ratio is more in favor of multiple processes if the 
energy of the proton is > M. If this should turn 
out to be true, the range of the fast protons would 
be diminished and the number of mesons in- 
creased. 

(iii) We have always assumed that the par- 
ticles in the nuclei of the matter traversed act 
independently. It might very well be that the 
nuclear fields of the nuclear particles overlap and 
interfere (as their Coulomb field does) producing 
an effect which increases with a different power 
of the atomic weight than the first power. This 
might lead to some changes of our theoretical 
results for heavy materials but hardly for air. 

(iv) So far we have considered only protons 
and neutrons as primary agents. It is to be 
expected that mesons also could create pene- 
trating showers themselves, by the very splitting 
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up process mentioned in (ii). Again no quantita- 
tive discussion is possible unless the calculations 
are performed for energies >M. A crude esti- 
mate on grounds of the non-relativistic calcu. 
lations shows indeed that the rate of occurrence 
of this process is of the observed order of 
magnitude. 

Mesons have also been found to occur in very 
large cascade showers. This is easily understood, 
According to our theory, a large cascade shower 
owes its origin to a very energetic primary proton 
emitting an energetic transverse meson which, 
by decaying, produces the very energetic elec. 
tron responsible for the cascade. Naturally, the 
primary proton will also produce numerous other 
mesons, including pseudoscalar ones, during its 
path through the atmosphere. These mesons will, 
of course, occur together with the cascade shower, 
In addition, energetic light quanta, have a small 
chance of creating mesons themselves by process 
(1) instead of producing an electron pair. The 
cross section for (1) was found in JI to be 
V2x°ef/é for pseudoscalar mesons, if «<M. Com- 
paring this with the cross section for pair pro- 
duction ~15Z?/137 we find that the chance for 
a light quantum producing a meson with e=10, 
say, is in air 1/200. Since in a big cascade shower 
at the point of its maximum several thousand 
light quanta may occur, quite a number mesons 
are to be expected. The mesons accompany the 
cascade shower but are much more slowly ab- 
sorbed; the ratio of the number of mesons to 
that of soft particles may therefore be quite 
appreciable at sea level. 


VII. NEUTRETTOS 


According to our theory neutrettos should be 
very frequent particles in cosmic radiation. The 
total number produced is half the total number of 
charged mesons produced. We know nothing, 
however, about their absorption and _ possible 
B-decay nor does it seem that they could produce 
any noticeable secondary effects. It is true, 
though, that in colliding with a nuclear particle 
a neutretto can be transformed into a charged 
meson. But, as mentioned in Section II, the 
cross section is very small at high energies de- 
creasing like e~* and it is unlikely that this effect 
will ever be observed. Slow neutrettos must be 
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very rare because there is nothing that would 
slow down a once fast neutretto. 

There is, however, one observable effect, by 
which the existence of neutrettos could be estab- 
lished experimentally, namely, the transforma- 
tion of a charged meson into a neutretto. At high 
energies €>1/f~3 again the cross section #° for 
this transformation is negligible (#°««~*) com- 
pared with the cross section for scattering ( « ¢~*). 
This is, indeed, an important result of our theory 
because it explains why all attempts at detecting 
this transformation™ have failed whereas the 
anomalous scattering of mesons has been ob- 
served. This is no longer true for energies of the 
order of magnitude e~3. The cross section be- 
comes then appreciable and, as we shall see, big 
enough to lead to observable effects. 

The method for calculating ° as a function of 
¢ is described in detail in JJ and the calculations 
are carried out up to a point where numerical 
evaluation is possible. The exact formula would 
be a rather lengthy and complicated expression 
because no approximation can be made in this 
case except that we can assume e to be small 
compared with M. We shall not give the formula 
here. & is, for very small energies, proportional 
to p?=e—1, rises to a maximum at e=2 and 
decreases rapidly for higher values. The value 
of # at the maximum (e=2) is in our units 


®rnax = 0.56. This refers to the transformation of 
a pseudoscalar charged meson into neutrettos of 
all polarizations (longitudinal, transverse, or 
pseudoscalar). 

We shall calculate the total probability for a 
meson to be transformed into a neutretto while 
traveling through matter, starting with a high 
initial energy «>1/f until it is stopped. If 
—de/dx is the energy loss by ionization this 
probability is 

A r® P%e) 
w=— -~———de, (27) 
2 1 — de/Ox 


where N is the number of nuclei per cm’, A the 
atomic weight. The factor } is due to the fact 
that a positive meson can only be transformed 
into a neutretto if it collides with a neutron and 


a negative meson only if it collides with a proton. 


"(Cf., for instance, Nishina and Birus, Sci. Pap, Inst. 
Phys. Chem. Res. Tokyo 38, 360 (1941). 
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The average number of ‘‘active’’ nuclear particles 
is therefore A/2. We have worked out (27) by 
numerical integration for Pb and found” 


w=0.79. (28) 


For any other materials w will not be very 
different because d¢/dx is nearly proportional to 
NZ which differs only slightly from NA/2. 
Practically the whole contribution to the in- 
tegral (27) arises from energies between «= #% 
and e=4. 

We see that, once a meson has reached an 
energy as small as 4, the chance is so great that 
it will be transformed into a neutretto. Only 21 
percent of the mesons reach the end of their 
range as charged mesons. It may very well be 
that this is one of the reasons why so very few 
slow mesons with e<2 occur in cosmic radiation. 

The length of path a meson has to travel in 
order to lose energy from e=4 to e=# is about 
1 km of normal air or 18 cm of Pb. While 
traveling through the gas of a cloud chamber the 
chance is so negligible that the effect could be 
observed. But, if mesons are stopped in a thick 
block of fead an appreciable fraction of the 
mesons should not reach the end of their range 
but should be transformed earlier into neutrettos. 
This is important for the type of experiments 
carried out by Rassetti.“ Rassetti found that if 
mesons are stopped in 10 cm of aluminum 
(equivalent to 2.8 cm of Pb) only a certain 
fraction (42+15 percent) give rise to decay- 
electrons. The effect was interpreted by Rassetti 
in the following way: Once a meson is stopped 
entirely and retains only an energy of a few 
electron volts it may either decay or else be 
captured by a nucleus giving its rest energy to 
the nucleus. Only negative mesons, however, 
can be captured in this way, since a positive 
meson with such a small energy could not come 
in contact with the nucleus. Calculations re- 
ferred to by Rassetti indicate indeed, that for 
such slow negative mesons the capture-proba- 
bility is larger than the decay probability. It may 
very well be, that transformation into neutrettos 
as well as nuclear capture is responsible for the 
fact that only a fraction of the mesons are really 


2 For this result the numerical values (2) for the coupling 


constants have been used. 
23 Rassetti, Phys. Rev. 60, 198 (1941). 
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decaying. In the above experiment 12 percent of 
the mesons should be transformed into neutrettos 
before being stopped. Of the remaining 88 per- 
cent half should be captured leaving 44 percent 
decaying into electrons. 

The experiments are not accurate enough to 
decide whether the number of decay electrons is 
actually appreciably less than 50 percent. The 
actual percentage should depend on the thickness 
of the absorber. 


VIII. CONCLUSIONS 


On the whole the theory was seen to give a 
satisfactory account of all the chief cosmic-ray 
phenomena. Although a number of points, es- 


pecially the cascade production of mesons (cf. , 


Section II), the east-west effect of the soft com- 
ponent and the penetrating showers, must await 
further investigation, it is probable that this 
theory is fundamentally correct. The conclusions 
that can be drawn from our results are twofold: 
(i) The theory of damping used as a basis for all 
our calculations will be a correct part of future 
quantum electrodynamics or at least a good 
approximation. (ii) Cosmic-ray mesons are in 
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fact identical with the quanta predicted by 
Yukawa which are responsible for the nuclear 
fields. It may perhaps be too early to say that 
the special form of the meson theory used in this 
paper (the one suggested by Mdller and Rosen. 
feld) is the correct form of the meson theory, 
We have carried out all the calculations also with 
a different form of the theory (assuming that 










only charged mesons exist) and found the agree. J 





ment with cosmic-ray experiments less good, 
though this form of the theory cannot be wholly 
excluded. (The energy loss of fast protons would 
be three times smaller.) It is satisfactory that 
the form of the meson theory which gives the 
best account of the nuclear forces also agrees best 
with cosmic-ray facts. Especially, the assumption 
of both pseudoscalar and transverse mesons is 
strongly supported by the fact that all com- 
ponents of cosmic radiation can be explained as 
owing their origin to one kind of primary par- 
ticles, i.e., protons. 

We are very much indebted to Dr. L. Janossy 
for a most useful comment on this paper and also 
for communicating some of his experimental 
results to us before publication. 
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PART I: DESCRIPTIVE ANALYSIS 
I. Introduction 


NE of the greatest difficulties encountered 
in the production of low temperatures by 
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On the Adiabatic Demagnetization of Iron Alum 
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The behavior of iron alum during adiabatic demagnetization to temperatures near the 
absolute zero is of particular theoretical interest because of the deviations from Curie’s law 
that arise from the perturbing actions of crystalline field and magnetic dipole-dipole coupling. 
The effect of these perturbations on the magnetic moment and the entropy is calculated 
exactly to second-order terms in the magnetic coupling and to third-order terms in the erystal- 
line potential. The calculations are presented for crystalline fields of either cubic or axial 
symmetry and are valid for the case of large applied fields in which saturation effects become 
important. Theoretical values of the adiabatic moment are found to be in satisfactory agreement 
with the experimental values determined by Casimir and de Haas. A true thermodynamic 
scale is established that enables the temperature to be calculated at any value of the magnetic 
field during demagnetization. The relationship of this scale to the temperatures determined by 
the magnetic method of de Haas and Wiersma is discussed. 
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the method of adiabatic demagnetization is the 
determination of the thermodynamic tempera- 
ture. Until recently, it has been the general 
practice of experimenters to set up a ‘‘magnetic” 
temperature scale based on the validity of Curie’s 
law. Deviations from this law, however, will be 
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produced both by the magnetic interaction be- 
tween the paramagnetic ions of the alum and by 
the Stark splitting due to the crystalline electric 
field. These effects become progressively more 
important as the temperature is lowered so that 
there may exist an appreciable difference between 
the ‘magnetic’ or “reciprocal susceptibility” 
temperature 7* and the true Kelvin tempera- 
ture 7. 

To overcome this uncertainty in temperature 
measurement, various procedures have been pro- 
posed, most of which base their determination of 
the true thermodynamic temperature on calori- 
metric measurements. De Haas and Wiersma,! 
however, have proposed a method of determining 
the Kelvin temperature entirely from magnetic 
measurements. Their procedure is based es- 
sentially on the experimental fact that when 
certain alums are subject to an adiabatic de- 
magnetization process, the magnetic moment is 
found to remain constant, within the limits of 
experimental error, as the field is decreased from 
its initial high value to some much lower value. 
When this is the case, one can show purely from 
the laws of thermodynamics that the true temper- 
ature T is a function of the field strength H and 
is given by an expression of the form 


T—T,)=c(H—A)), (1) 


where c is a parameter depending on the value of 
the initial magnetization, and TJ») and Hp» denote 
values of temperature and field, respectively, at 
the beginning of demagnetization. This equation 
was applied by de Haas and Wiersma to the 
particular case of caesium titanium alum as this 
material experimentally showed no change in the 
value of the magnetic moment in the range from 
24,000 gauss to 100 gauss. 

Similar experiments on a powdered specimen of 
iron alum by Casimir and de Haas? indicate that 
the magnetic moment during adiabatic demag- 
netization does not remain constant over any 
appreciable interval, but varies with the applied 
field strength H7, in the following manner: 


M = M)(1—6/H.’). (2) 


(1935) J. de Haas and E. C. Wiersma, Physica 3, 491 

?H. B. G. Casimir and W. J. de Haas, Physica 7, 72 
— The symbol capital B in this paper is the same as our 
small b. 
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Here M, is the original magnetic moment due to 
an initial magnetic field Ho of 18,100 gauss, and b 
is aconstant whose experimental value is 1.7 X 105. 
From this value of 5 it is readily seen from 
Eq. (2) that for sufficiently large values of the 
field strength (H,>10,000) the term in 1/H,? is 
negligible compared to unity. Hence for values of 
the field strength above this magnitude the 
method of de Haas and Wiersma for determining 
the true Kelvin temperature is clearly applicable. 
For fields between 10,000 gauss and 5000 gauss, 
the existence of the second term in Eq. (2) above 
will cause a small error to be made in the 
temperature determination, while for fields less 
than 5000 gauss the error becomes appreciable. 

It is not surprising that the magnetic moment 
of iron alum fails to remain constant during the 
process of adiabatic demagnetization. Rather the 
reverse would be the case as even a simple 
theoretical treatment shows the existence of a 
1/H* term in the expression for the adiabatic 
moment. Such a treatment has been given by 
Casimir and de Haas? who show that for low 
values of the initial field, the coefficient of the 
1/H? term [i.e., b in Eq. (2) ] is some function of 
the original magnetization and hence of the 
original field strength, but does not depend on 
the value of the final field. They were unable to 
evaluate this function for the general case, but 
for the special case of low initial fields that do not 
produce saturation, they found b to be 1.3 X 10° as 
compared to the experimental value of 1.7 «10°. 
The discrepancy between the computed and the 
experimental value is probably due to neglect of 
saturation effects. The magnetic fields used in 
their experiment were quite large (18,000 gauss) 
and for fields of this order of magnitude one 
would normally expect saturation effects—neg- 
lected in their treatment of the problem—to be 
very important. That such is the case will be 
demonstrated in this article for, contrary to the 
opinion of Casimir and de Haas, it ts possible to 
make an exact theoretical determination of the 
adiabatic moment for high fields. The method of 
solution is described in Part I of this paper but 
the mathematical details associated with the 
determination are placed in Part IT. 

Once the necessary theoretical treatment has 
been developed, it can be used not only for the 
determination of the adiabatic magnetic moment 
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a. Cubic field 


Fic. 1. Splitting of J=5/2 level by a crystalline field. In notation of Penney and Schlapp 
af 6.=3y-=72D¢q, 54=3y.=18Aa. 


(reference 


but also becomes immediately available for use in 
two other connections. First, it enables us to 
calculate the true thermodynamic temperature 
and to estimate the error that will arise if the 
temperature is determined by the simple mag- 
netic method of de Haas and Wiersma, i.e., by 
use of Eq. (1). Secondly, by comparing the 
theoretical expression of M with the experimental, 
we obtain information in regard to the nature of 
the crystalline field that acts on the paramagnetic 
ions of the alum. 


II. Method of Approach 


The procedure adopted is based upon a method 
of attack originally developed by Van Vleck? and 
used by him to study theoretically the adiabatic 
demagnetization of caesium titanium alum 
(CsTi(SO,4)2-12H2O). The indicated mathemat- 
ical treatment is applicable for the case of iron 
alum (NH,Fe(SO,)2:12H:O) as well, but there 
are some significant differences in treatment due 
to basic differences in the structure of the two 
materials. In caesium alum the Ti+** ion has a 
spin S=} and lies in a doubly degenerate ground 
state and hence, according to a fundamental 
theorem of Kramers,‘ is unaffected by a sur- 
rounding crystalline field of any type of sym- 
metry. The free orientation of the paramagnetic 
ions is thereby hindered only through the agency 
of the magnetic coupling, and possibly a little 
exchange coupling, between spins. In iron alum, 
on the other hand, the paramagnetic Fe+*+* ion 
has a spin S=5/2 and lies in a sixfold degenerate 
ground state. This sixfold level is subject to 
Stark splitting by the electric fields of the crystal 


3J. H. Van Vleck, J. Chem. Phys. 6, 81 (1938). One 
typographical error should be noted. In the expression 
for G, on p. 83 the sign before the (g8H/kT) term should be 
reversed. 

*H. A. Kramers, Proc. Kon. Akad. Amst. 33, 953 (1930) ; 
34, 965 (1931). 


therefore be written 
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b. Axial field 








as has been shown by Bethe,® Kramers,‘ and Van 
Vleck and Penney.® A field of axial symmetry 
splits the sixfold degenerate level of the ground 
state into three separate doubly degenerate 
levels, while a crystalline field of cubic symmetry 
separates the sixfold level into a fourfold and a 
twofold one. This is illustrated in Fig. 1 (see 
Penney and Schlapp’). 










Hamiltonian Function 





The first step in our theoretical treatment must 
be to set up an expression for the Hamiltonian 
that will correctly represent the conditions to 
which the paramagnetic alum is subject during 
the adiabatic demagnetization process. In addi- 
tion to the applied magnetic field acting on the 
paramagnetic ions of the alum, there is a mag- 
netic coupling of the three-dimensional array of 
dipoles and an electrostatic interaction between 
the various paramagnetic ions and their sur- 
rounding crystalline fields. The Hamiltonian may 















KR=-—gBH >; Sz.+Liri Wit Li Vi, (3) 


where g is the Lande splitting factor, 8 the Bohr 
magneton, and S the spin of the Fe*** ion. The 
first term thus represents the Zeeman energy due 
to a field H supposedly directed along the z axis. 
The second and third terms respectively represent 
the magnetic dipole-dipole coupling and _ the 
crystalline field interaction. The exchange inter- 
action between the paramagnetic ions has been 
neglected in Eq. (3) because, in iron alum, the 
ions are known to be spaced relatively far apart. 






















5H. A. Bethe, Ann. d. Physik 3, 133 (1929). 

*J. H. Van Vleck and W. G. Penney, Phil. Mag. 17, 
961 (1934). 

7™W. G. Penney and R. Schlapp, Phys. Rev. 41, 1% 
(1932) ; 42, 666 (1932). 
(1955) A. Beevers and H. Lipson, Proc. Roy. Soc. 148, 664 
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In the mathematical treatment that follows, the 
Zeeman term of Eq. (3) is taken as the unper- 
turbed-part of the Hamiltonian, and the other 
two terms are considered as small perturbations 
causing slight shifts in the unperturbed energy 
levels. 

The first of the perturbing influences, namely 
the magnetic dipole-dipole coupling, is explicitly 
introduced into the analysis by writing for w;; the 
expression® 


wis= 20 { (S¢+ $3) —3(S;- 743) ($;- Tis) ri | 


where s; is the spin angular momentum vector 
in units of h/2x of atom i, r;; is the radius vector 
connecting atoms 7 and j, and r;; is its modulus. 
As the lattice spacing is known, there are no 
unknown constants in this expression, and hence 
the computation of the influence of this term, 
though lengthy, is straightforward. 

On the other hand, when the crystalline field 
potential is considered, either one or more un- 
known constants are introduced depending on 
the nature of the crystalline field assumed. From 
x-ray analysis by Beevers and Lipson,® it is 
known that the field surrounding the paramag- 
netic ion possesses trigonal symmetry. The most 
general potential of this type is'® 


Vi= Fir) T+ Fi(r) T+ Fi (r)(TE+T), (4) 


where the F’s are coefficients depending only on 
r, the 7;"’s are spherical harmonics of the form 
T,\"=P,"(cos 6) exp (img), and the P,"’s are 
Legendre polynomials. In these expressions @=0 
corresponds to the trigonal axis, which coincides 
with one of the four-body diagonals of the cubic 
lattice. If the fourth-order terms of Eq. (4) are 
neglected then the crystalline field is one of axial 
symmetry and that part independent of r can be 
expressed by 


Vi=A (x+y? — 22’), (4a) 


where A is an arbitrary constant associated with 
the strength of the field, and the primed coordi- 
nates are rectangular coordinates in which the 
trigonal axis coincides with the z’ axis. This is the 
only type of second-order term compatible with 
trigonal symmetry. Probably the most important 


* J. H. Van Vleck, J. Chem. Phys. 5, 320 (1937). 
A. Siegert, Physica 3, 85 (1936). 
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fourth-order terms in (4) are those having cubic 
symmetry. In terms of rectangular coordinates 
referred to the cubic axes, these terms are given 
by 

Ve=D(xt+y* +24 — fr’), (4b) 


where D is an arbitrary constant measuring the 
strength of the field, and the 2r* term is included 
in order to have V, satisfy La Place’s equation. 

The Hamiltonian of Eq. (3) is used as the basis 
for a rigorous calculation of the partition function 
Z. Once this is established, the magnetic moment 
and the entropy, both of which are required in 
order to get the adiabatic moment, are obtained 
by simple differentiation. In the calculation of Z, 
we assume the dipole-dipole coupling and the 
crystalline field potential are both small com- 
pared to the energy of thermal agitation. As the 
partition function is expanded in powers of the 
ratio of the perturbing influences to kT, the 
convergence of this expansion will be rapid as 
long as our original assumption is valid. Fortu- 
nately, the ratio is small except at sufficiently low 
temperatures where the interactions become 
comparable to kT. The unperturbed energies 
arising from the first term in the Hamiltonian are 
not assumed small compared to kT as it is the 
form of M for large values of H that is desired. 

The perturbing influences are taken explicitly 
into account by using the well-known quan- 
tum-mechanical expressions from perturbation 
theory," and then calculating the matrix ele- 
ments of the perturbing part of the Hamiltonian. 
The calculation of the matrix elements and also 
of the partition function is greatly facilitated by 
use of the methods of group theory. These 
methods are particularly applicable to our 
problem where the experimental use of a pow- 
dered specimen instead of a single crystal 
necessitates the writing of the crystalline field 
and dipole-dipole contribution to the partition 
function in terms of generalized coordinates over 
which an average must be taken before the final 
expressions for the magnetic moment and the 
entropy can be obtained. It is precisely in the 
carrying out of this averaging process—otherwise 
a long-winded and tedious job—that the group 
theory method affords such a great saving. 


"E. C. Kemble, Fundamental Principles of Quantum 


Mechanics (McGraw-Hill, New York, 1938). 
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III. Adiabatic Moment to Second Order 
in Perturbing Influences 


The development of the partition function is 
carried out in Sections 6 and 8 of Part II. The 
expressions for the magnetic moment and the 
entropy, given by Eqs. (21) and (22) of Section 9, 
are functions of both T and H. To find the 
adiabatic magnetic moment, the procedure is to 
eliminate T between the two equations, and then 
since S is constant during the course of an 
adiabatic demagnetization experiment, one has 
an equation for the adiabatic magnetic moment 
that depends only on the field strength. In 
Section 10, this expression is shown to be 


M=M)(1—b/H’), (5) 


where 6 represents a function of the initial field 
and initial magnetization given by Eq. (29) of 
Section 10. Equation (5) is of the same form as 
the experimental relation as regards dependence 
of M upon H., It should be noted, however, that 
in the theoretical expression H is to be inter- 
preted as the true field acting on the paramag- 
netic ions and not as the experimentally applied 
field. 

To find the numerical value of the theoretically 
derived expression for b the unknown parameter 
or parameters that entered into the crystalline 
field potentials must be evaluated. This evalua- 
tion is made by comparing the theoretical ex- 
pression for the specific heat in zero applied field 
with the experimental. Both are of the form 
C,=A'’/T?, where the coefficient A’ includes the 
contribution of both the dipole-dipole coupling 
and the crystalline field interaction. If the 
trigonal field acting on the paramagnetic ion is 
taken as purely cubic or purely axial in character, 
then only a single unknown parameter enters 
into the expression for A’ and the constant is 
directly evaluated by comparison with the known 
experimental value of this constant for iron alum. 
This evaluation is made in Section 10 of Part IT. 
The experimental value of A’ has been de- 
termined by several investigators, some of whom 
used entirely different methods of attack. Con- 
sidering this, it is remarkable how well they 
agree.” 

2H. B. G. Casimir, W. J. de Haas, and D. de Klerk, 
Physica 6, 241 (1939); F. K. du Pre, Physica 7, 79 (1940); 


P. Teunissen and C. J. Gorter, Physica 7, 33 (1940); C. 
Starr, Phys. Rev. 60, 241 (1941); F. Simon, private corre- 
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After these unknown crystalline field constants 
have been evaluated, one can calculate the 
numerical value of the coefficient of 1/H? and 
thereby obtain the final expression for M correct 
to the second order in the two perturbing influ. 
ences. In Section 10, this expression is shown to 
be 

M=M,(1-—9.4X10'/H?). (Sa) 


The numerical value of 9.4 is calculated on the 
assumption that the trigonal field acting on the 
paramagnetic ion is predominantly cubic. This 
value should be changed to 11.5 if an axial field 
is assumed. In either case, however, a comparison 
with the experimental values of Casimir and 
de Haas shows that the theoretical values of 4 
obtained from (5a) are too high. Part of this 
discrepancy is accounted for by the fact that, 
because of demagnetization corrections, the true 
field H appearing in Eq. (5a) is somewhat less 
than the applied field 7, entering into Eq. (2). It 
is shown in Section 10 of Part II that this 
accounts for only a minor part of the difference 
between the theoretical values for M and the 
experimental and leaves the major part of the 
discrepancy still to be explained. 

Two possible reasons for the existence of the 
discrepancy present themselves. The first of these 
is the neglect of third and higher order terms in 
the expansion of the partition function and the 
second is the restriction that the crystalline field 
acting on the paramagnetic Fe*** ion is trigonal 
in character. Of these two, the former seems to be 
the more likely as the trigonal nature of the field 
is attested by the experimental x-ray methods of 
Beevers and Lipson.® Still, from a theoretical 
standpoint, the possibility of the existence of a 
more general field should not be excluded and 
hence we shall first consider how the removal of 
the trigonal restriction would affect the results, 
spondence to J. H. Van Vleck (1940). Casimir, de Haas, 
and de Klerk determined the specific heat constant by 
direct magnetic measurements on pure crystals. They find 
A’ to be 2200+40. Du Pre records a value of 2300450 
determined by experiments on paramagnetic relaxation. 
Indirect measurements by Teunissen and Gorter, from 
a ye on paramagnetic dispersion, tend to lower this 
value slightly while similar indirect measurements by Starr 
tend to raise it by about the same amount. Private corre- 
spondence from Simon indicates a value of the splitting 
distance of 0.193 determined from direct magnetic measure- 
ments. This corresponds to a specific heat constant o 
2350. All values whether determined by direct or indirect 


measurements thus approximate 2300 and this we take as 
the correct magnitude. 
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TaBLE I, Comparison of theoretical and experimental 
values of adiabatic moment of iron alum. 


DEMAGNETIZATION OF 








M(calc.) M calc.) 








Cubic field Axial field 

He H Mi(exp.) 2ndord. 3rdord. 2nd ord. 
18,100 18,000 -~ 4.70 4.70 4.70 
11170 11,040 4.77 4.70 4.70 4.70 
5'550 5,420 4.66 4.69 4.68 4.68 
3,130 3,000 4.61 4.65 4.64 4.64 
1,352 1,230 4.25 4.41 4.33 4.34 
912 807 3.7 4.02 3.75 3.87 
494 412 2.9 2.10 0.10 1.53 














Note: Values of M expressed in units of N@. 


Failure of Arbitrary Crystalline Field to 
Improve Agreement 


Let us replace the trigonal expression for the 
crystalline field given by Eq. (4) by the more 
general expression 


4 2 
Vo= Lo aul st+ LD beT2’, (6) 


s=—4 r=—2 


where the a’s and b’s are arbitrary coefficients. 
If the coefficients of this general field are so 
selected as to retain unrestricted trigonal sym- 
metry, it should be noted that (6) still represents 
an increase in generality over the cubic and axial 
fields thus far considered. It would be useless to 
include terms in 7; and 7; as matrix elements for 
odd order potentials must vanish in virtue of the 
LaPorte rule. The matrix elements have already 
been computed for the 7, and T; term and hence 
it is an easy matter to set up the partition func- 
tion for any combination of these terms with 
arbitrary coefficients. Thecontribution of this type 
of field to the specific heat and to the adiabatic 
magnetic moment is given by Eqs. (33) and (32), 
respectively, of Part II. Examination of these 
two equations discloses that, at least as far as the 
second-order terms are concerned, mo assigned 
values of the arbitrary parameters can bring the 
theoretical and experimental results in closer 
agreement than already found under the assumption 
of a trigonal field of axial symmetry. Hence the 
failure of the theoretical results to give closer 
agreement with the experimental data is defi- 
nitely not due to any restriction on the nature of 
the crystalline field. 


IV. Beneficial Effect of Higher Order Terms 


The most likely reason for the discrepancy 
between the theoretical values of Mf and the 
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experimental values is the neglect of third and 
higher order terms in the expansion process 
employed in the calculation of the partition 
function. Fortunately, it is possible to save much 
of the labor associated with the calculation of 
higher order terms by making use of some previ- 
ous theoretical work on iron alum by Kronig and 
Bouwkamp." For large fields and for arbitrary 
direction of the magnetic field relative to the 
crystal axes—just the required conditions neces- 
sarily imposed on us by experiment—they have 
calculated the energy levels of iron alum to the 
third order in the crystalline field parameter. 
Their calculations were made for a crystalline 
field of cubic symmetry and do not take into 
account the effect of the magnetic dipole-dipole 
coupling. This is not a serious drawback, however, 
as the neglect of the third-order term in the 
dipole-dipole coupling will not appreciably 
change the final result. This is so because the 
contribution of the crystalline field terms to the 
magnetic moment is about eight times as large as 
the dipole-dipole contribution. 

Once the energy levels are known, it is a 
straightforward process to determine the partition 
function, the entropy, and the magnetic moment. 
The details of this treatment are given in Sec- 
tion 11. The effect of inclusion of the third-order 
crystalline field terms is to add a new term in 
1/H® to the adiabatic moment. The complete 
expression for the magnetization then becomes 


M=M,(1—bH~—b'H-), (7) 


where b’ is independent of the final field but 
depends upon the initial conditions. For the ex- 
perimental conditions of de Haas and Wiersma b’ 
takes the numerical value of 30X10* while } 
itself has the value previously stated, vis., 
9.4 10*. If M is expressed in the form of Eq. (2), 
the correction to 6 due to introduction of the 
third-order term in the cubic crystalline field 
interaction is a factor (1+320/H). This factor 
becomes increasingly important as the field is 
lowered. The effect of the additional term in the 
expansion has thus been to increase the value of 
the coefficient of the H-* term and hence its 
inclusion has tended to lower the value of the 
adiabatic magnetic moment in the proper direc- 


13R. de L. Kronig and C. J. Bouwkamp, Physica 6, 290 
(1939). 
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tion in order to give better agreement with 
experimental data. 

The theoretical values of M based on Eq. (7) 
and Eq. (5) are compared with the experimental 
data in Table I. From this table it is readily seen 
that the agreement between the third-order 
theoretical and the experimental values of M is 
now quite satisfactory. One exception is the last 
row in the table, but here the true field strength 
is less than 500 gauss and for a field as low as 
that, the temperature obtained by assuming the 
law of de Haas and Wiersma to apply is already 
below the characteristic magnetic coupling tem- 
perature. When this is the case, the convergence 
of the expansion breaks down completely. 

The relation between the magnetic moment 
and the final field strength is graphically shown 
by Fig. 2, where the actual experimental points 
are indicated by the small crosses. No experi- 
mental value of M is recorded by de Haas and 
Wiersma in the neighborhood of the original 
field, and no mention is made of the order of 
accuracy of their experimental data. We are no 
doubt safe in assuming that the experimental 
error is of the order of 3 percent. For all values of 
the field above 900 gauss the greatest difference 
between the theoretical and the experimental 
value of M is less than this amount. Hence it does 
not appear worth while at present to carry out 
any further calculations on the effect of still 
higher order terms. 


Nature of Crystalline Field Present in Iron Alum 


It was thought likely when the theoretical 
calculations were first started that a comparison 
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between the experimental values of the adiabatic 
moment and the theoretical values, as computed 
on the assumption of crystalline fields of varying 
types, would enable us to determine the proper 
symmetry characteristics of the crystalline field 
present in iron alum. On the basis of second- 
order calculations alone, it would appear from 
Table I that a trigonal field of axial symmetry 
leads to slightly better agreement with experi- 
ment than a crystalline field of predominantly 
cubic symmetry. However, when third-order 
terms in the cubic field are taken explicitly into 
account in the calculations, the agreement be- 
tween theory and experiment is found to be much 
improved. 

Further examination of Table I might lead one 
to suspect that had third-order terms also been 
calculated for the crystalline field of axial sym- 
metry, the theoretical values of M might be 
drawn still closer to the experimental values. To 
take this statement as evidence in support of a 
trigonal field of axial symmetry in preference to 
one of cubic symmetry would, however, be 
definitely misleading since the difference between 
the two cases is probably less than the experi- 
mental error involved in making the measure- 
ments. Furthermore, it is now quite possible that 
a more general trigonal field consisting of a 
combination of the axial and cubic terms would 
also fit the facts, but as long as the agreement is 
so close as indicated in Table I there is no means 
of deciding on the exact form this arbitrary field 
should take. The one definite conclusion that can 
be drawn from the evidence presented is that 
values of M calculated to the third order in the 
cubic field perturbation are in somewhat closer 
agreement with experiment than values based 


TABLE II. Temperatures and corresponding field strengths, 











T 
[From 
Ha H T (calc.) Ea. (1)] % difference 

18,100 18,000 1.18 1.18 —_ 
15,330 15,200 1.00 1.00 — 
12,330 12,200 .800 .803 0.37 

9,270 9,140 .600 .604 0.67 

6,210 6,080 .400 .405 1.25 

3,170 3,040 .200 .206 3.00 

2,400 2,275 .150 .156 4.00 

1,630 1,510 .100 .106 6.00 

1,305 1,190 .080 .085 6.25 

1,140 1,030 .070 .074 6.25 
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only on second-order calculations with an axial 
field assumption. 


Determination of True Temperature 


To find the true temperature appropriate to 
any value of the field strength after the beginning 
of adiabatic demagnetization, the procedure is 
to solve the entropy equation [(22) plus (40) of 
Part II] for T in terms of H. Since the entropy 
retains a constant value during the process of 
demagnetization, the resulting equation im- 
mediately enables us to find T for any H or H for 
any 7, provided only that the temperatures are 
not so low that the expansion process employed 
in the calculation of the adiabatic moment 
breaks down. In Table II are shown the values of 
the temperatures appropriate to a number of 
different fields." 

In column 4 of Table II appear the tempera- 
ture values one would obtain by using the simple 
de Haas-Wiersma procedure discussed in the 
Introduction of this paper. In principle their pro- 
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cedure is never applicable to iron alum, but it is 
clear from inspection of the table that the temper- 
ature estimated by use of the de Haas-Wiersma 
equation [Eq. (1) ] differs very little from the 
exact temperature for all values of the field down 
to about 7000 gauss. Between 7000 and 1000 
gauss, the percentage error incurred by use of 
Eq. (1) rises from about 1 percent to about 6 
percent. For fields less than 1000 g the expansion 
process employed in the determination of the 
entropy, and hence of the exact relation between 
T and H, is no longer valid. 

On the basis purely of the experimental rela- 
tion between the adiabatic moment and the field 
strength, Casimir and de Haas* have noted that 
for fields down to 2500 gauss the error in tem- 
perature measurement caused by use of Eq. (1) is 
probably less than 0.01°K. Our calculations show 
an error in temperature measurement at 2500 g 
of 0.006°K. This appears to be the maximum 
absolute error involved at all values of the field 
for which Table II and our method of calculation 
are valid. 


PART II: MATHEMATICAL ANALYSIS 
V. Preliminary Details 


We require expressions for the magnetic moment and the entropy of an assembly of atomic magnets 
subject both to an applied magnetic field and a crystalline electrical field and coupled together 
magnetically. These are most easily obtained from the partition function by means of the equations 


C) 
M=kT—(log 2), 
0H - 


The partition function is given by 


a 
S=—(kT log 2). 8 
—, og Z) (8) 


Z=>) exp (— W,/kT), (9) 


where the summation is over all states of the entire crystal considering the crystal as a giant molecule. 
Since we are dealing with the partition function instead of the individual energy levels of each atom, 
we are concerned only with the energy W) of the crystal as a whole. 

The Hamiltonian function has already been described in Section II. We choose a system of repre- 
sentation which diagonalizes the first (or Zeeman) term of and then consider the spin-spin coupling 
and the crystalline field potential as small perturbations. The unperturbed energies are given by 


W,°= — gBH Li Mi, 


(10) 


where the M; are the characteristic values of S,, and the letter \ is used to represent the totality of 
quantum numbers for a given state of the entire crystal. The crystalline energies W, are given in 


“Values of H in Table II are not whole numbers, as it was found simpler to assume values of T and then determine 


corresponding values of H. 
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terms of the unperturbed energies by the expression 
Wy = W,°+ Wit W,, (11) 


where W, and W; are the first- and second-order energy corrections corresponding to the perturbing 
potentials. If now we replace these latter two terms by the well-known quantum mechanical ex- 
pressions from perturbation theory, substitute Eq. (10) into Eq. (11), and Eq. (11) into Eq. (9), 
expand the perturbing terms of the partition function in powers of kT, and retain terms to the second 
order in the perturbing potentials, we find the partition function to be 


Z=Z((1+A+B+--:), (12) 



















with 
A=—(1/RkTZo) © H™ (A; A) exp (— W°/RT), (12a) 
B=(1/2k?T*Zo) Dyn H(A; AIC (A’; A) exp (— Wy°/kT) 
+(1/kTZo) Drv: {I (A; 0’)5eM(N'; 4)/(Wyr — Wa)} exp (— Wy°/RT), °(12b) 
Zo= La exp (— Wi°/kT) = (Du; exp { — Wi(M;)/RT} }¥. (12c) 
In these expressions Zo represents the partition function in the absence of the perturbing potentials 
and #™(A; A), which represents the perturbing part of the Hamiltonian, includes both the dipole- 
dipole coupling >> js; wi ;(A; A) and the crystalline field potential >>; Vi(A; A). The equality and prime 
signs on the summations are in accordance with the notation of Van Vleck.® The equality sign indi- 
cates that only those values of \’ for which W,, = W,, or in the later Eq. (13b) those M;,’ Mj for which 
Wi(M,)+W(M;) = Wi(M)+W MZ), are to be included in the summation; the prime sign means 
that just these values are to be excluded. 
Now let 
Liri(wi)w=L ii dd; DL Mim; wi;(M:M;; M;M;) exp [— (Wi( Mi) +Wj(Mj))/kT ], 
Di (Vidw=Di di LM; Vi(M;; M,) exp [— Wi(M,)/kRT], 


d;=(Lim,; exp [— Wi Mi) /kT])“*. 






and 












where 










Then the A part of the partition function becomes 
(1) 


= —(1/RT)Hw = —(1/RkT) {Di>i (Wisdw tL (Vidad, (13a) 





while the B part is 
B= fk? 7 {Dd ier (Wis (Wedwt Diie (WijsWje)wt LD ik (Vida Vidar} si, me, i, 54%, na, — (13d) 





where 
(wi; jW jn) = dd dy Lo MiMjMy {2kT Du,’ w;i;(M;iM;; M;M;)w(Mj Mi; M;M,)/[Wi(M;,) on W,(M);) ] 


+w;;(M:M;; M;<M;)wj.(M;M:i; M;M;) } exp [- (Wi(M,) + W;(M;) + Wi(M;))/kT] 






and 
@=$k°T-*{ > >i didi Dau ayemy |Wis( MiM;; M/M;)\* exp [—(Wi(M,) +W(M)))/kT] 
+20 di Diniay | VMs; Mi) |?exp[— WM) /kT)) 
+kT'{ {Dipi did; Duujayeuy {| Wis MiM;; Mi Mj) |?/(WiA( Mi) +Wi(M;) 
— W(M;)—W;(M;))} exp [— (Wi Mi) +W(M)))/kT] 
+206 di Dovey (| VM; Mi) |*/(WA ME) — W(M))} exp [(— Wi My)/kT}}. 








Equations (13a) and (13b) are more general than the corresponding equations given by Van Vleck’ 
as provision is here made to take account of a cystalline field contribution to the partition function 
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as well as of a contribution from the magnetic dipole-dipole coupling. It may be noticed that only 
terms linear or quadratic in each interaction appear in these last equations whereas one would expect 
the operations of Eq. (12) to lead also to cross terms, i.e., terms containing a product of the two 
perturbing potentials. Such terms do arise but have been omitted from Eqs. (13) as it can be shown'® 
that for our application their contribution to the magnetic moment vanishes. This being so, the com- 
plete contributions of the two perturbing potentials to the partition function are directly additive 
and may be considered separately. 








VI. Partition Function to Second Order in Dipole-Dipole Coupling 






In order to determine the contribution of the dipole-dipole parts of Eqs. (13a) and (13b) to the 
partition function, the matrix elements of this coupling must first be evaluated. This evaluation has 
already been carried out in a paper by Van Vleck,’* henceforth called I.c., in which the partition func- 
tion is calculated to the second order in the magnetic coupling for the case of arbitrarily large applied 
fields and for arbitrary direction of field relative to cubic axes. We have'’ 








Z/Zo=1—4NR“TBy?Q0+-3 NRT —7007[ Bi (NN — 4) +4B 7B | 
+4Nk*T—*0,[3Bi'— 3B B+ { Bo? + fe(S?+S— Be)?— WB? +3 Bi(S?+S— Bz) (kT/g8H) J 
+1Nk?T~*0.[3Bi(S°+S—B2) — Bi +B,Bz \(kT/g8H) 

+iNk?T—0;[ 95 (S?+S— B2)?— 7 BY — }(S°+S—B2)Bi(RT/gBH)], (14) 








where the 2’s are functions of the coefficients of the product terms S¢iS¢ (q=xyz) that appear in 
the expression for the dipole-dipole coupling referred to a set of axes fixed relative to the crystal 
instead of to the magnetic field direction. It has been shown in I.c. that the 2’s can be expressed as'® 










Q = — Ng*B*h, %=2N*gZtB*Qo, Q2=FN*g*B4Qo, Qs=— N*g*B*Qo, (14a) 


where Qo=2N-? >> ri; This summation, varying as the inverse sixth power of the atomic distance, 
converges rapidly. The value of Qo for a face-centered lattice such as that occupied by the para- 
magnetic ions in iron alum is 14.4. The quantity ® is defined as the sum of the Lorentz local field 
factor (44/3) and the demagnetization factor (n). It takes account of the fact that the H of Eq. (14) 
which represents the applied field, differs from the true field acting on the paramagnetic ions. 

The B functions of Eq. (14) are Brillouin functions defined by 


B,=(1/Q)(d"Q/06"), 












(15) 






where 





Ss 
Q= > exp(M@) and 6=(g8H/kT). 


M=—S 
VII. Matrix Elements of Crystalline Field 


For simplicity we consider first the matrix elements for a trigonal field containing only second- 
order terms, i.e., a crystalline field of axial symmetry. Such a field is represented by Eq. (4a) of 













% J. A. Sauer, Doctoral Dissertation, University of Cambridge (1941). The cross terms are found to drop out regardless 
of the shape of the specimen. Essentially the reason is that these terms are linear in the representation coefficient and 
hence vanish when the spatial averaging is performed. 

16 J. H. Van Vleck, Phys. Rev. 52, 1178 (1937), henceforth I.c. 

Our Eq. (14) differs from Eq. (21) of l.c. only in the absence of a term in 24. This term was needed in the study of 
anisotropy but for our ene to demagnetization experiments performed on a powdered specimen this term can be 
omitted as it vanishes when averages over all directions are included. 

8 In the 0’s as defined by Eq. (37) of I.c., there also appear terms in the exchange integral. These have been omitted 
from - Eq. (14a) as exchange effects can be neglected in our treatment because of the known magnetic diluteness of 
iron alum. 
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TABLE III. Matrix elements of crystalline field of axial symmetry. 











My 5/2 3/2 1/2 —1/2 —3/2 —S/2 
5/2 | D®/V/6 —D3"/V/5 D*//10 0 0 0 
3/2 | D¥/VS5 —3D"//6 —4D2"/V/2 #D”/\/2 0 0 
1/2 | Ds*//10 gD2¥/V/2 —D//6 0 1D”//2 0 

—1/2 0 #D.-™/./2 0 —tD”/./6 #D3"/4/2 D*//10 
—3/2 0 0 $D.-™/4/2 —4#D"//2 —4D"//6 D:*//5 
—5/2 0 0 0 D:-*/4/10 —Ds"/V/5 D®//6 











2 
Note: All matrix elements to be multiplied by 7 +/10 or, in the notation of Penney and Schlapp, by 10Aa //6. 


Section II or by the equivalent expressions 


24/10 


V.=Ar?(1—3 cos? 6) = —Ar? T:, (16) 





where 7,° is a spherical harmonic of the second order normalized to unity. Since the experiments of 
Casimir and de Haas were performed on a powdered specimen in which all directions of magnetic 
field relative to the cubic axes are equally probable, it is necessary to transform from the x, y, g 
system relative to fixed axes over to a new coordinate X, Y, Z system relative to the arbitrary mag- 
netic field direction. Under this transformation the Tesseral harmonic of order two will transform 
into a sum of other Tesseral harmonics of the same order. In particular, 72°= }>m D2”’:°T.”’ where 
the D,"”’’s are the so-called representation coefficients of order /. The desired matrix elements are 
now found by use of the Kramers symbolic method*'* or by direct computation from the funda- 
mental quadrature. The matrix representing V, in the ‘‘M”’ system of representation is given in 
Table III. The symbol J, multiplying the matrix elements is given by J2=cy fr°R*dr, where cy is an 
undetermined parameter that enters in the Kramers group-theory method of calculation of the 
matrix elements and R is the radial wave function. The value of Jz may be determined by comparing 
the diagonal element for M=5/2 with the corresponding matrix element determined by Penney and 
Schlapp.’ Their value is applicable only for a pure crystal in which the magnetic field is directed 
along the z axis. Taking this restriction into account arid comparing, we find that our expression 
#I,A 4/10 is the same as 10Aa¥/6 in the notation of Penney and Schlapp. 
If the second-degree terms of the trigonal crystalline field are omitted and only the fourth-degree 
cubic terms considered, that part of the crystalline potential independent of r can be written 


2 2 2 
V.=—v2T,°+—(35)'T,'+—(35) 'T4, (17) 
15 105 105 


where the 7,”’s are fourth-order harmonics normalized to unity. In this expression @=0 refers toa 
direction along one of the cubic axes. In transforming over to a new coordinate system with polar 
axis in direction of magnetic field, the 7,"’s transform like representation coefficients of the rotation 
group of dimensionality 2/+-1 so that 


T," == > Dy'*T,". 


The matrix table for any one of the fourth-order Tesseral harmonics is now obtained by use of the 
Kramers symbolic method. The matrix elements for that combination of spherical harmonics given 
by the cubic crystalline field of Eq. (17) are given in Table IV. In this table, the symbol D,” stands 


1H. C, Brinkman, Doctoral Dissertation, University of Utrecht (1932). 
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TABLE IV. Matrix elements of crystalline field of cubic symmetry. 








My 5/2 3/2 1/2 —1/2 —3/2 
5/2 D§//70 —2D¢//70 3D2/+/70 —Dé/V/5 Dé/V5 
3/2 2D,°//70 —3D/+/70 D2/V7 —D2/V/14 0 
1/2 3Dy?//70 —DeoV/V7 2D8//70 0 —D?2/V14 

—1/2 Dy/V5 —Dy?/V 14 0 2D8//70 —D¢/V7 
—3/2 Dyo*/V/5 0 —Dy?*/V/14 Dyo/V7 —3D°//70 
—5/2 0 Do4*/V5 —Dy/V5 3Dy?/V70 —2D,"/V/70 














2 
Note: All matrix elements to be multiplied by og DE or, in the notation of Penney and Schlapp, by 60Dg. 


for the following combination of representation coefficients : 


Dy =(14/5)'!Dg’+Dye'+De— = and Tames f r'Rer. 


The switch to the Penney and Schlapp notation is carried out in the same manner as indicated above. 

Numerical Checks: One check on the numerical accuracy of both matrix tables that we have 
verified is that those elements which arise due only to even values of m in T," or T4" must agree with 
values previously calculated by Van Vleck and Hebb” in connection with a different problem. 

Secondly, since the diagonal terms of the cubic field matrix were used to change from our notation 
to that of Penney and Schlapp, a comparison of the off-diagonal terms in the two cases can be used 
as a check. All that is necessary is to give the representation coefficients their specific numerical 
value for the special case considered by Penney and Schlapp, and when this is done and the com- 
parison made, the two expressions are found to be alike. 


VIII. Partition Function to Second Order in Crystalline Field Interaction 


_ The partition function for the crystalline field terms is now obtained by substituting the matrix 
elements from Table III or IV into Eqs. (13a) and (13b). The resultant expressions are functions of 
the D,”™’’s that appear in the matrix elements. These representation coefficients, however, can be 
eliminated by performing the spatial averaging required in virtue of the powdered form of the test 
samples.” The averaging is easily carried out by using the group theory’relationship 


f f f DiDy "dQ / f f f d2= BnpSngd11/v. (18) 


Here 6 is the Dirac delta function and » is the dimensionality of the representation given by 2/+1. 
After carrying out the indicated procedure, we find the contribution of the axial crystalline field 
to log Z to be 


(N/45)y.2T | —2Q-°(5 cosh 5x —cosh 3x — 4 cosh x)?+@Q-'(25 cosh 5x+cosh 3x+16 cosh x) 
+3Q07(kT/g8H)(25 sinh 5x—3 sinh 3x—4 sinh x)}. (19) 


For the crystalline field of cubic symmetry, the contribution to log Z is 


(Ny2/21T*) | —2Q-*(cosh 5x —3 cosh 3x+2 cosh x)?+@Q-'(cosh 5x+9 cosh 3x+4 cosh x) 
+40-'(kT /3g8H)(25 sinh 5x+18 sinh 3x—11 sinh x)}. (20) 


2% H. Van Vleck and M. H. Hebb, Phys. Rev. 46, 17 (1934). 
*1 The averaging should be carried out not on Z itself but on log Z. The reason why Z can’t be used is clear when we 
recall that log (1+) where x is small compared to unity is equal to x—4}x* to the second order’. Hence, if average was 


2 
taken - Z itself instead of log Z, we would be dealing with (x)w instead of with (x*),, and these are by no means 
equivalent. 
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In these expressions x = (g8H/2kT) and +, and 7, have taken the place, respectively, of the crystalline 
field constants 6Aa and 24D,. The y symbols may be considered as characteristic temperatures 
associated with the crystalline field splitting pattern. If we set up and solve the secular equation for 
the one atom problem, we find that the y’s of Eqs. (19) and (20) are exactly } the é’s of Fig. 1. 


IX. Magnetic Moment and Entropy to Second Order in Perturbing Potentials 


Expressions for the magnetic moment and entropy correct to the second order in the perturbing 
potentials are now found by substituting Eqs. (12c), (14), and either (19) or (20) into Eq. (12) and 
using Eq. (8). We obtain: 


(M/Ngp) =Bit(1/T) fart (t/T)*( fare +Gu) + (/T)*2, (21) 
S/Nk=log Q— (g8H/kT)B1+ (1/T) fe +(1/T)*( for +Gs) + (¥/T)*u. (22) 


Here r is a characteristic temperature associated with the dipole-dipole coupling and for iron alum 
with S=5/2 it takes the value 


r= Ng?s?S(S+1)/k=0.0472°. . 
The f’ and f’”’ terms of Eqs. (21) and (22) are the first- and second-order effects of the local field plus 
demagnetization corrections while the G terms, entering only in the second order, represent the 
perturbations due to the dipole-dipole interaction that are not included in the ® terms. They are 
defined by lit 
fur =®S—'(S+1)-'B,(B.—B,’), (21a) 
fur =8°S-*(S+1)-*{B,(B2.—2B;*)(B2—B,*) +3B1*(B;— BiB2)}, (21b) 
Gu =%}Q0S-*(S+1)-*{ (Bs—BiB2)[—Bv2+B2—3(S?+S—B2) +3} (kT/g8H)B; ) Si 
+ (B2—B,*)[2B,*—2B,B.—}B,+4}(kT/g6H)(S?+S+B2—6B,) ] 
+(kT/g8H)* — $Bi(S?+S—B2)—$BiB2+ 3Bi*)}. (2c) mi 
f= —®S-\(S-+1)-1(g8H/kT)B,(Bs— By’), (22a) 7 
fur = ¥®2S-2(S+1)-2{ By — By2Bs — (g8H/kT)B1?(Bs— B,B2) fie 
— (g8H/kT)2B,(B2—2B,*)(B2—B,*)}, (22b) ” 
G.= — (Qo/10).S-?2(S+1)-?{ Bit —2B;°B.+ B.?+3(S°+S—B.)?— B+ $3Bi(B;s— BB) vn 


+$(B.—B,*)[(S?+S—B:) —6B,°+2B2]+(g8H/kT)[(B2—B,*)(4B;*—4B,B.—}B;) fiel 
+ (Bs— B,B:)(2B2—2B,*—}(S*+S—Bs2))}}. (22c) 


These equations are true for any value of the spin S, but a convenient check on the correctness of 7 
the algebraic manipulations involved in the derivation of the above equations can be had by re- pra 
stricting the atomic spin to the value 3 and then comparing these expressions with the corresponding Eq 
ones given by Van Vleck.’ After making the required restriction, we find complete agreement. tere 
The v and u terms of Eqs. (21) and (22) represent the perturbations due to the crystalline field. tot 
Consider first the field of axial symmetry and let es 
C*(a, b, c) =Q-"(a cosh 5x+5 cosh 3x+<¢ cosh x)", ni 


S*(a, b, c)=Q-"(a sinh 5x+5 sinh 3x+c sinh x)" — 
with 





Q=2(cosh 5x+cosh 3x+ cosh x). 
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Then 
45y.=4B:C2(5, —1, —4) —2C(5, —1, —4)$(25, —3, —4) —B,C(25, 1, 16) +38(125, 3, 16) 
+4(kT/g8H)[C(125, —9, —4)—2B,S(25, —3, —4) —2(kT/g8H)S(25, —3, —4)], 
45u,=2C2(5, —1, —4)—C(25, 1, 16) —3C(125, —9, —4)+3B,S(25, —3, —4) 
+(g8H/kT)[—4B,C%(5, —1, —4)+B,C(25, 1, 16) —}8(125, 3, 16) 

+2C(5, —1, —4)S(25, —3, —4)].° (22d) 






(21d) 











For a trigonal field of cubic symmetry, the v and u symbols are given by” 
21v-=4B,C2(1, —3, 2)—2C(1, —3, 2)S(5, —9, 2) —BiC(1, 9, 4)+4S8(5, 27, 4) 
+32(kT/g8H)[C(125, 54, —11) —2B,S(25, 18, —11) —2(k7T/g8H)S(25, 18, —11)], (21e) 








4 
21u.=2C%(1, —3, 2)—C(1, 9, ote. 18, — 11) —3C(125, 54, —11) 






— (g8H/kT)[4B,C2(1, —3, 2)+48(5, 27, 4) —B,C(1, 9, 4) —2C(1, —3, 2)S(5, —9, 2) ].  (22e) 








If one neglects the correction terms in Eqs. (21) and (22) arising from the dipole-dipole and crystal- 
line field interactions, the expressions for M and S reduce to 


M=N¢g6B,, 
S=Nk(log Q—6B,). 











(23) 
(24) 










Since S remains constant during an adiabatic demagnetization process, it can readily be seen from 
the above equations, remembering that 6=(g8H/kT), that the temperature is a linear function of 
the field of the form T=C,H+C, with C2 zero. The procedure of de Haas and Wiersma for deter- 
mining the true thermodynamic temperature now clearly applies. It has been shown by Van Vleck® 
that even if the correction terms in f’ and f”’ are retained Eqs. (23) and (24) are still valid provided 
only that the magnetic field H be no longer interpreted as the applied field but as the true local 
field H+4M. The linear equation between T and H still holds, but the coefficients now satisfy the 
relation C2/C,;=®M. 

Thus, even for specimens of non-spherical shape—for spherical specimens the ® terms are always 
zero—the f terms may now be omitted from Eqs. (21) and (22) though care must be taken when com- 
paring the results of theory with experiment first to convert the applied field over into the true local 
field by introducing the proper demagnetization corrections for the shape in question. 














X. Determination of Adiabatic Magnetic Moment 





To derive an expression for the magnetic moment valid throughout the entire demagnetization 
process, one needs only to solve Eq. (22) for T in terms of H and then substitute this function in 
Eq. (21). In this manner one obtains an equation involving only M, Mo, and H. We are primarily in- 
terested in the form of M for large initial fields and therefore assume, as a comparable approximation 
to those previously made, that the primary terms in M and S in Eqs. (21) and (22) can be expanded 
as a Taylor’s series in @—@ and only first-order terms retained. Furthermore, it is assumed that the 
argument of the v, u, and G terms, which are all small perturbations, has the same value as at the 
initial helium temperatures. This means that @ can be everywhere replaced in these terms by @. 












* One helpful check on the algebraic accuracy of Eqs. (21) and (22) which has been useful is: 
(0S/0H)r=(8M/8T)x. 
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After expanding the leading terms in M and S, we find 
(S—So)/Nk=0= —6[B2° — (By*)"](0— 6°) + (7/T)*G.°+ (y/T)?u°, (25) 
(M— Mo)/Ng8=[B2°— (B,?)°](0— 6°) + (r/T)*Gu°+(y/T)*, (26) 


where the superscript zero in the B,, G, u, and v terms means that these are to be evaluated at 
6= 6 = gBH)/kT». From Eqs. (25) and (26) we get 


tT \2/G,° vy \2su° Hy\? 
M—Mo= Ne8| (—) —+Gu')+(~-) (+) }(=) : (27) 
To 6° T) 9° \H 
In order that this equation apply for both axial and cubic crystalline fields, u° and v® should be taken, 
respectively, as u.°, v.° and u,°, v.’. 
Equation (27) is the theoretical expression for the adiabatic magnetic moment correct to the 

second order in the perturbing potentials. For purposes of comparison with the experimental relation 
between M and H found by de Haas and Casimir, it can be put in the form 


M=M,(1—6/H?), (28) 


H.? rT \2/G,° v¥\2/u® 
ad a EG *)} oo 
B,° To 6° To 6° 

Another convenient check on our calculations is now available. De Haas and Casimir? have given 
a theoretical derivation for 6 that is accurate only for the case of low initial fields. Our analysis is 
valid for fields of any nature and hence our expression should reduce to theirs if we restrict H to 
low values only. The procedure is to expand the various correction terms G,, ua, etc., that appear in 
Eq. (29) in powers of @ and then retain only the lowest terms in the expansion. When this is done we 
find b=3A'/2r and this is equivalent to the expression given by de Haas and Casimir.** The nu- 
merical value of 3A’/2r is 1.3 10°. This should not be expected to agree with the experimental value 
of b of 1.7105 since the demagnetization experiments were carried out at such high values of the 
initial field that saturation effects cannot be neglected. 

For the more important case of high fields, the numerical value of b must be taken directly from 
Eq. (29). The characteristic magnetic coupling temperature + has already been evaluated and 
numerical values of G,° and Gy°,* representing the effect of the magnetic coupling on the adiabatic 
moment, can be calculated directly from Eqs. (21c) and (22c) as these equations depend only on 
#, Qo, and S, whose respective numerical values are 2.04, 14.4, and 5/2. The v and u terms of Eq. (29) 
can also be directly evaluated™® from their defining equations, but the complete crystalline field 
contribution to the adiabatic moment cannot yet be found as this contribution contains an unknown 
crystalline field parameter, vis., y. or y-. Hence, it is first necessary to evaluate these parameters 
and this is done by calculating the theoretical expression for the specific heat in zero applied field 
and comparing it with the known experimental value. 


where?’ 





Crystalline Field Constants for Specific Heat Data 


The theoretical expression for the specific heat is obtained directly from the partition function by 
means of the equation 


0 i) 
oT oT 


* The bracketed part of Eq. (29) is always found to be negative; hence b itself is essentially positive. 

%* The expression actually given by de Haas and Casimir is A /2a where small a is the coefficient of H/T in the ordinary 
formula for M when saturation is neglected. However, since in our notation this coefficient is r/3 and since our A’ is 
identical with their A, the two expressions are equivalent. 

% The calculated values of these various quantities are stated explicitly in reference 15. 


an 


nat 
ex{ 
fur 


cal 
fro} 
tab 
me! 
fro 
cry: 
wit 


The 
















ADIABATIC DEMAGNETIZATION OF IRON ALUM 109 








We desire the specific heat in the absence of an applied field. Taking into account this condition and 
substituting Eqs. (14) and (19) into Eq. (30) we find 












































25) 

26) C.= Nk{2.4(r/T)?+1.55(yo/T)*}, (31) 
where the first term represents the contribution of the dipole-dipole coupling and the second term 
at . A ie ° ° . . 

that of an axial crystalline field. If the partition function derived for the case of a cubic crystalline 
field be substituted in Eq. (30), then 1.55y,? should be replaced by 2.0072. 
| Equation (31) is of the form A’T~ as found experimentally by various observers.” Since the value 
27) of r is known, this equation may be used for determining y, or y- by substituting for C, on the left- 
hand side of the equation the accepted experimental value of 23007~-*. The crystalline field parameters 
en, now become y-=0.063 and y.=0.072 corresponding to splitting distances (see Fig. 1) of 6,=0.190 
and 6,=0.215. Substituting these values and those of the other parameters into Eq. (29), we find 
the b=9.4X10' for a field of cubic symmetry and =11.5 X10 for a field of axial symmetry. 
‘io 
" Effect of Demagnetization Corrections on Adiabatic Moment 
28) For either type of crystalline field, the calculated values of the coefficient b are found to be too 
| low and hence the theoretical values of M are higher than the experimental. As mentioned in Part I, 
, some of the discrepancy between the theoretical and the experimental values is due to the fact that 
a the H of Eq. (28) is no longer the same as the applied field H. of Eq. (2). To find the true local field 
H when H, is known, it is first necessary to determine the proper demagnetization factor » for the 
ven shape of specimen used in the experiments. De Haas and Casimir used a powdered specimen of iron 
s is alum pressed in the form of a cylinder 90 mm long by 30 mm wide. Since the true demagnetization 
T to factor is calculable only for an ellipsoid, all one can do is to estimate this factor by assuming that 
rin the specimen can be adequately represented by an ellipsoid of revolution with an axes ratio of 3 : 1 
we The magnetic field was applied at right angles to the longitudinal axis of the ellipsoid. For this case 
nu- @=(4/3)r—n becomes negative and hence the local field H(=H,+4M) is less than the applied 
lue field H,. This has a tendency to lower the theoretical values of M, thereby bringing them in some- 
the what closer relation to the experimental data. However, a direct evaluation shows that while a minor 
improvement has been made the major part of the discrepancy between the theoretical values of M 
‘om and the experimental still exists. 
and 
atic Effect of a More General Crystalline Field 
ee In Section III of Part I, the thought was raised that possibly a crystalline field of a more general 
= nature than one of cubic or axial symmetry would lead to better agreement between theory and 
eld experiment. To examine this question more thoroughly, we calculate the contribution to the partition 
= function and the magnetic moment of the arbitrary crystalline potential specified by Eq. (6) of 
a, Section III. We need the matrix elements of this potential and can easily obtain them from the 
eld calculations already made. The matrix elements of the 7; terms of Eq. (6) can be obtained directly 
from the numerical values listed in Table II, and each of the 7, terms of Eq. (6) leads to a matrix 
table similar to that of Table III. The only change that must be made in these tables is the replace- 
ment of the original coefficients of the 7;’s that appear in Eqs. (16) and (17) by the new coefficients 
by from Eq. (6). On substitution of the matrix elements into Eqs. (13) the contribution of the arbitrary 
crystalline potential to the log of the partition function is found to be the sum of 7“ (19) and (20) 
30) with y,2/45 replaced 7 
4 
(S, (d2r)?/750) and y.*/21 replaced by: (20. (au)?/630). 
—2 4 
Vs The expressions for the magnetic moment and entropy previously given by Eqs. (21) and (22) can 


be made to apply to the general crystalline field, by replacing, respectively, the y*v and yu terms of 
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these equations by the following expressions: 
‘ 2 : 2 ‘ 2 : 2 
(v-/30) Dos das+ (302/50) Se a2, and (u-/30) Ss d4e+ (30/50) ¥, aay. 
a" —2 all —2 
The adiabatic moment and the specific heat found in the same manner as previously indicated, are 


M=M,(1—},/II®) 


by = — (H?/B »{ (2) (Ste en \+e Eicto(e ‘)| ; 
ONG) VOD 30T Ne 50 Te Mot “ 


1 2 7 2 | 
C,= NkT-*) 2.47°4+--— 3, 08,4+-— 5, 03, 
| MATS De tutes Le dar _ 6) 


with 








and 


These are the exact expressions correct to the second order in the perturbing influences which 
include both the magnetic dipole-dipole coupling and the arbitrary crystalline potential of Eq. (6). 
All the coefficients of this potential are at our disposal and the problem now is to give b, and C, 
their experimental values for iron alum and then determine whether a suitable choice of coefficients 
can be found that will lead to agreement in both equations. The results of this procedure have already 
been discussed in Section III of Part I. 


XI. Adiabatic Moment to Third Order in Crystalline Field Potential 


The problem of the determination of the energy levels of iron alum supposed subject to a large 
applied magnetic field of arbitrary direction and to a crystalline electric field of cubic symmetry 
has been considered by Kronig and Bouwkamp." They were able to calculate, correct to the third 
order in the cubic field splitting parameter, all six energy levels of the normally degenerate ground 


state. The six expressions are :*° 
2qy*? 4ny* 

+ tess, 
gBH 2°8*H? 





‘ 3 
W g 


39B8H —3 dary” er’ . 
=+328H—3py+ + t:::, (34 
| . ? gBH 98°? 





2 
W; 
W, 
2gsy" ss 4rsy’ 

+ a 
gBH °° H? 





= > eH +4p 
=+-2 yt 
W, 2 
where 


5 Fy 25y) aoa (22—75y) Sine (50—113y) 
=1-— ’ = = ’ 2=—+—¢ wx Y), —— IP), 
P wea ” Os 2 1932 96 


- 5 15 
r1=—-(196 — 16359+3125¢2), r2=——+——o(79—6159+1125y), (35) 
144 128 128 
-5 25 
r3=—————- (113 —705 9 +1075). 
128 1152 


In these equations ¢ is a specific function of the direction cosines a, a2, a3 of H with respect to the 


26 Our symbol y is the same quantity as that indicated by a in the Kronig-Bouwkamp paper. 
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cubic axes given by 
9 = 27013? + 3a? +172”. 






Let the six energy levels of Eq. (34) be substituted into the partition function 


Z=[>du exp (— Wau/kT))%. 





Expanding the exponential and retaining terms only to the third order in the cubic field potential, 





we find 





1 
£/Zo= | 1—(7/T)pC(1, —3, 2)+(y/ r)| pC, 9, 4) +4(kT/g8H)S(qs, g2, «| 





N 





1 
- @/T)| peu, — 27, 8) +8(RkT/g8H)°C(rs, r2, 11) +2(kT/g8H)S(pgs, —3p¢2, 2001) (37) 





In order that our final results be applicable to a powdered specimen, it is neeessary to calculate the 
average value of the various quantities 1, r2, g3, etc., that depend on the direction cosines. Performing 
this averaging,*® expanding log Z, and retaining terms only to the second power in y, we obtain an 
expression for the log of Z that is identical with our previous Eq. (20) and thereby provides another 
check on the algebraic manipulation involved in its derivation. Since the effect on the magnetic 
moment and entropy of that part of the partition function given by Eq. (20) has already been con- 
sidered, it is only necessary for us to consider now the third-order terms of Eq. (37). The contribution 
of these terms to the average value of (1/N) log (Z2/Zo) is found to be 











(1/3003) (y/T)*{| —6C(1, —27, 8) —48C*(1, —3, 2)+36C(1, —3, 2)C(1, 9, 4) 





+ (kT /gBH)[S(125, 453, —664)+2C(1, —3, 2)S(—125, 151, — 332) ] 










+3(kT/g8H)?C(12295, 9441, —2854)}. (38) 








The additional term in the magnetic moment, due to inclusion of the above third-order term in 
the partition function, is found by differentiating Eq. (38) with respect to H and multiplying by kT. 
This term is 





Ng8(y/T)*wa, (39) 
way = (1/7) way’ + (1/3003) wy,” 





where 







with 
wy’ = —(kT/gBH)\ (RT /g8H)S(.268, 1.05, 1.55) +2(kT/g8H)*C(4.78, — 3.67, 1.04) 













+ B,S(.268, 1.05, 1.55) — $C(1.34, 3.15, 1.55) 






+B,C(4.78, —3.67, 1.04) — $8(23.9, —10.0, 1.04) 





and 
wy’ =+6B,C(1, —27, 8)—3S8(5, — 81, 8)—72C*(1, —3, 2)S(5, —9, 2)+144B,C%(1, —3, 2) 
—72B,C(1, —3, 2)C(1, 9, 4)+18S(5, —9, 2)C(1, 9, 4)4+18C(1, —3, 2)S(5, 27, 4) 
. +(kT/g8H){ —2(kT/g8H)C(1, —3, 2)S(—125, 151, —332) 

—4B,C(1, —3, 2)S(—125, 151, —332)+S(5, —9, 2)$(—125, 151, —332) 














+C(1, —3, 2)C(—625, 453, —332)}. 
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The contribution to the entropy of the third-order terms in the partition function is found to be 
Nk(7/T)*w,, (40) 
w, = (1/7)w,’ + (1/3003) w,”” 










where 





with 
w,' = +(kT/g8H)[ —S(.268, 1.05, 1.55) +C(4.78, —3.67, 1.04) 
+48(23.9, —10.0, 1.04) ]+B,S(.268, 1.05, 1.55) —$C(1.34, 3.15, 1.55) 






and 
w,"’ = (g8H/kT)[ —6B,C(1, —27, 8)+3S(5, —81, 8) —144B,C%(1, —3, 2) 
+72C7(1, —3, 2)S(5, —9, 2)+72B,C(1, —3, 2)C(1, 9, 4)—18S(5, —9, 2)C(1, 9, 4) 
—18C(1, —3, 2)S(5, 27, 4) ]4+96C#(1, —3, 2)+12C(1, —27, 8) —72C(1, —3, 2)C(1, 9, 4) 
+4B,C(1, —3, 2)S(—125, 151, —332) —S(5, —9, 2)S(—125, 151, —332) 
—C(1, —3, 2)C(—625, 453, —332) —2(kT/g8H)C(1, —3, 2)S(—125, 151, —332). 














Complete expressions for the magnetic moment and entropy to the third order in the crystalline 
field interaction and second order in the dipole-dipole coupling are obtained by adding Eqs. (39) 
and (40) to the previous Eqs. (21) and (22) developed in Section X. The adiabatic magnetization, 
found by solving the entropy equation for T as a function of H and substituting into the expression 
for the magnetic moment, is 


wands eA) Eee) +2) Co AAC) Cee) 


where the sub- (or super-) script indicates that the quantities are to be evaluated at 0= & = g8Ho/kT, 
= 2.04. Placing numerical values of third-order correction terms—obtained directly from their 
defining equations—into Eq. (41) along with the proper values of the second-order correction terms 
and of the characteristic temperatures 7 and ¥, we finally obtain 


M = M,(1—9.4 X 10*H-*—30 X 10°H-*). 

























(42) 








This expression is exact to the second order in the dipole-dipole coupling and to the third order in 
the cubic crystalline field interaction. The agreement between it and experiment is discussed thor- 


oughly in Section IV of Part I. 
The writer wishes to thank Professor Van Vleck for suggesting this problem and for freely offering 


his assistance on many occasions 
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The electrical conductivity of cuprous oxide is studied from room temperature to 700°C. 
It is observed that whereas the conductivity may be represented by the formula ¢=A,~*/*T 
at the higher temperatures, below 300°C the exponential law is not obeyed. The maxima 
and minima of the conductivity curves observed in the latter region are presumed to result 
from a combination of a saturation of clusters and a rapid decrease of mobility with tempera- 
ture. Hall effect measurements up to 450°C show that the number of current carriers does 
not obey the exponential law of temperature dependence. An anomaly in the mobility is 
observed around 150°C where the mobility changes from a T-* to a T~* dependence. 





INTRODUCTION 


HE semi-conductor, cuprous oxide, has been 

the object of many investigations in the 
past decade because of its peculiar electrical 
properties which are of theoretical and technical 
importance. However, the results and interpre- 
tations of these investigations have differed so 
markedly as to prevent the formulation of a 
unique theory of conduction in cuprous oxide. 
Englehard' and Angello* have studied the de- 
pendence of the conductivity and Hall coefficient 
on the impurity content and past history of 
various specimens. Their results indicate that 
the conductivity increases rapidly with rise in 
temperature and with an increase in the concen- 
tration of excess oxygen in the lattice. The Hall 
coefficient is positive, indicating that the current 
carriers are positively charged. Vogt*® reached 
the same conclusion as a result of measurements 
of the thermoelectric effect. Between 100°C and 
1000°C Jusé and Kourtschatow! and Diinwald 
and Wagner® have investigated the electrical 
conductivity. Jusé and Kourtschatow determined 
the excess oxygen content of their specimens at 
room temperature and attempted to correlate 
this with the conductivity. They concluded that 
at the high temperatures the electrical con- 
ductivity was an intrinsic property of the cuprous 
oxide and depended solely on the temperature. 


1E. Englehard, Ann. d. Physik 17, 501 (1933). 
2S. J. Angello, Phys. Rev. 62, 371 (1942). 
*W. Vogt, Ann. d. Physik 7, 183 (1930). 


‘W. Jusé and B. W. Kourtschatow, Physik. Zeits. 
i 2, 453 (1932). 

*H. Diinwald and C. Wagner, Zeits. f. physik. Chemie 
22B, 212 (1933). 


Diinwald and Wagner, on the other hand, 
investigating the effect of varying the external 
oxygen pressure on the specimen at constant 
temperature observed that the conductivity 
varied as the 1/7 power of the oxygen pressure. 
This indicates that the conductivity is a function 
of the amount of absorbed oxygen. The only 
data on the Hall effect at high temperatures are 
those of Schottky and Waibel,* who interpreted 
their data initially as indicating a change in the 
sign of the Hall coefficient. In a later paper 
Schottky concluded that this interpretation was 
erroneous. Thermoelectric measurements made 
at high temperatures by Diinwald and Wagner 
showed that the sign of the current carriers 
was positive. 

The evidence of the above works indicates 
that cuprous oxide fits into the band theory of 
solids in the following manner: Pure cuprous 
oxide is an insulator which possesses a gap of at 
least one ev between the top of the highest 
filled band and the bottom of the lowest unfilled 
band. When cuprous oxide is in equilibrium 
with oxygen gas it absorbs oxygen so that the 
composition is no longer stoichiometric. The 
absorbed oxygen produces empty discrete energy 
levels which lie within the forbidden region but 
close to the top of the filled band. Thermal 
excitation of electrons to these levels results in 
electron “‘holes”’ in the lower band. The resulting 
conductivity is then apparently due to positively 
charged particles. In addition, it is strongly 
dependent upon the concentration of the ab- 


*W. Schottky and F. Waibel, Physik. Zeits. 34, 858 
(1933); 36, 912 (1935). 
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sorbed oxygen. At low temperatures (less than 
300°C) the rate at which equilibrium is estab- 
lished between the absorbed oxygen in the solid 
phase and the gas phase is so low that it is 
impossible to control the concentration of ab- 
sorbed oxygen by varying the pressure of the 
oxygen gas. Thus, the electrical properties at 
low temperatures depend upon the past history 
of the sample, such as the manner in which it 
was cooled from elevated temperatures. At high 
temperatures, where equilibrium may be as- 
sumed, the following mechanism is postulated 
for the absorption of the oxygen.’ Oxygen 
molecules are absorbed on the surface of the 
cuprous oxide, and copper ions and electrons 
diffuse to the surface to form additional cuprous 
oxide with the oxygen molecules. This diffusion 
process leaves copper ion vacancies in the 
lattice. The process may be summarized by the 
reaction 


2Cu*++402(gas)—>2Cu,+++O0-+2Cut, 


where Cu,** represents a copper ion vacancy in 
the lattice with an electron “‘hole’”’ bound to a 
neighboring Cut ion. This copper ion vacancy 
is the source of a discrete energy level which 
lies in the forbidden region. By niinimizing the 
free energy of the above system one obtains 
the relation for the number of copper ion 
vacancies 


n= NoT—*'*p} exp (— W/kT), 


where m» is the number of copper ion vacancies, 
p is the pressure of the oxygen gas, 2W is the 
energy needed to place an oxygen atom in the 
lattice and produce two copper ion vacancies, 
No is a constant, and k is Boltzmann’s constant. 
It is assumed that m, is small compared with the 
number of particles in the solid and also in the 
surrounding gas. 

Using this expression for m, we may now 
determine the number of free electron holes to 
be expected at a given temperature. It is assumed 
that the ionization of a hole may be treated as 
if equivalent to the reaction 


Cu,*+@Cu,*+ +h. 


Minimizing the free energy of the reaction with 


7C. Wagner, Trans. Faraday Soc. 34, 851 (1938). 
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respect to a change in the number of free holes, 
h, we obtain 


h=hon'T? exp (— E/kT), 


where E is the energy necessary to free a hole, 
Using the expression for nm, we have 


h=const. p'/877"6 exp (—€/kT), 


where «= W+E. 

With the assumption that the free electron 
holes may be treated as a perfect gas, their 
contribution to the electrical conductivity may 
be calculated by simple classical theory.* The 
conductivity under these conditions is given by 
o=neyu, where n is the number of free holes per 
unit volume, e is the charge of a hole, and yu the 
mobility of the carriers obtained from the rela- 
tion p= 4elo/3(2m*kT)', where Io is the average 
mean free path and m* the effective mass of the 
holes. Using the equation above for the number 
of free holes we obtain for the conductivity 


o=const. p!/®7—"!8], exp (—€/2kT). 


This equation can be expected to be valid only 
when the temperature is sufficiently high so that 
equilibrium conditions are established at all 
times. At low temperatures when equilibrium is 
not established the number of oxygen atoms 
absorbed will be independent of the temperature. 
and one would expect the conductivity to show 
an exponential dependence on temperature with 
an activation energy E. Of course, the assump- 
tion is made that the copper ion vacancies in the 
lattice do not interact with each other. If, 
however, the density of copper ion vacancies is 
large enough for this to happen, one might 
expect a decrease in the activation energy.® In 
the preceding discussion it has been assumed 
that the mean free path is either constant or isa 
slowly varying function of the temperature and 
impurity content. This should be true at the 
higher temperatures, since the theory assumes 
that the mean free path is determined by the 
interaction of electrons with lattice vibrations. 
In any case, a measurement of the Hall 
coefficient in conjunction with conductivity 
measurements allows one to determine both the 


8 F, Seitz, Modern Theory of Solids (McGraw-Hill, New 


York, 1940), P. 190. 
°P. H. Miller, Jr., Phys. Rev. 60, 890 (1941). 
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Fic. 1. Schematic diagram of furnace, specimen holder, and 
specimen in position for measurements. 


number of carriers and the mobility independ- 
ently. The Hall coefficient R is given by the 
relation R= —3x/8nec and the product of the 
Hall constant and the conductivity gives the 
mobility 

u= —8Ro/3xc. 


If one compares the simple theory with the 
results of the investigations on cuprous oxide 
the following conclusions are reached: At high 
temperatures the conductivity can be expressed 
rather well by an equation of the form 


o=o0,exp (—U/kT), 


where U is apparently the same for all samples. 
Therefore, Jusé and Kourtschatow‘ conclude 
that this high temperature conductivity is 
intrinsic. On the other hand, the results of 
Diinwald and Wagner® indicate that the con- 
ductivity depends upon the excess oxygen con- 
centration. Furthermore, the work of Waibel 
and Schottky® on the Hall effect at high temper- 
atures indicates that the number of current 
carriers is not a simple exponential function of 
the temperature, and that the mobility is a 
very rapidly varying function of the temperature. 


EXPERIMENTAL PROCEDURE 


The specimen holder as shown in Fig. 1 is a 
lavite block which is milled to accommodate 
the specimen. The Hall and conductivity con- 
tacts, as indicated, are made of Nichrome wire 
which offers the advantage of being able to 
retain its spring at the higher temperatures. 
Enclosing the specimen holder is a rectangular 
transite furnace in which are embedded Nichrome 
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heating elements. These are covered with a 
mixture of Alundum cement and sodium silicate. 
Narrow strips of transite separate the two fur- 
nace walls and the entire unit is clamped to- 
gether with screws. 

The furnace and hokler are mounted between 
the pole faces of an electromagnet which is 
made of two cylinders of Armco iron. With 12 
amperes through the windings and an air gap of 
+’ the field is 7150 gauss at room temperature, 
as calibrated with a bismuth spiral. All electrical 
measurements are made with a Leeds and 
Northrup type K potentiometer. The tempera- 
ture of the furnace was controlled manually, 
and temperature measurements of the specimen 
were made with two Chromel-Alumel thermo- 
couples placed at opposite ends of the sample. 

The cuprous oxide was prepared from strips 
of Chilean copper, 2}’’X}”" Xe" which were 
suspended in a furnace at 1000°C for 16 hours.'® 
This time was sufficient to effect a complete 
oxidation of the specimen. After oxidation the 
sample was placed in an annealing furnace for 
two hours at a predetermined temperature and 
then quenched in cold tap water. 

Gold contacts, about 0.01” in diameter, 
evaporated on the specimen, were used for most 
measurements although it was observed that 
aluminum gold combinations responded as well. 
A schematic diagram of the specimen is given in 
Fig. 1 which also indicates the position of the 
probes for measurements. 

For the Hall measurements two B batteries 
were used to supply the sample current. It was 
found advisable to keep the current below one 
milliampere while fresh samples were being used, 
since the use of larger currents caused aging,’ 
or a decrease in the number of carriers. The 
e.m.f.’s which are determined initially with the 
potentiometer must be corrected for the contri- 
bution due to electrode misalignment and tem- 
perature gradients. Both of these effects are 
eliminated by taking two readings, one with the 
magnetic field in a fixed direction, and the other 
with the field reversed. From the theory of the 


1° The samples of cuprous oxide were kindly furnished by 
Mr. C. C. Hein of the Westinghouse Research Labora- 
tories, East Pittsburgh, Pennsylvania. A typical analysis of 
the copper and also the processing technique are given in 
reference 2. 
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* Fic. 2. Conductivity (curve A) and Hall constant 
(curve B) as a function of temperature for a fresh sample of 
CuO, quenched from 450°C. 


Hall effect the Hall coefficient is given by the 


equation 
R=End/TH, 


where d is the thickness of the specimen, H is 
the magnetic field in gauss, and J is the specimen 
current in amperes. 

Immediately following the Hall measurements 
point-probe conductivity determinations are 
made. The method consists of measuring the 
potential drop across a fixed portion of the 
specimen while a known current is flowing in the 
specimen. The advantage of the technique is 
that no current is drawn by the probes. The 
potential drop measured across the probes must 
be corrected for the contribution due to the 
thermal e.m.f. Around 100°C the latter represents 
about a 2 percent correction, and continues to 
increase till at 900°C the thermal e.m.f. becomes 
comparable with the true potential drop across 
the probes. This, of course, depends upon the 
current flowing through the specimen. For cur- 
rents of 10 milliamperes the correction at the 
higher temperatures is about 8 percent, while 
for one milliampere it approaches 20 percent. 


RESULTS AND INTERPRETATION 


In Fig. 2 are shown two curves (A) giving the 
variation with temperature of the logarithm of 
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the specific conductivity for a specimen of 
cuprous oxide quenched from 450°C. The upper 
curve was obtained with rising temperature and 
the lower curve with decreasing temperature, 
The various points were taken at time intervals 
of twenty minutes. It may be noted that on 
the lower curve corresponding values of the 
conductivity at a given temperature have de. 
creased considerably. This is indicative of a 
definite aging process in which an appreciable 
number of the current carriers is depleted. 
Attempts were made to obtain a reversible 
curve for a fresh sample by making the measure. 
ments rapidly. However, the rapidity with which 
aging set in above 60°C made this impossible, 
Also, the temperature to which the sample was 
raised during the run, and the speed of cooling 
may have influenced the amount of absorbed 
oxygen and removed all significance from the 
results. 

In curves (B) of Fig. 2 are given the Hall 
constant curves for the same specimen. These 
curves bear a one to one correspondence relative 
to the conductivity curves. The regions at which 
the conductivity shows a hump are those at 
which the Hall coefficient has a plateau, and 
the regions at which the conductivity rises 
rapidly the Hall coefficient decreases rapidly. 
It might be remarked that the curves of Fig. 2 
are representative of the initial behavior of 
specimens that were quenched from different 
temperatures. The only distinction is in the 
magnitude of the Hall constant and the specific 
conductivity. Since the Hall constant decreases 
and the conductivity increases with an increase 
in the number of carriers, we should expect that 
at a given temperature the latter will be large if 
the former is small, and vice versa. 

In Fig. 3 (curve A) are given the conductivity 
curves for a well-aged specimen run in air. The 
specimen had previously been quenched from 
250°C. The curves exhibit the same sort of 
characteristics as those in Fig. 2; however, we 
shall consider them in more detail. At the higher 
temperatures the conductivity increases in a 
manner described by a Boltzmann relationship, 
but below 300°C a marked deviation occurs. 
Part of the deviation can be explained by 
assuming that true equilibrium is not attained at 
the lower temperatures, that is, that oxygen is 


cul 














frozen in the lattice. It was possible to quench a 
given sample with a blast of cold air from a given 
temperature to room temperature and obtain 
conductivities which varied by a factor of 20. 


The higher the quench temperature the higher 


was the room temperature conductivity. This 
feature accounts for the many branches in the 
low-temperature region of the curve. These 
differences in the conductivity disappear around 
300°C for above this temperature all curves 
follow the same pattern. Presumably, the rate 
of diffusion of oxygen is sufficiently large in the 
neighborhood of 300°C to bring about equi- 
librium in the time allotted between measure- 
ments. The slope of the curves in the high 
temperature region turns out to be 0.78 ev and 
remains the same for all specimens regardless of 
past history. This slope agrees quite well with 
that obtained by Jusé and Kourtschatow‘ who 
found an activation energy for the high tempera- 
ture part of their curves of 0.72 ev. 

An experiment similar to that of Diinwald 
and Wagner® was performed in order to test for 
the pressure dependence of the resistance of 
cuprous oxide. A small slab of cuprous oxide 
was placed between two platinum contacts which 
were forced together by a spring, and the 
specimen was then placed in a furnace at 850°C. 
Resistance measurements were made for various 
oxygen pressures, and the resistance was found 
to vary markedly in the range of 10-* to 150 mm 
oxygen pressure. In Fig. 4 is shown the varia- 
tion of the logarithm of the resistance with the 
logarithm of the oxygen pressure for two speci- 
mens. The results of several samples give slopes 
ranging from 1/6.8 to 1/7.2 which agree quite 
well with Wagner’s value of 1/7. However, the 
departure from a 1/8 law as evidenced in all 
specimens appears to be a real property of 
cuprous oxide. Wagner suggests that an incom- 
plete dissociation of the ions and electrons might 
be the cause of the disagreement. 

Figure 3 (curves B) show the variation of the 
Hall constant with temperature for a well-aged 
sample. The surprising result obtained from 
these curves is that m, the number of carriers, 
cannot be represented by the relation 


n=ny exp (—€/kT), 


where mp and ¢ are constant. Data taken from 


CONDUCTIVITY AND HALL EFFECT 





117 


the curves show that ¢ varies from about 1.0 ev 
at 450°C to about 0.1 ev at room temperature. 
This agrees qualitatively with the Hall effect 
data of Waibel and Schottky;* however, no 
evidence was found for a change in sign of the 
Hall coefficient. 

Figure 5 shows the behavior of the mobility 
as a function of temperature. Again, qualitative 
agreement is obtained with the data of Waibel 
and Schottky. The mobility varies approximately 
as T-* from room temperature to 150°C, and it 
varies approximately as 7-’ from 150°C to 
450°C. This means that the conductivity, alone, 
cannot give accurate information concerning the 
variation of the number of carriers with tempera- 
ture. Furthermore, it seems difficult to explain 
this rapid variation of the mobility by assuming 
that the electrons are scattered only by lattice 
vibrations. One wonders whether some scattering 
or trapping process involving the number of 
impurity centers is not responsible for the 
observed behavior. It is difficult to reach any 
definite conclusions at present because little is 
known about the interaction of slow electrons 
with lattice vibrations. 

In conclusion the results may be summarized 
as follows: 
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Fic. 3. Conductivity (curves A) and variation of the Hall 
constant (curves B) of a well-aged sample of cuprous oxide 
as a function of temperature. The slope of the curve A at 
the higher temperatures is 0.78 ev. 
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Fic. 4. Resistance as a function of oxygen pressure for 
two specimens of cuprous oxide measured at 850°C. Slope 
of curve A is 1/6.8; curve B, 1/7.1. 


(1) There is no evidence that an intrinsic 
conductivity exists in cuprous oxide. 

(2) There is no evidence for a change in sign 
of the Hall coefficient. 

(3) The simple theory of semi-conduction 
needs considerable extension before it can fully 
explain the data. However, the following general 
interpretation is suggested. At low temperatures 
a frozen equilibrium exists in cuprous oxide, and 
the number of absorbed oxygen atoms is inde- 
pendent of the temperature. The extremely low 
activation energy observed in this range of 
temperature is to be associated with the forma- 
tion of clusters of centers, in the sense suggested 
by Miller? for the case of zinc oxide. The maxima 
and minima of the conductivity curves observed 
in the low temperature region result from a 
combination of a saturation of the clusters and 
the rapid decrease of the mobility with tempera- 
ture. Above 300°C an equilibrium exists between 
the absorbed oxygen and the external oxygen gas. 
As a result, the number of oxygen atoms ab- 
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Fic. 5. Mobility as a function of temperature for a well- 
aged specimen of cuprous oxide. The data are taken from 
Fig. 3. The low temperature part of the curve shows a T-* 
—. and changes to a 7~’ dependence beyond 


sorbed and the fraction of holes ionized depend 
exponentially upon the temperature. Apparently, 
we must believe that one or both of the activation 
energies, E and W (E is the energy to free an 
electron hole and W is the energy to place an 
oxygen atom in the lattice as an ion and produce 
two copper ion vacancies), are functions of the 
temperature, or the concentration of absorbed 
oxygen. 

Finally, the author wishes to express his 
sincere appreciation to Professor F. Seitz for 
his helpful advice and criticism, and to Drs. A. 
W. Lawson, R. J. Maurer, and S. Pasternack for 
many stimulating discussions on the interpreta- 
tion of the data. 
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A grounded tungsten wire (0.0078-mm diameter) was mounted coaxially with a cylinder 
which was made negative relative to the wire, and a difference of potential was applied between 
the cylinder and the wire. The procedure consisted of measuring the resistance of the tungsten 
wire first without an electric field and then with a field. With a vacuum of about 10-* mm Hg 
and with fields of about 10* v/cm, measurements showed that: (1) An increase in resistance 
always resulted with the application of the field; (2) the change was different for different 
temperatures of the filament and for different vacuum conditions; (3) no measurable change 
occurred unless the field exceeded a certain value; (4) a small leakage or ionization current 
was associated with the changes in resistance. The observed changes might have been caused 
by: (1) leakage current, (2) ionization current, (3) radial stress in the wire resulting from 
the high field, or (4) other reasons. Calculations showed that reasons (1) and (3) produce 
changes much smaller than those observed and that the results are explainable wholly by 





reason (2) if it is assumed that most of the ionization occurs near the surface of the wire. 





I. INTRODUCTION 


ERKINS' in 1921 and again in 1925? observed 
that gold and bismuth exhibited an increase 

in resistance in the presence of an electrostatic 
field. M. Pierucci,? M. Perucca,‘ and R. Deaglio® 


have observed changes in resistance of thin, 
metallic films deposited on quartz or glass upon 
application of a field. Their results indicate that 
there is a decrease in resistance upon applying a 
field. In work on the effect of high electrostatic 
fields upon the vaporization rate of metals A. G. 


Worthing® and G. B. Estabrook’ observed 
changes in resistance upon applying high fields. 
Worthing observed that there was a sudden 
decrease in the resistance of a tungsten wire 
upon applying a field and a sudden increase upon 
removing the field. Estabrook also observed this 
effect for the case of platinum and molybdenum, 
but the effect was opposite to that found by 
Worthing. 

The present research was undertaken to obtain 
quantitative results on the effect of high electro- 
static fields on the electrical resistance of small 
tungsten wires. 


1H. A. Perkins, Phys. Rev. 18, 131 (1921). 

?H. A. Perkins, Phys. Rev. 25, 584 (1925). 

3M. Pierucci, Nuevo Cimento 9, 33-42 (1932). 

*M. Perucca, Comptes rendus 198, 456 (1934). 

°R. Deaglio, Nuevo Cimento 11, 288 (1934). 

* A. G. Worthing, Phys. Rev. 17, 418 (1921). 

7G. B. Estabrook, Ph.D. Thesis, University of Pitts- 
burgh (1932). 


Il. APPARATUS 


The tungsten® used was in the form of a wire 
0.00078 cm in diameter. This value measured by 
a microscope checked well with that obtained 
with the value 0.18 mg per 200-mm length 
furnished by the manufacturer. Figure 1 shows 
a cross-sectional view of the apparatus used. 
Approximately a 10-cm length of this wire was 
mounted coaxially with respect to a cylinder. 
The two springs S made of 5-mil Kovar were 
used to keep the wire under slight tension. To 
facilitate centering, holes (0.29-mm diameter) 
were drilled along two diameters approximately 
at right angles to each other at the ends and 
also at the center of the cylinder. The wire was 
centered by means of the adjusting screws A at 
each end. The rings N, holding the adjusting 
screws, were grounded and served as guard 
rings. The ends of springs S were threaded 
through short glass capillaries encased in the 
plates at each end, and the tungsten wire was 
soldered onto the ends of the Kovar springs. 
Although solder does not adhere to tungsten, 
it made a good mechanical and electrical joint 
as long as the filament was not heated to a very 
high temperature. The cylinder and the rings 
were supported by glass rods and the entire 
assembly introduced into a glass tube connected 


’ Furnished by the General Electric Company, Lamp 
Development Division, Nela Park, Cleveland, Ohio. 
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To ION GAUGE 
Fic. 1. Cross-sectional view of apparatus. A, adjusting 
screws; B, glass beads; C, cylinder; G, glass capillary; 


K, Kovar tubes; N, guard rings; P, adjustment plate; 
R, glass rods; S, springs; L, brass plate (} normal size). 





to a mercury diffusion pump.® For ease in 
replacing filaments, the glass tube rested on a 
rubber gasket in a groove of the brass plate L. 
The electrical connections were brought to the 
outside through two Kovar tubes K soldered 
onto the brass plate. The ends of the tubes were 
sealed with glass beads B which made vacuum 
seals and insulated the wires from the plate. 
A No. 46 radio tube served as an ion pressure 
gauge. A Kohlrausch bridge was used for meas- 
uring resistances. Changes in resistance were 
measured by the galvanometer deflection from 
the balanced position, galvanometer deflections 
being previously calibrated in terms of resistance 
changes. All lead wires and instruments were 
shielded and grounded. The galvanometer of the 
Kohlrausch bridge circuit was introduced into 
the ground circuit in order to detect and measure 
any “‘leakage’”’ current. 

High voltages for obtaining the fields between 
the cylinder and wire were produced by an 
x-ray transformer and a Kenotron KR-3 rectifier. 
The negative side of the line was connected to 


® Built by R. Nestler at the Bartol Foundation. 
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the cylinder in the tube and the positive side 
was grounded at the transformer. The tungsten 
wire was also grounded. 


III. PROCEDURE AND RESULTS 


The tungsten wire was well centered along the 
axis of the cylinder, and the entire assembly 
(see Fig. 1) was then inserted in the tube 
connected to the vacuum pump. The tube was 
outgassed by heating it locally with a gas flame, 
and all measurements were taken with a vacuum 
of about 10-* mm Hg. Previous to making 


-measurements with the wire in the electrostatic 


field, the variation of resistance with current 
through the wire was measured with increasing 
and decreasing currents (Fig. 2). This was 
necessary in order to calculate changes in 
resistance resulting from a change in current in 
the circuit in which only the resistance of the 
tungsten wire would vary as occurred in the 
course of measurement. The relation is practi- 
cally linear. Measurements were made after the 
filament had been heated for about 10 minutes 
at the highest temperature always below incan- 
descence. The wire was not “well aged ;’’ never- 
theless, it behaved in a reproducible and con- 
sistent manner. Throughout, the current was 
measured to the nearest 0.0001 ma by means of 
a carefully adjusted potentiometer. The temper- 
ature of the wire was calculated from a tempera- 
ture-resistivity calibration. 

The deflections of the galvanometer in the 
Kohlrausch bridge from a balanced position 
were calibrated in terms of changes in resistance. 
For the most sensitive conditions, a deflection 
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Fic. 2. Variation of resistance with current for an 
unaged tungsten filament. Diameter 0.00078 cm, length 
10.38 cm. 
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of one division represented a change in resistance 
of 0.0015 ohm in 500 ohms. A deflection of 0.1 
division could be estimated so that it was pos- 
sible to detect a change of 3 parts in 10,000,000. 

The negative side of the high voltage source 
was then connected to the cylinder, and the 
tungsten wire was grounded at the galvanometer 
connection of the bridge. Because of the small- 
ness of the wire, a relatively small difference of 
potential gave rise to a high field strength at the 
surface of the wire. For example, a difference in 
potential of 10,000 volts between the wire and 
cylinder gave rise to a surface field strength of 
3.35X10® v/cm. Measurements of changes in 
resistance with applied fields were made on 
numerous tungsten wires. The results for one 
particular wire are shown in Fig. 3. It shows 
that there was always an increase in resistance 
with the application of an electric field provided 
the field exceeded a certain value below which 
no change could be detected. The magnitude of 
the change for a particular field strength de- 
pended on the temperature of the filament and 
upon vacuum conditions. The results on the 
other tungsten filaments were similar. The maxi- 
mum observed changes are shown in Table I. 

No change of resistance was observed until a 
surface field strength of 0.75X10° v/cm was 
reached for the filament at 20°C and the value 
at which any change was observed increased as 
the temperature of the filament increased. The 
field strengths at which a change was first 
observed for the various temperatures and 
vacuum conditions are shown in Table II. 

When the galvanometer of the bridge was 
moved to the ground circuit, it measured any 


TaBLE I, Maximum observed change in resistance of 
tungsten wire at various temperatures with field strength 
at surface of wire. 

















Field 
Temp. of wire Max. observed Percent a. 
> change, ohms change 106 v/cm 
Vacuum 3X 10-* mm Hg 
20 13.8 8.8 2.90 
280 aun 0.71 3.97 
0.37 0.072 4.77 





Vacuum 8X 10-7 mm Hg 


20 13.5 8.7 3.18 
280 2.5 0.77 4.77 
580 0.34 0.066 5.31 











EFFECT OF HIGH ELECTROSTATIC FIELDS 
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Fic. 3. Change in resistance with applied field of unaged 
tungsten. Curves A, C, E for pressure 3X10-* mm Hg. 
Curves B, D, F for pressure 8 X 10-7 mm Hg. Curves A and 
B—R=157.3 ohms, T =20°C. Curves C and D—R=323.9 
ohms, 7=280°C. Curves E and F—R=514.7 ohms, 
T =580°C. 


leakage current passing to ground. The measured 
results are shown in Fig. 4. The maximum 
deflection was never more than 2 divisions, at 
100 cm and represented a current of approxi- 
mately 3X10-* amp. The deflection appeared 
immediately as soon as the field was applied. 
As Fig. 4 shows, the current varied with the 
temperature of the filament and also with the 
vacuum much in the same way as did the 
changes in resistance. If the field was kept 
constant and the temperature was changed from 
one value to another, there was a change in the 
leakage current. There was an increase if the 
change was one from a high temperature to one 
of a low temperature, and a decrease if the 
reverse change was made, the increase and 
decrease always being the same. Also the final 
leakage current was the same as that found by 
suddenly applying the field. However, when the 
temperature of the wire was changed suddenly 
from a low temperature to a high, there was a 
sudden decrease in leakage current, but when 
the temperature was changed from a high to a 
low value, the current increased very slowly, 
possibly in consequence of the very slow adjust- 
ment to a new temperature equilibrium. 


IV. DISCUSSION 


Measurements with high electrostatic fields at 
the surface of the tungsten wire showed that 
the wire exhibited increases in resistance with 


applied fields. The magnitude of the changes 
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“ Loeb," and Druyvesteyn and Penning” that . 
» positive ions of sufficient energy striking the tl 
aw . . . 
: cylinder would liberate electrons from it. These 7 
°e electrons would then move toward the wire and ti 
produce more ions by collision. Regarding elec- sl 
& tron emission from the cathode by positive ions, w 
e Loeb quotes: ‘Concerning this process much e 
3 information is at hand, though none is satis- a 
3 factory owing to the impossibility of producing de 
5 or studying reproducible surfaces.” A steady ct 
4 5 é e 7 4 current will occur (neglecting ionization by out- gi 
POTENTIAL BETWEEN FILAMENT AND CYLINDER (IG*VOLTS) side sources) when the total number of positive fil 
Fic. 4. Variation of leakage current with potential. Curve ions produced by a Agee cemen a st 
designations same as for Fig. 3. one electron from the cylinder. The number of of 
extra electrons liberated from the cathode per at 
depended on the temperature of the wire and ion formed between the electrodes is not known di 
on the vacuum conditions. Associated with these for different electrodes at such high ratios of or 
changes was a leakage current which varied with _ field strength to pressure as used in the course of on 
the temperature of the wire and vacuum condi- these measurements, so that a quantitative tr: 
tions much in the same manner as did the changes discussion of this is not possible here. A complete re: 
in resistance. and detailed discussion of what occurs between fer 
The changes in resistance that resulted with concentric cylinders at low pressures and high II 
the application of high fields might have been field strengths is complicated by impurities and 
due to several reasons: general conditions in the tube.'® Quoting from ane 
1. Leakage current passing through the fila- Loeb: ‘‘In tubes at very low pressures, 0.01 mm, tio 
ment, thus heating it and increasing its resistance. the character of the gas becomes somewhat ot 
2. Ionization current in the space between the doubtful, for adsorbed gases coming from the & 
wire and the cylinder passing to the filament, walls and electrodes, whose nature is unknown, 514 
thus heating it and increasing the resistance. begin to play an important role. Below 10-' mm 
3. Radial stress in the wire resulting from the _ the gas present is completely indeterminate.”* 
high electrostatic field. An ionization current would therefore consist 
4. Others. of a migration of electrons or negative ions 
Simple calculations show that changes due to toward the wire and positive ions toward the 
reasons 1 and 3 are much too small to account cylinder. The energy of the electrons in passing 303 
for the measured changes in resistance. The fact through a potential difference in the cylinder 
that the leakage current measured varied with agre II. Field strengths at surface of tungsten wire at 
the temperature of the wire and with vacuum which a change in resistance was first observed. 
conditions suggests that it might be anionization ~~ 
current. That this is possible at pressures of *%¢” (Vac. 30070 ee Ra S18 on tl 
10-* mm Hg or better is the basis of all ionization 20 0.75 0.80 
pressure gauges. 280 1.3 1.7 
If an electron is generated in the gas in the 89 2.5 3.3 aw. 
cylinder by photons, y-rays, or cosmic rays, it 
| will move toward the wire and at some appro- 4 Leonard B. Loeb, Fundamental Processes of Electrical 
| priate distance will start to ionize by collision. 1939) 377. 490° (John Wiley and Sons, New York, 
The negative ions formed will move in toward 2M. J. Druyvesteyn and E. M. Penning, Rev. Mod. 
) the wire and the positive ions will move outward >> ity 11, p. 378. a. 
to the cylinder. It has been shown by Dempster,’° a - . J. Druyvesteyn and F. M. Penning, reference 12, p. . 
= 16 Leonard B. Loeb, reference 11, p. 451. = 


10 A. J. Dempster, Phys. Rev. 46, 728 (1934). 16 Leonard B. Loeb, reference 11, p. 480. meast 
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EFFECT OF HIGH 





would be transferred to the wire, heating it and 
therefore increasing its resistance. To calculate 
what changes in resistance such assumed ioniza- 
tion current would produce, Fig. 5 was drawn 
showing the variation of resistance of the wire 
with power (J°R) input. The slope of this curve 
was taken at the experimental values of re- 
sistance. The product of this slope and a power 
determined as dissipated because of ionization 
currents and transferred to the filament would 
give the expected change of resistance for the 
filament on the application of the high electro- 
static field. Since the ionization current consists 
of a stream of negative ions in one direction and 
an equal number of positive ions in the opposite 
direction, depending on their place of origin, 
one might expect something on the order of 
one-half of the calculated ionization power to be 
transferred to the wire. The changes that would 
result if exactly half of the power were trans- 
ferred to the wire are shown as AR, in Tables 
III and IV. 


TABLE III. Comparison of measured changes in resist- 


ance with computed changes assumed to be due to ioniza- 
tion power. Vacuum 3X 10-* mm Hg. 








R: Temp. E AR; AR» 
ohms <— 10° v/cm ohms ohms 


514.7 580 0.0068 
0.013 
0.026 
0.069 
0.14 
0.18 
0.28 
0.37 
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R:=Resistance of wire; E=Field strength at surface of wire; 


AR: = Measured change in resistance due to field; ARp =Calculated 
change in resistance from ionization power; (ARp/AR:) =Ratio of 
measured change to calculated change in resistance. 
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Fic. 5. Variation with power input of resistance of an 
unaged tungsten filament. Diameter—0.00078 cm, length— 
10.38 cm. 


Comparing the calculated with the observed 
changes in resistance, it is seen that the computed 
effects for the observed ionization currents 
account for more than the observed changes in 
resistance. The ratios of the calculated changes 
to the observed are shown also in Tables III 
and IV. 

The ratios of calculated changes to the ob- 
served changes for a particular temperature of 
the filament remain practically constant at 
different field strengths. The ratios, however, 
are different for different temperatures of the 
filament. These results are compatible if the 
assumption is made that there exists a region of 
greater density of gas molecules near the surface 
of the wire than anywhere else inside the 
cylinder and that the density of the gas near the 
surface of the wire is decreased when the temper- 
ature of the filament is increased. This could 
well be in the form of adsorbed gas on the 
surface of the tungsten. This assumption is 
reasonable because of the high field strength and 
the large gradient of the radial field existing 
there. Because of this greater density near the 
surface of the wire, a large part of the ionization 
by collision is assumed to occur in this region. 
Hence, the electrons or negative ions would not 
pass through the potential difference applied, 
but would pass through a smaller potential 
difference, and the energy transferred to the 
wire per second would be smaller. 

By taking the requisite potential difference to 
give the measured changes in resistance, the 
distance from the surface of the wire giving this 
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TABLE IV. Comparison of measured changes in resistance 
with computed changes assumed to be due to ionization 
power. Vacuum 8 X 10-7 mm Hg. 











R:z Temp. E AR: ARp ARp 
ohms ~ 106 v/cm ohms ohms AR; 
513.4 580 3.18 0.011 0.046 4.2 

3.44 0.034 0.10 2.9 
3.97 0.093 0.17 1.8 
4.24 0.130 0.25 1.9 
4.50 0.162 0.26 1.6 
4.77 0.246 0.49 2.0 
5.31 0.339 0.69 2.0 
322.6 280 1.99 0.018 0 0 
2.38 0.067 0.12 1.8 
2.65 0.14 0.26 1.9 
3.18 0.35 0.47 1.3 
3.71 0.71 0.73 1.0 
4.24 1.40 1.5 1.1 
4.77 2.50 2.8 1.1 
156.1 20 0.93 0.14 0 0 
1.06 0.47 0.4 1.4 
1.59 1.94 4.8 2.5 
2.12 4.48 11.4 2.6 
2.65 9.13 23.7 2.6 
2.91 13.5 31.2 2.3 








Rz=Resistance of wire; E=Field strength at surface of wire; 
AR; = Measured change in resistance due to field; ARp =Calculated 
change in resistance from ionization power; (ARp/ARs) =Ratio of 
measured change to calculated change in resistance. 


potential difference can now be calculated. The 
potential at a distance r from a grounded 
cylinder of radius a having a charge q per unit 
length inside another cylinder of radius D is 
given by the relation 


V,= —2q In (r/a). (1) 


If the potential is to be a fraction of the applied 
potential, say ,Vs, then 


Ve —2qIn (b/a) am 1 _In (b/a) 





—= =-=——,_ (2) 
sVp —2qIn(r/a) f In (r/a) 
from which is obtained 
r=a(b/a)’. (3) 


Substituting the values for a and 5 one obtains 


8.2 \/ 
r=0.0039 mm (——) ; 
0.0039 


The average value for f when the wire was at 
20°C, pressure 3X10-* mm Hg, is (1/2.8) and 
leads to 0.06 mm for the value of r. This distance 
is 7.3X10-* of the radius of the cylinder or a 
distance 15 times as great as the radius of the 


wire. This value of r would indicate that most of 
the ionization occurs in a volume of 1.110- 
mm per unit length of the wire or 5.4X 10-5 of 
the total volume. The value of r when the 
filament was operated at 280°C, vacuum 3X 10-8 
mm Hg, was 1.4 mm or 0.17 of the radius of the 
cylinder. The ionization occurred inside a volume 
6.2 mm? per unit length of wire, or 0.029 of the 
total volume. This value for r at 280°C was 
larger than the value at 20°C. This is in agree- 
ment with the assumption made, namely, that 
as the temperature of the filament is increased, 
the density of the molecular cloud around the 
filament decreases. The ionization would be 
expected to occur, therefore, at the higher 
temperature throughout a larger volume and to 
be smaller for a given voltage applied between 
wire and cylinder. The value for r at 580°C 
should be larger than at 280°C, and the ionization 
should be less for any given voltage applied. 
However, the ionization current at the highest 
temperature was found, as expected, to be 
smaller than at the lower temperatures, but the 
value of r, contrary to expectations, was smaller 
than the value for r at 280°C. To check with 
what was expected, the leakage current should 
have been very much smaller. A larger ionization 
current than expected might have resulted from 
several reasons: (1) The filament might have 
liberated adsorbed gases at the higher tempera- 
ture; (2) A noticeably greater part of the 
ionization current at the high temperatures was 
leakage surface current because of the higher 
voltages at which the data for the high tempera- 
tures were taken (see Fig. 3 and Tables III 
and IV). 

The values of r calculated for various condi- 
tions are shown in Table V. As has been men- 
tioned, the results depend considerably on the 
impurities present and the nature of the cathode 
surfaces. 


TABLE V. Values of radii inside which most of the 
ionization occurs. 











Rz Temp. r r 
ohms °C mm mm 
(Press. 3X10-* (Press. 8X107 
mm Hg) mm Hg) 
514 580 0.04 0.22 
323 280 1.4 2.3 
156 20 0.06 0.08 
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There are reasons to suspect!” that the 
ionization may not occur in such large volumes 
as indicated in Table V. The results can be 
explained then if it is assumed that only a 
fraction of the measured ionization current 
passes through the potential difference applied 
between cylinder and filament. The mechanism 
can be viewed in this fashion. It is assumed that 
there exists adsorbed gases and/or adsorbed 
tungsten oxide at the surface of the tungsten 
filament. Suppose that by some means or other 
an electron is liberated at the inner surface of 
the cylinder, the cathode. The electron will be 
accelerated towards the filament and will ionize 
only when it strikes the adsorbed layer. There 
equal number of electrons and positive ions will 
be produced; the electrons passing through a 
very small potential difference to the filament 
while the positive ions are accelerated towards 
the cylinder by passing through the entire 
potential difference. The positive ions arriving 
at the cylinder will liberate more electrons and 
the same process repeated again. A steady 
current will result when the number of positive 
ions produced by one electron at the anode or 
filament will produce one electron at the cathode 
or cylinder. 

Let the ratio AR,/AR, as previously defined 
be designated by “f.” To account for the 
measured change in resistance, AR,, the number 
of electrons passing through the applied potential 
difference per second would have to be a fraction 
“f” of the measured number. One electron, 
then, on the average would produce (1—f)/f 
positive ions at the filament and (1—/f)/f positive 
ions would in turn produce on the average one 
electron at the cylinder wall. The number of 
electrons per positive ion liberated at the cathode 


W. P. Reid, Phys. Rev. 63, 359 (1943). 
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is known as an “alternative second Townsend 
Coefficient.’ The value for this coefficient is 
given by the relation 


y=f/1—-f. 


Using the average values for 1/f=AR,/AR, 
found in Tables III and IV we find the values 
for y at a pressure of 3X10-* mm Hg to be 0.40, 
3.3, and 0.55 for temperatures 580°C, 280°C, and 
20°C, respectively, and at a pressure of 8X10-* 
mm Hg, they are 1.1, 2.5, 0.66 at the same 
temperatures. These values are averages at a 
field strength to pressure ratio, E/P, of the 
order of 10" v/(cm mm Hg). Other experimental 
values for y are not available at such high ratios 
of E/P as used here so that no comparisons can 
be made here. It must be remembered that these 
coefficients apply only to the surfaces used for 
this investigation. 

Variations in resistance with applied field 
might havé resulted because of other reasons 
such as a change in emissivity or a structural 
change in the wire itself. These changes are 
entirely possible, but an explanation of the 
results wholly on these changes would not 
account for the effects of the ionization currents, 
and it is seen that these can account for the 
observed changes. To detect changes due to 
other assumed reasons, one would have to reduce 
ionization currents to negligible amounts, and 
this can be done only in a vacuum of 10-® mm Hg 
or better. 
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18 Leonard B. Loeb, reference 11, p. 687. 
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Equivalent Circuits to Represent the Electro- 
magnetic Field Equations 
GABRIEL KRON 
General Electric Company, Schenectady, New York 
August 21, 1943 


T is known that for certain electromagnetic waves some 
of the characteristics of the wave propagation may be 
correctly described in terms of distributed inductances, 
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capacitances, and resistances. This is commonly done, for 
example, in conventional transmission line analysis. Here 
it is the integrated effects (voltage and current) that are 
represented by the circuit and not the fields. For static 
fields, another example, the field distribution itself may be 
related to current or voltage distribution in a network of 
similar circuit elements. The general relation between 
electromagnetic waves and distributed constant networks 
has apparently not been investigated. The extent to which 
such networks may represent characteristics of various 
waves or indeed space in general as a medium for the 
propagation of any wave distribution, is apparently not 
known and is of considerable importance. 
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Fic. 1. Transient equivalent circuit of the field equations of Maxwell for charge-free space. 


Max 
netw 
whicl 
dielec 
magn 
of th 
repre: 
ment: 


LETTERS TO 


z: % hs au®au Kp 
h, ou c 


a ated Hp 


4u 


na aces ee 


2)€,,-MODE (H,:€, £5) 


EDITOR 


hghs autau® o 


h, au c 


hyhp au'au® o 


hs au> 


h 
Sanb 
z hphsaut au 


: HP 
h 4u c 
Tihg Bula jp 


Le 


4) Hyp MODE (E,:H,: Hs) 





Fic. 2. Two sets of two-dimensional networks. 


This has suggested that attempts be made to represent 
Maxwell's equations directly by equivalent circuits. A 
network which does this is shown in Fig. 1. In this network, 
which represents a charge-free space of permeability y, 
dielectric constant ¢, and conductivity ¢, the electric and 
magnetic field components along the three directions 1, 2, 3 
of the orthogonal curvilinear coordinates 1%, 2, “3, are 
represented by voltages and currents in impedance ele- 
ments as indicated. (i, hz, h; are the metrical coefficients 


and p=0/dt.) The values of these elements depend only 
upon yu, ¢, ¢, and the size of the infinitesimal cube of space 
being replaced. Thus the network represents space rather 
than any particular wave mode or frequency; the latter, 
as in reality, being determined by the boundary conditions. 

For two-dimensional cases (8/8u,=0), the three-dimen- 
sional network separates into two independent networks 
lying in one plane (Fig. 2). All networks shown are merely 
examples. Other equivalent networks have been obtained 








which together with the derivations are expected to be the 
subject of a later paper. 





On the Branching Ratio of Na”? 


Hans WELTIN 
University of California, Berkeley, California 
August 2, 1943 


ASLETT?! has shown that Na* may be formed by the 
bombardment of Mg with deuterons. 
Mg*+ D—Na®-+ He‘. (1) 
Na emits positrons, and on the basis of the Fermi 8-decay 
theory it can also decay by K capture. Lamb* has extended 
the calculations of Yukawa and Sakata to Uhlenbeck- 
Konopinski coupling. His calculations of the ratios of the 
number of nuclei excited to the number of positrons 
emitted are quoted here: 


Fermi Uhlenbeck-Konopinski 
Ai=0 1.2 4.7 
Ai=1 2.1 30.4 


Lamb suggested that the theory might be checked by 
measuring the amount of He produced in (1) and com- 
paring it with the number of positrons emitted. Alvarez* 
indicated later that the theory might be checked equally 
well by comparing the amount of Ne with the number of 
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positrons emitted according to 
Na*®—+>Ne®+e?*. (2) 
Na®+e-—Ne”. (3) 


The latter method was used to obtain this ratio, A 
procedure for micro-gas analysis very similar to that 
developed by Paneth and Peters‘ made it possible to 
observe the pure Ne spectrum in a fine capillary and then 
to measure the gas. The number of positrons emitted was 
determined by comparing a known fraction of the total 
sample with a Segré standard of uranium oxide, U;0s, by 
means of Geiger point counter and scale of 64. The average 
amount of Ne measured for a decay period of 58 days was 
found to be 3.8X10-* cm’? N.P.T. The amount of Ne 
calculated from the number of positrons emitted was 
1.2<10-* cm* N.P.T. Hence, the ratio of the number of 


- nuclei excited to the number of positrons emitted js 


3.8X 10-*/1.210-*, or 3.2. A comparison between this 
result and Lamb's calculations supports the Fermi theory, 
provided that the change in nuclear spin is Ai=1, as 
Lamb has suggested. In order to distinguish convincingly 
between Ai=0 and Ai=1, a stronger Na™ sample than the 
one available for these experiments would have to be 
employed. 

I am indebted to Dr. Emilio Segré for his very helpful 
discussions during the entire course of this research. 

IL. Ji. Laslett, Phys. Rev. 50, 388(A) (1936). 

2 W. E. Lamb, Jr., Phys. Rev. 50, 388 mo 


IL. W. Alvarez, Phys. Rev. 54, 486 (19 
‘F. Paneth and K. Peters, Zeits. f. physik. Chemie 134, 353 (1928), 
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